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L. S. PONTRYAGIN'S MAXIMUM PRINCIPLE IN OPTIMAL 
SYSTEM THEORY - II* 


L. I, Rozonoér 


(Moscow) 


Questions are treated which are related to the proofs and applications of 
L, S, Pontryagin's maximum principle in optimal system theory, The work also 
contains certain new results, The first portion of the work deals with the opti- 
mization problem for the case of a trajectory with a free right end. In thesecond 
portion, the maximum principle is formulated for a more general type of bound- 
ary conditions, The third part establishes the connection between the method of 
dynamic programming and the maximum principle, provides a method of solving 
the optimization problem for linear discrete systems, and also presents a number 
of considerations in the use of the maximum principle for solving a definite class 
of problems, related to the theory of dynamic accuracy of control systems, 


3. The General Problem 





We consider the problem, formulated in Section 1 (Part I) of optimization,given arbitrary boundary con- 
ditions for the trajectory's right end. The principle special feature of this problem, as compared with the problem 
considered in Section 2, is this: that inthe present case the requirement that the trajectory's right end lie in set G, 
imposes, by its very nature, an additional limitation on those controls from among which the optimal one is to 
be chosen, Indeed, the optimal control must here be chosen only from the class of those controls which guarantee 
that the system will fall in the set G. Naturally, this circumstance gives rise to additional difficulties, It turns 
out, however, that the action of the maximum (minimum) principle in this case is completely analogous to the 
principle formulated in Section 2, and the special feature of the problem under consideration here, as cited above, 
is expressed in the greater complexity of the boundary conditions for the corresponding differential equations. We 
first give the fundamental result (theorem 3) which asserts the applicability of the maximum principle, and we 
then give a method of determining the initial conditions, 


In the sequel, we shall always assume that the set G in which the end of the trajectory x(t) must lie is closed 
and convex* . 


In the proof of the fundamental theorem (Cf, the Appendix), it is also assumed that the set G has at least 
one interior point in phase space X, i.e., a point in a sufficiently small neighborhood of which there are contained 
no points which do not lie in G, It is easily proven that, with this, there are interior points in any arbitrarily small 
neighborhood of any point of G. 


a. The Maximum Principle 





We now turn to the exposition of the basic result, Let the initial system of equations have the form: 


ay = fey, ++ 25 ni Uy, ++ Uri t),zi(T,)= 29 (i=1,...,R), 


* Part I of this work was published in Automation and Remote Control (USSR) 20, 10 (1959). 
** A set G is said to be convex if the straight line-segment joining any two points of G, lies completely in the setG, 
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We consider the set of vector functions p(t) which satisfy the relationships 


. of (z, » t) , 
rn=— Dre (i=1,...,m), (2) 
1 


Equation (2) (in contradistinction to the problem of Section 2) does not define the functions pj(t), since the 
boundary conditions for p(t) are still not fixed. Therefore, even if the control u(t) and the trajectory x"(t) cor- 
responding to u(t) are given, the vector p(t) is still not uniquely defined by (2). 


n 


As before, we introduce the function H (z, p,u,t) = >) Pili (Xy, . + +» Xpi Uy, . ~ «» Ur, t), by means of which 
Eq. (1) and Eq, (2) are written in the form 


. . 0H j 
A= 3) B= aT) =H Conky... 5m), (3) 


Let u(t) be some admissible control, x"(t) the trajectory corresponding to it, and p(t) some one of the vector 
functions which satisfy relationships (2). By substituting the values of the variables x(t), pj(t) in the function H 
(x,'p, u, t), we obtain the quantity 


K (t,u,,...,Ur) =A [x(t), p(t), u, tl, 


which, for each fixed moment of time t, is a function of the point u = (uy, . . ., up) which lies in the set U of the 
space R, We say that control u(t) satisfies the maximum (minimum) condition with respect to the vector p(t) if, 
at any fixed moment of time t (T, =< t < T), the function K(t, u) attains an absolute maximum (minimum) on the 
set U for given variables, equal to the values of the control at the same moment of time, i.e., for u, = u;,(t) 
ees ® 


As was shown in Section 2, Part I, the geometric meaning of the maximum (minimum) condition amounts 
to this: that the control which satisfies this condition tends, at each moment of time, to “accelerate” maximally 
the representative point in the direction defined by the vector p(t). 


We again emphasize that vector p(t) is, om the one hand, not arbitrary and, on the other hand, is still not 
fixed (in contradistinction to the problem of Section 2, Part I). Therefore, if some permissible control u(t) is 
given, it is not immediately clear whether one may choose such a vector p(t) that, relative to this p(t), the con- 
trol u(t) would satisfy the maximum (minimum) conditions, It turns out that an optimal control possesses precisely 
this property. This fact also constitutes the basic content of L. S, Pontryagin’s maximum principle. 


In order to formulate exactly the theorem which expresses the maximum principle, we make one more re- 
n 


mark, We consider the function (x) = >) cixi of points x = (x;, . . . Xp) on the set G, The set of points 
1 
x* (x*;, ..., Xn), at which the function ~(x) assumes a value which is a minimum with respect to all other points 


n n 
of G, we denote by G*. If the set G* exists, then the following inequality always holds: ¢° = 7; oti < yc 
1 1 


n 
where x* is an arbitrary point for G* and x is an arbitrary point of G), The functional S = dens (7') can not, 


obviously, assume values less than ¥*. Then, an arbitrary control which takes the system over into set G* is min- 
optimal, since for such a control the functional assumes the least of all its possible values, equal to ~*, There- 
fore, if the time, T—T), during which the system passes in the necessary optimal way to set G, is given in such a 
way that the set G* turns out to be admissible, the problem of minimizing the functional S then reduces to the 
problem of determining some control which takes the system from the fixed point x® to the fixed set G*. We call 
this a degenerate problem, It is not considered in the present paper. 
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n 
Analogously, we define the. case when one deals with a maximum of the function § = >} cia4 (T). as a 
1 
degenerate. 
The nondegeneracy of the problem may be perceived in practice by comparing the given time T—Ty with 
the minimal time T* necessary for translation of the system from point x° to the set G*. 


We note further that the problem considered in Section 2, Part I can not be degenerate for any given time T. 
n 


Indeed, in the given case, the function > c;z; is not limited to the set G (occupying the entire space) and, con- 
1 


sequently, assumes an extremal value at no point. The matter is analogous in those cases when the set G is open* 
or is not bounded, 


We now formulate the fundamental theorem, 
Theorem 3, If, in a nondegenerate problem, the control u(t) is min-optimal (max-optimal) with respect to 


n 


S = Dy c,x;(T), then there exists a vector function p"(t) such that u(t) satisfies the maximum (minimum) con- 
1 


dition with respect to p(t). 


Theorem 3 is valid, certainly, in that special case which was considered in Section 2, Part 1. Theorem 1 
asserts that, in the case of the problem with a free trajectory right end, the vector p“(t), the existence of which 
is guaranteed by Theorem 3, is subject to the boundary condition p"(T) =—c, It is now necessary to formulate 
an analogous condition for the vector p4(t)** in the general problem being considered. 


Theorem 3 and the boundary conditions formulated below can be obtained on the basis of the results of work 
[1], in which the maximum principle with the corresponding transversality conditions was proven under more 
general assumptions, 


b. The Boundary Conditions 





Let the point x! = (xj, . . ., x4) be the end of an — maT x(t). That portion of set G for which 


Sea > cyl, we denote by G’; the portion for cniceThee > » a we denote by Gt, The sets G and Gt 
1 1 n 
siden, are closed and convex and have the hyperplane C, given by the equation > (4, — z}) == 0, for 
1 
their common boundary, Point x’, lying in both G~ and G*, is a boundary point of both sets. For definiteness, we 
n 


shall now consider the case when one seeks a minimum of the functional S = deus (T). 


If the control u(t) is min-optimal, then there do not exist controls from the admissible class which would, 
during the time T— Tp, translate the system into set G_ since, for interior points of set G’, the value of the quan- 
n 


tity >) cr is less than for the point x, 


1 


With this, it turns out (Cf,, the Appendix) that the vector p(T), at the 


moment of time t = T, must be orthogonal to some hyperplane A, which is a bracket*** for set G’ at the point 
x' and, moreover, must be directed to that side of hyperplane A where set G™ lies, Therefore, if the equation of 


*A set G is said to be open if none of the boundary points of set G lies in the set, 

** In the sequel, the superscript u will be —- from p(t). 

*** A hyperplane A, containing ‘the point x’, is called a bracket to the convex set G™ at the point x’ if set G 
lies completely to “one side” of A. It has been proven that a bracket hyperplane may be passed through each 
boundary point of a convex set, 
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hyperplane A has the form >a (x; — 2!) = 0, and the signs of the coefficients aj are so chosen that set G~ lies 
1 


n 
on the side of the hyperplane where >} a (x; — x!) 0, then 
1 


pi(T) = —a, (i=1,...,A) (4) 


or p(T) = —a, if we introduce the vector a = (a),..., a.,). Equation (4) generalizes the corresponding conditions 
for the problem of Section 2, Part I, In fact, two cases of the placement of point x! = x(T) are possible: either 
inside G or on the boundary of G, In the first case (Fig. 1) there may be passed the unique hyperplane A, bracket- 
ing G™ and coinciding with hyperplane C, through point x', Any other hyperplane A through point x’ will neces- 
sarily “intersect” the set G’ and, consequently, will not be a bracket, Therefore, if the end of the trajectory lies 
in the interior of set G, the aj = c; and pj(T) = —cj. In the problem considered in Section 2, Part I, this case was 
always realized, since the set G occupied the entire space and in general, did not have boundary points*. 











Fig. 1 Fig. 2 


If the second case holds, and point x! is a boundary point of set G, then the bracket hyperplane to G" is not 
uniquely defined (Fig. 2). We shall assume that, at point x’, there exists a unique hyperplane B which is a bracket 
to set G, The sign of the coefficients b; of the hyperplane B are so chosen that set G lies on that side of B where 


nn 
yo (1; — z}) <0. It is then easily seen that no hyperplane A which is a bracket to set G~can lie in that region 
1 n 


n 
of space where the inequalities > (4% — z}) < 0 and yb: (x; — z}) < 0. hold simultaneously, By now taking 
1 1 


into account that G~ lies on that side of A where ja (x; — x!) < 0, we obtain the following expression for the 


1 
coefficients a,: 


a; = he; + pd; (i= 1,..., 0), (5) 


where , 20 and p 20 are certain nonnegative numbers which do not reduce to zero simultaneously, For A= 0, 
hyperplane A coincides with B and, for uy = 0, with C, IfX> ) (or p > 0) we can, without loss of generality, set 
A = 1 (or p = 1). Of both constants are positive, then either of them can be set equal to unity. 


This follows from the fact that, to define hyperplane A, the coefficients aj need only be given to within an 
arbitrary positive constant (this constant is required to be positive in order that there be no change in the sense 


of the inequality Da(ai—2)< 0). Therefore, essentially,only one constant enters into conditions (5). 
1 


* This case is also always realized in the problem in which the set G is open, For these problems, all the results 
of Section 2, Part I are completely valid. 
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Let the set G be given by the inequality F(x;,...,Xp) = 0. Then the coordinates of the boundary points of 
set G (including the coordinates of point x’) must satisfy the condition 


F (x, 00+52%n) = 0. (6) 


For each boundary point we can find coefficients b;(x;, . . ., X,) which define the hyperplane B, which is 
a bracket to G at the point x, In particular, if the function F(xy,.. .,Xn) is differentiable, then the bracket is 
the hyperplane which is tangent to surface (6) and therefore 


OF (21, .. +4 %q) 


os, (i= 4i,...,”), 


bi (a1, ee +9 Zn) = 





By uniting conditions (4)-(6), we convince ourselves that the coordinates of vectors p(t) and x(t) at time 
t = T must satisfy the set of relationships 


pa(t) = — dey — phi (2, (T),...,2n(T)] (i=1,...,M), 


F (x, (T),...,%n(T)] = 0. 


(7) 


We thus obtain n + 1 equations, (7), in which, as was stated above, one additional constant (Aor #) must 
enter, Therefore, (7) gives essentially n independent relationships. Together with the n relationships 


x (T,) = 29 (im i,...,n) (8) 


we obtain 2n boundary conditions for the solution of system (3) of differential equations, also of order 2n, The 
control u(t) is eliminated from system (3) by the same method as in the case of the problem of Section 2, since 
the maximum condition allows one, generally speaking, to express u(t) at each moment of time in terms of x(t) 
and p(t). 
n 
In the case when functional 5 = >} cjx; (7), attains a maximum, boundary conditions (7) and (8) are re- 
1 
tained if the inequality giving the set G is written in the form F(x;,...,%,) =0. To eliminate u(t) from system 
(3), it is necessary to use the minimum condition, 


c. Determination of the Boundary Conditions in Certain Particular Cases 





If some coordinate x, does not enter into the inequality which defines the set G, then the corresponding 
coefficient bg in (7) is equal to zero (G is a “cylindrical"body with its generator parallel to the x, axis and, con- 
sequently, the bracketing hyperplane is also parallel to the xz axis), 


In many cases of practical importance the set G lies on some linear manifold of phase space X and, con- 
sequently has no interior points in X (since in any arbitrarily small neighborhood of any point of a linear manifold 
one can find a point which does not lie in the manifold), In particular, this circumstance occurs if some of the 
coordinates are exactly fixed at time t = T but the set G is given by the system of equations 


z,(T) =z (s=1,...,93 <n), (9) 


We consider this case, Since, in practice, a sufficiently small violation of conditions (9) is always admitted, 
one could require that instead of (9), the following inequality hold 


> (x, aad xi)? < e*, (9") 


1 


which, for small ¢€ , guarantees that the violation of conditions (9) is sufficiently small and that | X,—Xxgl = €. 
Inequality (9°) gives a convex and closed domain with interior points in X. However, it turns out to be possible 
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to write simpler boundary conditions than those which follow from (9°). Conditions (9) can be directly considered 
as the boundary conditions for system (3) by replacing by them the conditions from (7) for the corresponding q. 
Since, with this, the equations defining G do not depend on the remaining n—q coordinates, then b 
= by = 0 so that there remains in the result the following set of 2n boundary conditions: 


q+i +e. 


x; (T,) = 2%, imi,...,M, 
x, (T) = x1, e=i,...,9, (10) 
Pm(T) =—Cm,™ =9+14,...,2. 


We now consider the frequently arising problem when it is required to determine the control which gives 
T 


an extremum to the integral \ @ (xz, u) dt, where, for t = Ty and t = T, the values of the coordinates must be 


fixed: - 


%(T,)= 29, x(T)=2) (i=1,...,2), 


t 


The new variable xp 4 ;(t) =| @ (x,u) dt does not enter either into definition of set G or in the right 
Te 


members of the original system of differential equations, By virtue of this latter fact, py 4 ,(t) = 0. Since cy = 
=Ce=...= Cp = 0, Cy 4 3 and Dp 4 3 = 0, then pp y 3(t) = pn 4+4(T) =—1. The system of differential equations 
and the boundary conditions for them take the following form: 


. . OH od 
an ae Pia? or: ae H =)} pfs (x, u, t) — (2, u), 
mt : 1 (11) 


x; (T,) = 2°, 2; (T) = 2 (i=1,...,m), 


The equation X, , 1 = (x, u) is excluded, since it gives no additional condition, 


d. Optimization Problems in Which the Time of Motion is not Fixed Beforehand 





Until now we have considered problems where the time required to implement the optimal process in the 
system was fixed. On the basis of the results obtained, we can now easily formulate the solution for the problems 
in which the amount of time taken by the system to transfer from point x° to the set G is given beforehand, but 
must be so chosen that the functional takes a least (greatest) value in comparison with the values which would be 
obtained for all other possible amounts of this time. 


We note first of all that if the problem is degenerate for at least some values of T, then these values of T 
also provide a minimal (maximal) value of the functional in comparison with all other values, such that the 
optimization problem with nonfixed time reduces to the degenerate problem, We shall therefore assume that 
the problem, with arbitrary fixed values of T, is nondegenerate. 


It is obvious that all the results given in this section must remain in force for problems with “free” times, 
Indeed, if the time sought exists and is equal to T, one can consider the problem with fixed time T = T, the 
solution of which would be the solution of the initial problem, However, since the magnitude of T is unknown, 
one additional condition which would define T is still required, It turns out that such a condition is the require- 
ment in accordance with which the function H(x,p,u,t) must equal zero at the end of the motion: 


Ypi(T) fila (7), u(T), 7] =0. 
1 


If u(T), in correspondence with the maximum (minimum) condition, is expressed in terms of x(T) and p(T) 
and substituted in (12), the equation obtained gives the additional boundary condition which, in conjunction with 
boundary conditions (7), allows one to find the solution of system (3) and the optimal time T. 


1410 











Fig. 3 Fig. 4 


Condition (12) has a simple geometric meaning. Indeed, the function H is the scalar product of the vector 
p(t) by the vector x(t) for the velocity of the representative point, The scalar product equalling zero means that 
at the end of the motion the velocity vector must be orthogonal to p(T) and, consequently, must (for example, 
in the min-optimal case) lie in hyperplane A, which is a bracket for set G’(Fig. 3), This requirement becomes 
particularly clear in the case when the end of the optimal trajectory lies in the interlorofG (Fig. 4). If the 
velocity vector were directed to the interior of G™ then, by increasing the time T, one could so arrange matters 
that the end of the trajectory would lie inside of G~ which ts impossible if the control u(t) is min-optimal, [If 
the vector X(T) were directed from G™ “outward", this would mean that for a smaller value of T the represent- 
ative point would be still found interior to G’ which, again is impossible. 


e. Problems of Maximal Speed of Operation 





We can now turn to the problem of the minimum time of transition from the given point x® to the set G, 
For its solution, it suffices to add the equation %, 4 , = x (where x is an arbitrary positive constant) and to solve 
the problem of minimizing xp 4 ;(T) for a time not fixed beforehand, It is easily seen that with this we have 


n 
in boundary condition (7), ¢y =... = Cp = 0, Pp 4 (T) =— 1 and, in (12)H = >) p (T) fax (7),u(T), Te = 
= 0. Moreover, Pp 4;(t) = 0 and pp 4;(t)=— 1. Since is an arbitrary positive constant, we can introduce the 
n 
function H = >} pifi (x, u, t),, equal to H+ & and, finally, we can write the system of differential equations 
1 


in the form: 


r= p= —Z (Hts...) 


‘ (13) 


and its boundary conditions in the form 
x (Tq) = 29, pi(T') = —pbi[z (T),---,2n(T)) = 4s --- om) 


< 14 
F(x, (T),.<-,%n(T)=0, Dipe(T)folz(T), 2(T), 71>. “ 


If the set G is a fixed point x! then, by substituting (x, u) = a > 0 in formula (11), we easily obtain the 
following boundary conditions for system (13): 


a(T,) = 29, 4 (T)=a2i(i=1,...,n), >) p(T) fstx(T), u(T), T1> 0. (15) 


f. Linear Systems 








The maximum principle, as was already mentioned for the case of the particular problem of Section 2, 
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gives jn general,only necessary conditions for optimality, However, for linear systems, described by equations 
of the form 


a 
2, = Dix (t) te + 9; (Wy, . 5 Ur) (i=1,...,m), (16) 
1 


sufficient optimality conditions can be established which show that the maximum principle, as applied to such 
systems, is “almost® sufficient. We formulate the corresponding theorem, 


ry before, we will weaete by bi(x}, . « «» X4) the coefficients of the hyperplane bracketing set G at the point 
xis (xi, oo xi) if point x' is a boundary point of G. To solve the problem of a minimum of functional S, the 
* 
signs of the coefficients by are so chosen that, for each point x(x;, . . .. X,) ¢ G, the inequality > 5: eae oi 
1 
(xj- x) 0, holds, and for solving the problem of a maximum of the functional, the inequality holds with re- 
versed sense. We do not me now that there is a unique bracketing hyperplane to G at the point x4, so that 
we denote by by(x}. « « «» XA) the coefficients of any bracketing hyperplane, 


If point x 14s interior to G, we should. set by(xt, . . os xi) = 0, Then, in the case of a minimum of the func- 
tional, the following inequality always holds for x' ¢G and x ¢G, >) bi (zi, ..., 21) (ai — 2}) < O(and the in- 
1 


equality with reversed sense for the problem of the maximum of the functional), 


Now let x(t) be the trajectory corresponding to control u(T), and let x(T) = xi, 


holds* , 


Then the following theorem 


Theorem 4, If in system (16), control u(t) satisfies the maximum (minimum) condition with respect to 
vector p(t), whose coordinates at time t = T take the value, 


pi(T) = — de; — pb; (xj, ..., zi) (i =14,...,2), 


nn 
where p is nonnegative and) is positive, then the control is min-optimal (max-optimal) for S = > ca; (T). 
1 


The conditions of Theorem 4 are very “close" to the necessary conditions for optimality expressed by 
Theorem 3 and formulas (4) and (5) of this section, The sole difference amounts to this that, unlike formulas 
(4), it is required in the conditions for the theorem just formulated that the constant A be strictly positive. 


The proof of Theorem 4 is given in Appendix II, Section 2°, We remark that the proof rests solely on 
Theorem 2,on the sufficiency of the maximum principle in linear problems with free-trajectory right ends and, 
in particular, does not assume that the set G must have interior points. 


g. Continuity and Piecewise Differentiability of the Function H 





We now formulate a theorem which expresses a very characteristic property of optimal processes, Let the 
control u(t) be given on the segment [T), T] and let x(t) satisfy system (1), Then the following theorem is valid. 


Theorem 5, If control u(t) satisfies the maximum (minimum) condition with respect to p(t) then the func- 
tion M(t) = H[x(t), p(t), u(t), t] is continuous and piecewise-differentiable on the interval (T,, T) where, at the 
continuity points of control u(t) the derivative exists and equals 


aM __ 0H |x (t), p(t), u(t), t) 


dt ot 





It follows from Theorem 5 that , if the right members of differential Equations (1) do not explicitly depend 


Of; - u) 


n 
on time,and 4 a dP pi —— = 0, then M = 0 and (with continuity taken into account) the function 





* In the problems ileus time T is not fixed, the theorem analogous to Theorem 4 apparently does not hold, 
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H[x(t), p(t), u(t)] is constant, If,in particular, H[x(t), p(t), u(t)] reduces to zero at time t = T as, for example, 
in formula (12), then the magnitude of H is identically zero for any t, 


The constancy of function H, analogous to the Hamiltonian function of analytical mechanics, corresponds 
to the law of the conservation of energy in conservative mechanical systems (the Hamiltonian function, as is well 
known, expresses the system's total energy). 


h. Singular Controls 





The recommended procedure for eliminating the control u(t) from Eqs, (3) by means of the maximum 
(minimum) condition may turn out to be impossible if the function H, for definite values of x and p, does not in 
fact depend on u, This may easily be followed in the simplest example of determining the control which, during 
fixed time T, takes the system described by the equations 


zt=u, |ul<M, 


T 
from point x(0) = x to point x(T) = x' and minimizes the integral ‘la \ x*dt, In correspondence with formulas 
0 


(11),we can write 


H = pu— 5-2°, p=. 











u(t) It follows from the maximum condition that, for p # 0, the con- 
PSK Se em = ——— trol equals u = M sign p, However, p = 0, the function H does not 
7 —F u depend on u, In addition, as may easily be seen, the optimal control, 
_ 2 coordinates and pulses, for the condition T > (x° 4 x) / M, must be 
x(t) changed as shown by Fig. 5. Therefore, on the time interval [t,, tg], 
r the function p(t) = 0, and the control u(t) is not defined by the maxi- 
X. ese aee mum condition, Such a control can be called singular, A singular 
¥ t control, if such exists, can be easily detected, Indeed, by first determin- 
pt} t, 7 ing for which values of x and p (in the given example, for p = 0) H 
Si "4 t does not depend on u, we can seek such a control for which this independ- 
Tat “ T ence would subsist in time, In the given example (since p = 0 =const), 
“tlie \ p= 0 and, consequently, x = 0, from which we obtain the singular con- 





trolu = 0, The function p(t) and, consequently, the control u(t), on the 
Fig. 5 entire segment [0, T], can be defined by the condition of continuity of 
the singular solution with the nonsingular solutions (in ghe given example, 
the solution for u = 0 with the solutions for u = 4 M), Theorem 5 can 
also be used, In the given case, for the singular control H = 0 and, for the nonsingular control, H = M| p | —x*/ 2. 
Since H = const = 0, we have 


- 


|P(O)| = sy (2, | p(T) = age (28, 


and the signs of the quantities p(0) and p(T) are quickly determined, 


i. The Synthesis Problem 





We now consider the methods of using the practice the maximum principle of L, S, Pontryagin for solving 
the problem of designing optimal systems, 


The problem of constructing optimal systems may be posed in two aspects. In the first case it may be suf- 
ficient to calculate beforehand an optimal control in the form of a function of time for definite, fixed and previ- 
ously given initial conditions, and thereafter to design a program-controlled system, In the second case, the sys- 
tem's initial conditions may be unknown beforehand, but, nevertheless, it is required that the system realize an 
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optimal process with arbitrary initial conditions, In other words, in this case,the optimal control must be expressed 
in the form of a function of the current value of the system coordinate and, possibly, of time, u = (x, t): 


Ux = Pe \Z,-+ +» In, ) (k= f,..-5 7), (17) 


which allows one, by the use of a regulator with nonlinear connections, to design a system which is optimal for 
any initial conditions, This problem bears the name of the optimal-system-synthesis problem. 


The maximum principle essentially solves the first problem immediately, since the finding of the optimal 
control in the form of a function of time amounts to the integration of a definite system of differential equations 
with the corresponding boundary conditions, It should be mentioned, however, that the technical difficulties of 
integrating the differential equations which flow from the maximum principle require, as a reule that computer 
technology be called into play, Integration of the equations in closed form turns out to be possible, naturally, in 
rare Cases only, 


The synthesis problem ordinarily arises in those problems in which the time of motion is not fixed before- 
hand, The following procedure for determing the functions (17) derives from the maximum principle, From the 
maximum (minimum) condition one determines the control in the form of the function u = #(x, p, t), With 
fixed initial conditions, x(T,) = x’, there are given various initial values of p(T,) = p’, so chosen that the solution 
of system (3) of equations, with initial cauiiiions x(T,) = x° and p(T,) = p’, is satisfied together with the required 
boundary conditions, If the value of p° is not thereby defined uniquely, one then chooses the set of values, 

Pi Ci = 1,. - « «» M), for which the functional's value turns out to be extremal, For each fixed set of values of 

xi (i ~ i. - « «, M) at the fixed moment of time Tp, ane is thus obtained a single set of values of pli =1, «i 
which provide the optimal process. The function p° =x (x°, T,) is thus defined, Therefore, at time Tp, the fol- 
lowing equation holds 


u(T,, 2°) = $[z°, x(x", 79), T,] 


But it is easily seen that this relationship must also hold at any arbitrary moment of time, which can always 
be taken as the initial moment, by virtue of which the control u is defined in a form of the type of (17): 


= [z, x (z, t), t]. 
In those cases when time does not enter explicitly in the right members of the differential equations, the 


function ¢, in (17) also turn out to be independent of time. 


We now turn our attention to the circumstance that among the variables x; there may be, not only the co- 
ordinates of the controlled object, but also artificially introduced quantities (for example, in the case of an in- 
t 


tegral extremum, the quantity xp 4 ;(t) = \ @ (x, u) dt). However, in those cases when the coordinates xs does 
T. 

not enter explicitly into the right members of the differential equations or into the relationships defining the 
domain G, the quantity p,(t) is constant in time and does not depend on the boundary values of the coordinates 
Xy, ++ +»Xp- In addition, x, does not enter into the function H and, consequently, does not enter into the expres- 
sion for the control u in terms of x, p, andt, Therefore, Herat (17) will also not contain xs. In particular, 
in the case of the optimization of the integral xp , ,(T) = -{ ® (zx, u) dt, only the object's coordinates, and not 
the quantity xp 4 1, will enter into (17). 

The synthesis problem becomes much more complicated if the artificially introduced coordinate, although 
not entering explicitly in the equations’ right members, does figure in the relationships which define region G, 

T 

Examples of this are the cases of "integral limitations of the type \ @ (z,u)dt <A, imposed on the control, 


T, 
Here, the introduced parameter In+1(t) = { D(z, u) dt enters into the definition of region G: Xp 4 a(t) <A. 


Therefore, functions (17) generally turn out to depend on x, , 3. In this case, in order to realize an optimal 
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system, one must either have a computing device which determines x 4 ;(t) at each moment of time, or one must 
previously compute x, , 4 in the form of a function of the current, xi = 1,...,m), and initial, xte coordinate 
values, 


The synthesis problem, to an even greater degree than the problem of determining the control in the form 
of a function of time, requires the use of computing technology. Even in those cases when it is possible to deter- 
mine functions (17) in closed form, the form of these functions is ordinarily very complicated, Therefore, a fore- 
most role is played by the problem of synthesizing systems which are sufficiently-simple and, in addition, are close 
to optimal, 


Another synthesis procedure, based on the solution of partial differential equations, derives from the dynamic 
programming method of Richard Bellman. 


We will speak of this in the following part of this work [Cf., Automation and Remote Control (USSR) 20, 12 
(1959) ]. 
APPENDIX 


1°. The Proof of Theorem 3 





We define aclassf of modified controls of a special form, On the interval (Ty, T), we arbitrarily choose 
j = 1 moments of time t = ty and we construct the system of j nonintersecting segments I, = [t, ~ T, t,] of length 
t, where the number r > 0 fs sufficiently small(k=1,...,j). We introduce j r-dimensional piecewise-con- 
tinuous vector functions ok (z) = ez), . + +» %,(Z)], defined on the segment (0, 1], on o(z) GU (k=1,.. J). 


In addition,to control u(t) we consider the modified control v(t; T, ty, . «» tfi Peics op ~) such that 
v = u(t) everywhere outside of the segments |, and v = oX ((t},.—t)/r jon kk. Control v(t) is obviously admissible, 
The points t, and the functions ¢(z) are called, respectively, defining points and defining functions of the con- 
trol v(t). The class © of controls includes all the controls v(t; T, ty, . « «» ti Finds #) with arbitrary j, t, 
oX(z) é U(k=1,..., 4) and r sufficiently small, 


Let x(t, [v]) be the trajectory corresponding to control€ I and let 6 x(t, [v) = x(t,[v)—x(t, [u), where 
x(t, (uJ) is the optimal trajectory. For x(t, [v]) we have the system of equations 


82, = f(z + dz, v)—f,(z,u) (i= 4,..., n), (18) 
which, by using Taylor's formula, we can rewrite in the form 


n 


a 
bz, = yo bz, + f; (2, *) —f,(z, u) + ay (t, [o]) (i= 4,.- +5) (19) 
1 § 


where 








a, (t, [o) = > Of, (z, v) _ Of, (2, 2) bs, +4 > Of, (x 4+ Sta, v) sae 


; Oz, oz, aed 0x ,02, .@ (20) 


and 0 < 9 < 1, By integrating (19) between the limits of Ty, and t(ty, < t< ty 417 T), we get 


tn m 
bz, (t, [v}) “ sae ba, dt + S\JE (lol) + By (t, (Dt Elms bmg — 7 (21) 
1 ° k= 1 


where 


tk 


7 ey =( [1 (=, e*(2—*))— Ae, « ©) ar, 


tht 

















( 5 88/, (2 + 92, »v) 


t2.t2.di + 
0x ,0%, . Ss 





T, ©9=1 


(- —t 
Of, | 2, o( : —*)) Of, (z, u) |, . 
zat 


+3 af — Se 











By taking into account the boundedness of df ;/0xz and 0°fj/ @x,xq on [To, T] and also the estimate for 
5 Xs 


| 32, (t,[e) |< Mt, tE&[T>, T] 

(where the constant M does not depend on 1), we easily obtain the estimate for 8{: 
1B, Ct, [*)|<SN*, tE[T,, T), 

where N also does not depend on t and r, by virtue of which, uniformly in t, 


lim Bi (¢, (el) (t, (ob) 


0 t site (22) 


We now consider the integrals i ([v}) (21). If we make the change of variables z = (t,— t)/r, we shall 
have 


1 
Fy ({v}) = * Ni (x(t, — 72), e* (2)) — f, (w(t, — t2), u(t, — t2))] ds. 
0 


As r-~ 0, by virtue of the continuity of x(t) and the continuity of u(t) to the left, i.e., u(t—0) = u(t), we get 


k 1 
“ = =(ue (ty)» @* (2) — fy (@ (tq), 4 (t))] as. (28) 


0 


R; (t,, (e*}) = lim 
tT 0 


It is easily shown that the limit y;(t[v) = lim, 5 x4(t,[v])/r (i =1,..., 1), which we call the variation 


of trajectory x(t) at the point t, exists for all t € [T,, T] different from t\%. Then, by dividing (21) by r and by 
taking (22) and (23) into account, we get 


tn m 
af, (x, 
y; (t, (el) a > a ae + DS) Riltys (2D, #€ ne bmg): (24) 
1 ° k=1 


It follows from (24) that, on the intervals (t,,,, tr 4 1), the functions y;(t[v]) satisfy the system of differential 
equations in the variations 


5 oh (@, w 1) 


1 on, . 


¥;= (i= 1,..., a), (25) 


and, at the point t = t,., experience the jumps 


Y; (ty + 0) — y(t, —0) = Ry (t,, [e*}) (i=i,... 








1 








where, obviously, y;(ty—0) = 0 (since 5 x,(t{v]) = 0 for t€[T,, y— T]). We denote by yet > t,) the solutions 
of system (25) with the initial conditions y;(t,) = Ri(tk,[ ¢]}. Due to the linearity of the equations, the solution 
of (25) with conditions (26) is obtained in the form 


Yj; (t, [v}) — >} y* (t), t € (tins bn-3)° (27) 


k=1 


We now determine the yi. If ¢Ury=[c4t),.. > 649) (q=1,..., Mm) is a fundamental system of 
solutions of (25), and the matrix || 9) l| is the inverse of |] ¢ ;*(t) ||, then 


vit) = >) CF) 98 (ty) v5 (ty)- 
8, q=1 


By taking the initial conditions into account, we get 


vi (t) = SAjy (t, ty) Ry (ty Ce"), 


s=1 


with the notation 


n 
Aj, (t, ty) =D) CF (¢) 98 (t,). 
q=1 


By substituting the expression for yk (t) in (27) and introducing the additional notation By.(t,) = A(T, t), 
we obtain, for t = T 


j n 
yy (Lr) = (7, (21) = >) DY) Buslty) Ry ty lo“) (i= 4,.... 9. 
k=] 8= ] (28) 
The matrix B,,(t,) is defined only by Eqs. (27) in the variations, and does not depend on the choice of the 
functions o*(z). The vector y([v]) = (y3, . . +» Yn) is linearly expressed in terms of the vectors R" = (Rh, ee oo RK) 
(k=1,..., 4). 


The points of phase space of the form x’ + y([v}, where ver’, form a certain set Il, We now prove that Il 
is convex. Let y([v"]) and y([v")}) be arbitrary variations of the trajectories obtained with the modified controls 
v’ and v" with defining points and functions t,, ¢“(z)(a=1,...,3") tB, g™(z) (8 =1,.. ., J”) respectively, 
By definition 


rea2i+yi([er’"Jel,  =z'+y((e")eEn. 
We show that any point of the form 


eS yn’ + (1 —p)e* = 2! + py fo’) +(1— vy (lo"), OSH <1, (29) 


which lies on the line segment joining s and w" also lies in I. For this, obviously, it is necessary to show that 
there exists a modified control v¢f such that 


y ({2}) = vy ([e’) + (4 — ») y’ (Le")). (30) 
We define control v in the following manner. As the control's defining points we take the ynion of the 


points t', andt%. The number of defining points of control v thus equals j' + j"—h, where h is the number of 
pairs of coincident points t' and t"g. The defining points of control v are divided into three classes: 1) the 








points which coincide with the defining points of control v’ and do not coincide with the defining points of con- 
trol v"; 2) the points which coincide with defining points of control vy" and do not coincide with defining points 
of control v’; 3) the points which coincide simultaneously with defining points of controls v' and v", For the 
points of the first class, t) = t'.,;(k and a, run through j'—h vaiues), we give the defining functions of control v 
in the form 

ets) O<:<p 


Hi—| ; 
u(t, ) p< 2zgi. 


For the points of the second class t2 = t"g, (I and 6, run through j" —h values) 


o® ( - ) O<s<i-—p 
Ha)=| i=? 
u (t,.) i—p<sq1. 





For points of the third class, try = t'yz = t"g2(m, %, 8 run through j values) 
@* (=) 0<s<ap 
p» 


G=4) <r 


9” (s) = 





By computing Ri(t,,[*J, Ry( ty » [ ¢'Dand Ri(tm, [¥™) by formula (23), we find that 


Ry (ty, (9"}) = wR, (t,, [9D 


Ry (ty, (e")) =(1— w) Ry (tg, fe), 
Ry ltenr 19") =P R(t, (91) +(1—v) Ry (te, 9°”). 


If we now determine y;([v]) by (28), we convince ourselves of the validity of relationship (30), We also 
note that x'€Il since, by giving the modified control v(t; T, ty... os tki Gis « o g*) such that; at each defining 
point tk, the defining function gX(z) = u(t,), we obtain y,([v]) = 0. 


Thus, II is a convex set which includes point x’, For definiteness, we consider further the case of a mini- 
mum of the functional S, The sets G™ and I have the common point x’, We now show that, in the nondegenerate 
problem, none of the points g~€ G™ (in particular, the point x4) can be an interior point of II, Indeed, let, on the 
contrary, point g~ be interior to Il, If the problem is nondegenerate, the set G~ will have interior points, Other- 
wise, hyperplane C would be a bracket for set G, and for any point x = (x,,.. .,X,) © G, the inequality 


n 

ye (z,— zl) > 0, would hold, this being the condition for a degenerate problem, Therefore, without any limit- 
1 

ation on generality, we can assume that g™ is an interior point of G” since, if g’ were a boundary point of G~, 
one could always find another point interior for both G™ and I (by virtue of the fact that, in any arbitrarily small 
neighborhood of any point in G~ there are contained points which are interior for G~), We then find some neigh- 
borhood E of point g~, all the points of which are interior in both 1 and G . We choose the modified control 
v¢r sucit that x’+ y([v]) = g (this is always possible, since g~€I1), The end of the trajectory x(t,[v]) is the point 
xi4 5x(T,[v), where 5x(T,[v}) =ry([v) + Te(T). Since lim, €(r)= 0, we can always choose a r sufficiently 


small so that the point g = x) + y(iv) + €(r)=g + €(7) lies in E, But points of the form (1- p)xt+ ug, 

0< pw <1,g€E, are interior points of G-. By choosing yw = 1, we are led to the conclusion that point x! + 
+ 5x(T,[v]}) is also an interior point of G, i,e., in other words, control v(t; T, ty, . « «» t; . ¢) takes 

the system to the interior of set G . This contradicts the optimality of control u(t) and point g~ €-G~ cannot be 
an interior point of II. We show analogously that point m€ Il cannot be interior of G~. 


It thus follows tnat the convex sets [1 and G can have in common only boundary points, one of which (or 
the only one of which) is the point x’, where one of the sets G~ has interior points in phase space X, Therefore, 
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hyperplane A, defined by the expression >} (x, — z!) =0,, may be passed through point x’, thus separating sets 
os 

G and II. The signs of coefficients as can be so chosen that the set II is located on that side of hyperplane A 
n n 

where’ p a, (z,—z!)>0. (set G” is then completely located on the side of the hyperplane where d\% (z,—2)< 
1 tf 1 

=0). Therefore, for any control v(t; T, t,.. «» ty: Fo Sas eér, the inequality De. ({v]) >0, will hold, t,e, 

1 


n 

>) a,3z, (7, [r}) 
lim + > 0. 
™0 t 


n 
% > a,bz, 


We now consider the functional S= >) a,x, (T). Then, the requirement that lim, _____$._ > 0 means that 
4 





lim “° >0, (31) 
T+0 


As a modified control, we select the control y€Ir with only one defining point ty, such that 5u(t) = v(t)— 
— u(t) may differ from zero only on the segment [ty— T, ty]. For the control v(t) = u(t)+6 u(t), inequality (31) 


n 


holds (as for all controls in’), By now introducing the vector p(t) and the function H (z, p, u, t)= Pili 
1 





(Xq, 5 + +» Xi Use « o» Ups Ct), which satisfy the relationships 
oH : 0H : 
, = , =— =1,..., ”), 
i“ ap,’ dz; ") 


n 
and by setting py(T) = —aj (i=1,..., 1), we can use, for the functional § = >} as*, (T), the formula for the 
i 
increment of functional values (Appendix to Part I), By repeating virtually word for word the argument in the 
proof of Theorem 1, we can convince ourselves that condition (31) can hold only in the case, when at any moment 
of time t, the function K(t,u) = H[x(t) p(t), u, t) assumes its maximum value with respect to u, 


n 
The case of the maximum of functional S = >} c,x, (7) is treated analogously, Thus, Theorem 3 is proved, 
1 


2°. The Proof of Theorem 4 





We carry out the proof for the case when one seeks a minimum of the functional S. The case, where the 
maximum of the functional is sought, is treated analogously, 


Let Wr be the set of points in phase space X which are attainable from point x® at time t = T for all the 
possible controls u(t)€ U. We consider the intersection Q of sets G and Wr: Q= GW. The control u(t) and 
the trajectory x(t), corresponding to it, are min-optimal if, for each point x(x, . . ., Xp) EQ, the inequality 


n n 
>) (7) <>}, 2;-holds, By denoting x;(T) = x}, we get 
1 1 


n 
dy % (z, — z}) > 0. (32) 


1 


To prove Theorem t, it suffices to show that, if the conditions of the Theorem hold, then inequality (32) is 
satisfied for each point x €Q. 
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We use the notation 


a; = re, + ph, (i = peers n). 


If the control u(t) satisfies the maximum condition with respect to the vector p(t) which, at time t = T, 


assumes the value 


p(T)=—a, a= (a),..., ay), 


then, in accordance with Theorem 2 of the present work, control u(t) is min-optimal with respect to S = 
h 


= D7 (T)in the problem with free trajectory right ends, i.e., for each point x¢ Wy, the following inequality 


is valid 


n 


yal z!)>0. 
1 


On the other hand, for each point x ¢G the following inequality holds, by definition, 


n 
Sel. (2) D0. 
1 


Therefore, for each point x“¢Q = G!]W 7, inequalities (34) and (35) hold, 
But it follows from (33) that, for A > 0, 


We now choose an arbitrary point x ¢ Q and set up the expression 


n 4 n 

. . 
% (z,—2z)= xD" (z,;— 2}) — ‘ 
1 1 


\ b; (x, — 23). 


»Ma 


By virtue of inequalities (34) and (35), inequality (32) follows from this. g.e.d. 


3°. The Proof of Theorem 5 





For definiteness, we consider the case when the control satisfies the maximum condition, The case of 


the minimum condition is handled analogously, 
We introduce the difference 


AM =H [y(t + At), u(t+ At), t+ At]—A y(t), u(t), t). 


Now, to shorten the formulas, we introduce the vector y = (yy,.. ., Yen) (¥1 = Xu; Yn 41 = Pps 1 = 1; 


By virtue of the maximum condition 
H ly (¢ + At), u(t + At), t+ At] >A [y(t + At), ult), ¢+ Ad}, 
A ly(), u(t), t] > H [y(t), u(t + Ad), ¢). 


From whence follows, for AM, 


H |y (t+ At), u(t), ¢+ At] —A [y (t), u(t), tS AMS 
<H [y(t + At), u(t + At), t+ At] —A [y(t), u(t + As), t]. 


o oan 


(33) 


(34) 


(35) 


(36) 








Because of the continuity of y,(t) (s = 1, . . ., 2n) in time,and of the continuity of H(y, u, t) over the set of 
arguments, for At + 0, the extreme terms in inequality (36) have the limit zero, Consequently, lim 4M = 0, which 
proves the continuity of the function M(t), 


Now, let t be a continuity point of control u(t), Then, Afra att + At) = u(t) and, moreover, the derivatives 


ys(t) (8 = 1, ..., 2m) exist and are continuous, By dividing (36) by At and letting At approach zero, we convince 
ourselves that the limits for the extreme terms of the inequality exist and equal 


an 
OH {yit), u(t).t] 3 0H [y (t), u (t), t] 
> oy, ¥.(t) + ot y 








1 


By virtue of the relationships, y, = 0H/dy, 4 1, yn+1= ~OH/dyj(1=1,. . .» m), the sum in the pre- 
vious expression reduces to zero and the limit of both extreme terms equals 0 H/ét, thanks to which, dM/dt = 
= 8H/dt, Consequently,the theorem is proven, 
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ON THE SYMMETRIC PERIODIC MOTIONS OF MULTISTAGE 
RELAY SYSTEMS 


Yu. I, Neimark and L, P, Shil'nikov 
(Gor’kil) 


Many works [1-15] have been devoted to the exposition of various methods of seeking, and investigating 
the stability of the simplest symmetric-periodic motions of linear systems, Systems with two relays [13, 16] and, 
in particular, a two-stage relay system [17], were created with these methods, However, the complexity in the 
use of these methods increases with an increasing number of stages in the relay system under consideration, 


In the present work, we solve the problem of finding, and investigating the stability of, the simplest sym- 
metric motions of multistage relay systems, The solution is attained by the use of the so-called phase character- 
istic [15] which is of independent interest, 


1. The Phase Characteristic 





Let a relay with the symmetric characteristic,represented in Fig. 1 be connected to the output of a linear 
link, If an infinite series of symmetric rectangular pulses with half-period T is applied to the input of a linear 
link of this system then, at the relay output, there is produced some series of rectangular pulses. For some inter- 
vals of values of the half-period T, the series of output pulses will also by symmetric and periodic, with period 
2T, but with the moments of switching shifted with respect to the moments of switching of the input pulses by 
the quantity 9(T), The curve @(T) will be called the phase characteristic of the open-loop relay system under 
consideration by us. 


As is well known, the output of a linear link may be expressed in terms of its transfer conductance function 
@(t) in the form [12-13] 




















co 
y(t,T) = 2X (— 1) Ip (t—(&—1) 7) — 9 ¢— AT) + 90. al) 
k=0 
The moment 6 (T) of relay switching from —1 to 1 which is closest 
I to time t = 0 is found as the least positive root of the equation 
y(t, T) =8, (2) 
7 for which the conditions of switching hold, In the case when y; is continuous, 
e y the switching condition can be written in the form 
y, (t,T) > 0. 
Fig, 1 


In seeking the phase characteristic 6 (T) in practice, one should use 
formula (1) if the series in (1) converges rapidly and the necessary accuracy 
can be attained after not too many of its terms have beensummed, If the transfer constant K(p) of the linear 
link can be represented to a good enough approximation in the form 
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A te 





K (p) = J, (3) 
7 


then, as is well known [12], the linear link's output can be written in finite form: 
c; h.T c 
. j rts 
y(t, T) =— a ej th —-~ + Ye (eit —1). (4) 


Finally, it is well known [14] that the sum in (1) can be transformed to the form 





pena $5 eleva) a 
m=1 
x sin| (2m —1)F t+ arg K (ie *), 
P r eee 


which it is natural to use if series (5) converges 
sufficiently rapidly, i.e., if the quantities 


| 

i 

' 

1 7 3 

K(1 =) K(1 Terre tend to zero with sufficient 
T Tho, Inj, rapidity. 











r) @ We give the easily found phase character- 
a istic of standard linear links of automatic control 
‘ f systems (Fig, 2) with relays for which 6 = 0, 
7 | 7 We note that when a lag r appears in the 
‘ r(9 linear link, the phase characteristic © (T) is re- 


placed by 6 (T) + r—[(r/2T2T, where [r /2T} 
denotes the greatest integer less than r/ 2T. 


We now consider the open chain consisting 

of m series-connected linear links and relays, 
Let 6,(T) be the phase characteristc of the first 
linear element and the relay following it, 6,(T) the 
phase characteristic of the second linear element 
and the relay following it, etc, Obviously, to with- 
in a multiple of the period, the phase character- 
istic 6 (T) equals 6,(T) + 99(T) +... + 9_(T) 
or, if we understand by @(T) the moment of relay 
Fig. 2 switching which is shifted in time by less 

' than 2T: 











o. 
Ss 
4 
\ 
Be 





~ 
+ 
~ 
3 
% 
x 
~ 


0:(T) +... +,,(T) 
a 27. (6) 





(7) =0,(7) +... 4+%(7)—| 
where [ ] denotes the integral part of its argument, 


2. Seeking the Simplest Symmetric-Periodic Motions of a Multistagé Relay System 





We consider a multistage relay system consisting of m series-connected linear links and m symmetric 
relays, The phase characteristic 6(T) corresponding to the open-loop system can be found in terms of phase 
characteristics of the constituent single-relay systems in the form of (6), Obviously, the closed-leop relay system 
will admit the simplest symmetric motion with half-period T if only 


6(7) =0, (1) 
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and, consequently, the search for the simplest symmetric motions reduces to finding the points of intersection 
of the curve 9 °(T)= 6,(T) +... + ©,,(T) with the lines 6 *(T) = 2kT, where k= 0,1,2... 





























Transmission , 
xa | factor of the Output equation of the linear link [phase characteristic fee 
linear link in the interval O<t<T 6 (t) oh No, 
ee | 
l I T 
1 — Tig iof- —— 2, a 
Pp 2 2 
, t t 
Ppp : , 
2 Te +1 —€¢ th Ty eid |- T, In me a 2, b 
ite ‘ 
; T 
e t TT. 
| ah Th ite & 
T 2 
4 oa Big de T 
Pp , 2 
i 1 AT Ls . § T Anal l oh 
5 pA Tat SS + ne. ae a es es nalytical phrase 
P(p +A) h (tn ag t)e +X 2+! lig evidently absent » g 
Pp 2e"T—H (8T sin @,f — sin m, (7'—?)) | 4 Sin 7 
6 er =< aT TS “Garo See 
PP+2bp + ov w, (1 + 2e”* cos w,7' + e*") 4 + cos w,7 e 
For w? = 02 —b?> 0 
2e%T— (6°T sh Ot — sin Q(T — 2) | 1 th <hr? . 
(1 + 2eF char + eT) iy Ore Ty cha®r | “ 
For 0? =t?—w?>0 
—pt au» Po a] ‘ 
For 2 = 0, the equation of the phase characteristic will be T/(1 +e 6 T), 


8. Investigation of the Stability of the Simplest Symmetric Motion of a Multi- 
stage Relay System 
The investigation of the stability of the simplest symmetric motion of a multistage relay system may be 


easily carried out by the method described in [12, 13], We denote by 4 (9 the stimulated moments of switching 
of relay i, and we write the linear link's output in accordance with (1) as 








—0o 
ya(t) = Dy (— 1) Lei (¢ — 1.) — i (tt) + (— 1)" (¢ — 1), (8) 
k=n j 
where ¢ is the transfer conductance function of the i'th linear link and n is the ordinal number of last relay 
switching up to time t. Let nT + o,..., 0T + Gp, be the switching moments of, respectively, relays 1, 2,... 
. « «» m which lie in the interval [t— T,t], The following relay switching in the interval (t, t + T] will occur 
in correspondence with the order of increasing «;, Therefore, the order of relay switching may be presented in 


the form of the substitution (' Bees = By virtue of this, linearized Eq. (8) is written in the following form: 


r Tg... «lm 


ait? Dd (— 1) Ig (a + 1) T — (kK A) T + igs — ou) — 


k=n+€; 


— 9, ((n + 1) T —AT + a4, —a4)] + 


1424 








+ (— 1)" — 8p (1 — ©:) T + oi4s — oH) + 
+ 2(—1)*t8 —*p, ((2 — e:) T + a4, — m4) dtl es + (9) 
+2(— 4)ntes 89 ((83— 2) T + ais — a) dt 1p t+--=0 
(i == 1,2,...., m), 


where €, equal zero or one depending on whether rj > rj + 1 Or ti < Tj 4 4. Equation (9) allows one to determine 
(i+ 4) 


the variations dt’, , ;° of the moments of switching of the ( i + 1)'st relay as functions of the stimuli ar”) +e’ 


atDise jeeee of the i'th relay, From the conditions for relay switching in the steady state, the following 
relationships are easily established: 


0; for 6<T 
1 = (1—e8) 7 + a4, —a, = (0,—T for 06,>T7 (10) 
T for 0; —0. 


By substituting ar’) = off), nth in Eq. (9) and then eliminating the unknown of) of?), oe cA™) | 


we get to the characteristic equation in the form 


x (2) = x1 (2),- + +s Xm(z) —1 = 9, (11) 


where 


+00 : 
>) (—1)* ey (kT +4) 2* 
k=0 


xi (2) = —S 
> (—1)*-* 9; (kT + 1) 


k=0 





The stability condition for the periodic motion under consideration consists of this, that all but one of the 
roots of Eq, (11) lie outside the unit circle |z| <1 [12, 13}. 


If the transfer constants of the linear links of the multistage relay system under consideration are written 
in the form of (3), then the characteristic Eq. (11) can be written in the explicit form: 


Avi —1 


~ C5 ze ti V4 Ci; ¢ 4 
x (2) = [] 3 — iF (x iF ) —i =0. (12) 


i=x] 1+ ze‘ i+e 





In case the linear links contain lags 1;, the functions x ;(z) entering into characteristic Eq, (11) will have 
the following form: 


8, +1 Aj jBi - 594 —l 


ui 
_ S$ %j * ¢ Sy My 
x (8) = 2 + se (2 44 @tit ) ? (18) 





where s, is the least positive integer for which sjT +yj— T; > 0. 


If the series entering into (11) converge sufficiently rapidly for | z| = 1, then numerical investigation of 
stability may be carried out, for example, by means of the reductions described in [13]. 


In case the frequency characteristics Kj(i w) of the linear links are given, the characteristic Eq, (11) by the 
use of the relationship [14] 
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- ape m ti 
>> 9; (2kT + t) ete = 5 — yy» e K,(m+.i—2 i), (14) 


k=—©co n-=—o 


may be given in the following form: 


x (e*) = x (e'*) .. . xm(*?) —1 = 0, (15) 
where 
xj (ce?) = i-¥) y gil  K, or —f) x 
m=—oo 
x eo Paamial i K{ (2m - ‘ys )) , (16) 


If series (16) converge sufficiently rapidly, then investigation of stability is conveniently carried out by the 
numerical construction of the boundaries of the stable regions which are obtained as ¢ varies from 0 to 27. 


4. Examples of the Construction of the Phase Characteristics, the Finding, and the 
Investigation of the Stability of the Periodic Motions 








1, Let the graph of the transfer conductance of the linear link be given, In this case, the phase character- 
istic can be constructed by the following method, For each fixed value of the half-period T we find the output 
of the linear link, for which we geometrically sum the reactions of the linear link to the rectangular pulses in 
accordance with formula (1) and we determine the moment of switching of the relay placed after the linear link*, 
By constructing the outputs of the linear links and determining the shifts of the moments of switching of the relays 
for various values of T, we can find the graph of the phase characteristic, Thus, for the transfer conductance 
function, represented on Fig. 3,a, the phase characteristic found by this method has the form shown on Fig, 3,b. 














a 
a, b 

a ) b Hy 0 

3 as 
Miwy A 45 — - 
fee Te ; ; 0 ews ws wt 

t T 
Fig. 3 Fig. 4 


2. If the frequency characteristic K(iw) of the linear link is given then, to find the phase characteristic, 
it is convenient to use formula (5), In this case, the construction of the linear link's output reduces to the geom- 
etric summing of the reactions of the linear link to the harmonic frequencies w/t, 3n/T,54/T,,.., the amp- 
litudes and phases of which are directly determined from the frequency characteristics as the modulus and argu- 
ment of the quantities K(it/T) , K(3 mi/ T). By determining the shift 6 for different values of T, we can con- 
struct the graph of the phase characteristic, Thus, for the frequency characteristic shown in Fig. 4,a, the phase 
characteristic has the form shown in Fig. 4,b. 


3. To illustrate the theory presented in this paper, we consider the m-stage relay system whose linear links 
have one of the following transfer constants, e" P7/(T,,p + 1), e"P™/(T, p— 1), e-P7/ p._In accordance with 
Section 2 of this paper, the half-periods of the possible periodic motions are found from the condition that the 
phase characteristic of the open-loop relay system intersect the lines 9= 2kT, k= 1,2,... Thus, for the four- 
stage relay system whose linear links are identical and have the transfer constant e~P7/(T;, p + 1), the half- 
periods of the symmetric periodic motions are found from the condition that the phase characteristic of the linear 
link intersect the lines 0* = kT/2,k =1, 2, 3, 4 (Fig. 5). 

















If the number of stages is greater than two, all the symmetric motions will be unstable, Indeed, the char- 
acteristic equation for such a system, after some simple transformations, can be written in the following form: 





EIN — 4) 5 ay(ENt™—1 1) +... + ay ym (E—1) = 0, te 


where a; > 0. Dividing the left member of Eq, (17) by ¢ —1, we find that 


ENtm—1 +. + bys mae + 1+ B=, (18) 


where B > 0. Since the product of the roots of Eq, (18) is equal to 1 + B, there are roots which do not lie within 


the unit circle, and, consequently, the symmetric periodic motion is unstable, By virtue of this, in the multi- 
stage relay systems considered here for which m= 2, only more complex periodic motions are possible, 


\ 


(W) 





Fee 
i 


; 





Fig. 6 


We now consider an m-stage system whose linear links are also 


identical and are described by the transfer constants al we - 
p + 2 + 


—5*=wi > 0, In this case, the half-periods of the possible periodic 
motions may be determined from the conditions that the phase character- 
istic of the linear link, the form of which is shown in Fig. 2,e, intersect 
the lines 6° = 2kT/m,k=0,1,... Thus, fork =0, it is clear that 
the periodic mode with half-period #/ wy, is possible in the system, 


The question of the stability of the periodic modes in a system 
with identical linear links leads to the investigation of the roots of the 
equation 


(x (2))" — 1 =0. (19) 


In our case, each individual factor of this equation for the periodic mode with half-period # / w, is written 


in the form 


where z = 1/€ and W takes the values e 











~~ . 

ihe 1 

are = =W, (20) 
i_e 

“120t/M »2400,....m. 


By substituting € = e'!? in (20) and constructing the graph of the regions of stability for the parameter W, 
we find that for any number of links there will be stability, since the unit circle lies completely inside the region 











of stability (Fig. 6), The stability condition for the periodic solutions for k # 0 leads to the requirement that the 
roots of all the equations 


B 
g (— oO sin @, (7' +-0)+ cos , (T +0)) +o-* (~ =. sin 0,9 + cos w,0) 








£2 + 2e-*! cos wTE + eF lig 
x ; 1+ 2e—*T cos wT +e—*T ; = (21) 
— >, sino (T + 9) + cos , (7 + 6) + eT (_ — sin 010 + cos @,0 ) 
= 
=e ™ (7 =1, 2,...,m) 


of the second order in € , except one, equal to unity, in the equation with r = m, lie inside the unit circle, 
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APPROXIMATE DETERMINATION OF THE AUTO-OSCILLATION 


PARAMETERS IN RELAY SYSTEMS WITH LAGGED FEEDBACK 


N. A, Korolev 


(Moscow ) 


An approximate method is presented for determining the parameters of the 
auto-oscillations of one class of relay systems on the basis of a modification of 
the harmonic balance method, 


The harmonic balance method [1, 2] is the most widely used method of determining the auto-oscillation 
parameters in relay systems of automatic control, An exact method of determining the periodic modes was 
developed for relay systems [3, 4], It turns out that an entire group of relay systems can not be subjected to 
investigation by the harmonic balance method [5], namely, systems with lagged feedback (LF), We now give 
an elementary example wherein the harmonic balance method does not give the correct answer to the question 
of auto-oscillations, For this, we consider the simplest saw-tooth generator, the block schematic of which is 
given in Fig, 1, a, It consists of a linear portion, in the form of an aperiodic link, and a relay element with 
hysteresis kK». The solution by the harmonic balance method does not give an intersection of the frequency char- 
acteristic W(jw) with the curve of the nonlinear impedance 1/ (A) (Fig. 1,b), i.e., it shows the absence of auto- 















































oscillations, 

r J jf 2(t) 

Af 0 9 *t> 
Zz y “Zo & 
Xo Yi(A) W, 
ptt b (je) Te) 
a 
Fig. 1 


The solution by the frequency (exact) method [4] gives the exact value of the frequency and "amplitude" 
of the auto-oscillation (Fig. 1,c). In the given case 





Im J () = —2(=) = — ky th => = — x, (1) 
whence 
7 
oe 2 
” ig — %o 


The maximum deviation z(t) from the equilibrium position is. 
[z (t)]max = Xo» (3) 
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(3) 





i.e,, the auto-oscillation occurs inside the relay's hysteresis 



























































f ro] pg y loop (Fig. 1,d). 
y { By a relay system with LF we mean a system with an 
- LF » i’ arbitrary linear portion (LP) (Fig. 2, a) in which the relay 
: -%,| 0; {2 element (RE) is shunted by an internal feedback path (LF) 
: LP be in the form of an aperiodic or integrating link [5], One 
; b can sometimes reduce to such a system a relay system in 
. which the derivative of the time characteristic ho) >0 
Fig. 2 (for h(0) = 0) or, what amounts to the same thing, the dif- 


ference in the degrees of the denominator and the numerator 
of the linear portion's transfer function W(p) equals one, For the concrete examples in [5] there were given the 
qualitative and quantitative differences of the results of computing the auto-oscillations in a system with LF by 
the exact method of Ya, Z. Tsypkin [4] and the approximate method of L. S. Gol'dfarb [1], The relays were 
chosen without dead zones (Fig. 2,b). For convenience in comparing the results of this paper and those of [5], we 
consider the same numerical example of computing the auto-oscillation parameters in a relay system with LF: 


W,(p)= = C 
. P(Tip + 1)(T ap + 1) p (0.05p + 1)(0.03p + 1) ° 
1 





(4) 





ke 
W3(P)=7o41 ~oulp tt’ 


where W,(p) is the transfer function of the linear portion of the system without feedback and W,{p) is the feed- 
back's transfer function, 























JO sol’ 

we 

a 3 
b 
Fig, 3 


on . om 4k 
On Fig. 3,a there are weustracted the frequercy characteristic Wofjx') + \vy/te ) + WA) = ‘e Wo (iw) and 
the reciprocal of the relav’s nonlinear gain KA), 


Figure 3,b grves the ¢ harecteristic} 9(‘) of the relay system as the sum of the characteristics of the inherent 
system J,(w) and of the internal feedback path J,(u’) (the “left"characteristic), 


Figure 3,c gives, in the same scale, the characteristics Wy(jw) and Jo(w) = J,(w) + Jw). As we see, these 
characteristics are very remote from each other, 


By using the well-known formula [5] 








nO) 4 wy Vim —1\o} 
T(0) = by {0 + 4 Su em—tyo 4 + YS (5) 
m=] m=l1 
and by limiting ourselves to the first harmonic, we can write 
—(—) 
J (w) = — ky *-. 4 ate W (jo) = : W (jo) —k, hO) -- W' (jo). (6) 


It follows from this that a first approximation (in the sense of a first harmonic) to J(w) differs from the given 
characteristic Wjw) by the quantity Kk, h(0) / w, completely defined by the ordinary parameters of the lagged feed- 
back, since h(0) = kp/ Ty (5). 
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Figure 4 gives the characteristic WGw) of the system 
j described by Eq. (4). There is no difficulty in obtaining it from 
-2 -f the frequency characteristic, and the results of determining the 
(~ - get 0 auto-oscillation frequency by its use are significantly closer to 
, “sat sido aT Tl the truth, By remarking that the entire curve of W~)(jw) lies in 
30 the left half-plane, we can assert that there will be auto-oscil- 
0 lation in the system even when the relay's hysteresis kK» is such 
Jo 00 4-4 that a > Ky > 0 (a is shown on Fig. 3,c), In these limits of Ko, 
w 20 f the harmonic balance method shows an absence of auto-oscillation 
(since the semiline —1 MA) =— Kg lies completely in the left 
half-plane), However, in those cases when the Gol'dfarb method 
Fig. 4 does give some result in determining the auto-oscillation frequency, 
this result is not changed by replacing Ww) by the characteristic 
W Guy, since the auto-oscillation frequency is defined by the characteristics’ imaginary components, and these 
are identical (6), 











The quantitative difference in the solution by the exact and the approximate methods, as before, remain 
essential, and this is explained by the sharp difference of the form of the auto-oscillation from the sinusoidal, i.,e., 
by the low filtering properties of systems with LF. 


It was shown in [4] that the determination of the auto-oscillation frequency by the harmonic balance method 
in relay systems of the type in (4) without internal feedback (ky = 0) gives quite accurate results, This is alsoclear 
from Fig. 3. By determining the auto-oscillation frequencies in the points of intersection of the line — Ky with 
the curves W,(j w) and J,(w) respectively, we obtain close values of w»' and w"», This fact shows that the linear 
portion of system (4) itself is a good low-frequency filter. The inclusion of a lagged feedback changes the picture 
siginificantly, Higher harmonics in the periodic curve Z (t) appear due to the feedback path (this assertion is 
easily verified by the construction of Z ,(t) and Z,(t) by the well-known rules [2]), since it is described by a first- 
order differential equation, It therefore seems advantageous to take into account, in the approximate solution, 
the first harmonic of the signal Z ,(t) passing through the fundamental circuit and the exact form of the signal 
Z(t) passing through the lagged feedback path. The function 


F (jo) = W, (jo) + J2(w) = J, (o) (7) 


will satisfy the problem posed, Its construction is much simpler than that of Jy(w), For a lagged feedback in the 
form of an aperiodic link as in (4), the construction of the characteristic 


4 , 
Ja (@) = —keghp jap (1 — th 55) + 7th 5-5 —| (8) 
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poses no great difficulties, For a lagged feedback in the form of an integrating link, W,(p) = ky/p, the character- 
istic 





k ky + 
J,(@) = — = (1475) (9) 
is a straight line, When the transfer function of the system's linear portion, W;(p) is unknown, and only its fre- 
quency characteristic is available, W,(jw), the characteristic F(jw) is determined by the direct graphical addition 
of W(jw) and Jxw). 


The characteristic F(jw) for system (4) was constructed on Fig, 4. As we see, it differs slightly from the 
exact characteristic J9(w), and the error in determining the auto-oscillation frequency is insignificant for any Kk» 
or, what is equivalent, for any gain ky in the internal feedback path, 


We also note that the harmonic balance method can give erroneous results even in computation for those 
relay systems in which h(0) = h(0) = 0. As an example, we take the system with a transfer function W,(p) as 
given above in (4) and with an internal feedback transfer function 
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wo 40 so = Figure 5 shows the components, W,( jw) and 
W jw), and the total, We(jw) of the given system 
frequency characteristic, It is clear from the figure 
that, for certain values of w, the characteristic 

W (jw) enters into the right half-plane, i.e,, the 
harmonic balance method shows, for certain values 
of the relay element's hysteresis x 9(k» * 2), an 
absence of auto-oscillations in the system, In fact, 
“Jt however, there are auto-oscillations for any Kg> 0, 











" as may easily be verified by constructing the relay 
system's characteristic Jy(w), The qualitative dif- 
Fig. 5 ference between the results in the present case, as 


previously, flows from the essential difference of the 
auto-oscillation form from a harmonic (sinusoidal) one, In this case one can also use the suggested method of 
determining the auto-oscillation parameters, It is necessary to compute J,(w) on the basis of (10): 
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after which Eq, (7) is used for constructing F(jw). Figure 6,a shows Jw) and F(jw) (where we have taken k, = 1) 
and, on Fig. 6, b, F(jw) and Jg(w) are shown, As is clear from Fig. 6, b, the difference in the determination of 
the auto-oscillation frequency wp» for any K, is very small, 
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SUMMARY 


1, The harmonic balance method does not give reliable results for the group of relay systems with lagged 
feedback, 


2. The suggested method of taking the exact characteristic of the lagged feedback into account allows 
one to make direct use of the linear portion's frequency characteristic in investigating relay systems with lagged 
feedback. 
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A MULTICHANNEL AUTOMATIC OPTIMIZER FOR SOLVING 
VARIATIONAL PROBLEMS 


R, I, Stakhovskii 


(Moscow) 


The block schematic, as well as several functional schematics of the 
component blocks, are described for an electronic variant of a 12-channel 
optimizer, Results of testing it as an automatic approximater of a given func- 
tion of time are provided, A search method is given when “false” extrema, 
of the "ridge" and “saddle” types are present, 


It is well known that there is great interest currently in the development of devices which permit the find- 
ing (and retaining, if the external conditions thereafter change) of the most advantageous mode of operation of 
productive aggregates, There arise the problems of synthesizing the best, in some sense, control system and of 
best adjusting the aggregates, The automatic solution of such problems is provided by automatic optimization 
systems [1]. These consist of the object to be optimized and a computing device (an automatic optimizer) which 
carries out a search of the object's input parameters in such a way as to obtain an extremal of the possible values 
of some indication, With this, the presence of additional limitations may be taken into account [1], It wasshown 
[2] that the most efficient algorithm (program) of search is a union of the methods of steepest déscent and gradi- 
ent, this giving the possibility of joining the rapidity of the initial phase of the search with the accuracy of its 
final phase, Also demonstrated [1, 3] was the effectiveness of dividing the optimizer into two basic blocks, 
namely, an operational (analog block) and a controlling (digital) block, and a number of simple schemes for 
the controlling block, implementing the chosen search algorithm, were developed [3]. This paper is a further 
development of the work cited, and, in particular, is devoted to the description of, and test results for, an elec- 
tronic variant of a 12-channel automatic optimizer, 


Short Description of the Block Schematic 





The instrument to be described (Fig. 1) is designed to find an extremum of a function of n variables 
(n < 12) given m limitations (m <4), It is assumed with this that the function whose extremum is to be found, 
Q(x,, . . -» Xp) possesses one exact extremum, and on its form, as on that of the limiting functions Hj(X1, - - +» Xp 
there is imposed the condition that the first-order discontinuities (if there are such in the function) do not exe 
ceed a given magnitude [1]. The basic task of the optimizer variant to be described is to operate with high- 
speed models (frequency of solutions up to 30-50 cycles, and thereafter up to 100 solutions per second) while all 
its elements are to be purely electronic, Certainly, the variant to be described can operate with the ordinary 
models (for example, model EMU-6). 


The search algorithm used in the optimizer is capable of seeking extrema of functions of, generally speak- 
ing, arbitrary form, 


In fact, the method of search is based on the exploration of the neighborhood of the working points, 
X1, .. «» Xn, by means of small tentative increments (steps) of all parameters in turn, and an extrapolation of 
further motion on the basis of the exploratory data, The results of the explorations are the magnitudes of the 
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Fig. 1, Block schematic of the 12-channel optimizer 










partial derivatives (in practice, the first partial derivatives) which permits the making of a prognosis of the most 
effectivé temporal motion of the working step (directed along grad Q), The search process consists of an alterna 
tion of exploratory operations (gradient determinations) and working steps in the direction of the gradient, where 
the order of alternation depends on the search method employed. The instrument consists of two basic blocks, 
namely, an operational (analog) block and a controlling (digital) one. 








The purpose of the operational block is to carry out the determination of the magnitudes of the partial 
derivatives, to remember them and to carry the working steps for each independent variable, the steps being 
proportional to the corresponding partial derivatives, In accordance with this purpose, the operational block has 
12 parallel channels, at whose common input is located the AQ device which determines the first difference of 
two values of the input signal, Each channel consists of a memory device Mj(i = 1,.. ., 12) and an integrating 
device I;, In addition, there is a determination circuit which gives the output voltage of integrator lj as the 
tentative increment 5x at the proper moment of time. The outputs of the integrator channels are the object's 
inputs, The outputs of the object of optimization are the function Q whose extremum is sought, the limiting 
functions Hy, Hg, Hs, and Hy and the total limiting function Hy [1]. 













The channel elements, as well as each channel with the A.Q block, are joined by means of electronic keys. 






The operation of the apparatus is carried out with a definite succession of closings at certain times (cycle 
times) of the corresponding keys (or groups of keys) in the operational block. With this, the given order of switch- 
ing corresponds with a definite order of the operations to be carried out on the quantities Q, Hj and xj, i.e., is 
determined by the search algorithm, Thus, for example, closure of key ky in the AQ block leads to the remember- 
ing, in this block's memory cells, the value of Q(x, . . ., X29). Closure of k; (kp open) with simultaneous ap~ — 
plication of 5x; leads to the remembering in Mj of the i'th partial derivative, Closure of all k, keys leads to 
the change of xj by the quantity Ax;, proportional to 8Q/@x;. The sequency of operations is described in detail 
in [1, 3] and also in [4, 5]. 





































Fig, 2, Photograph of the 12-channel optimizer, 


To implement the necessary order of switching (operations), the controlling block is used, several block 
schematics of which are given in [3], In the variant described here, the controlling block is constructed on the 
basis of a refinement of the “two-loop ring” circuit [3] where, as the ring elements, triggers (flip-flops) were 
chosen, The circuit is a two-loop closed chain (register) of triggers (Fig. 1, near the bottom), each of which 

is controlled by the corresponding key (or group of keys) in the operational block, If the excited state is induced 
in one register trigger then, at each following cycle (clock time), the excited state will transfer to the succeed- 
ing trigger. Thus, it will circulate in a closed ring either in the large loop of Ty, Ty, . . ., Tig, Os, Thi, depend- 
ing on the state of the keys kgr and Kg,. 


Since each register trigger is controlled by the corresponding keys in the channel (trigger excitation cor- 
responds to key closure), the propagation of the excited state around the large loop, Ty, Ty, . . ., Ty, implements 
the operation of determining the gradient, and localization of the excitation in the small loop Ty, O,, Tj, im- 
plements the motion in the direction of the gradient [4]. 


The circuit refinement is in closing the ring, not immediately, but by means of a reset of the entire register 
to the null state at the end of a working step, The null state corresponds to the nonexcited state of all the trig- 
gers except the blocking trigger, Tj}, which is excited in preparing the register for the next cycle, With such a 
circuit, a random error in the ring (for example, the simultaneous excitation of several triggers ) would disrupt 
optimizer operation only during one cycle of register operation, 


The controlling block also includes additional devices for processing the signals necessary for controlling 
keys Kor and k,,, operating as a transfer contact and also kop and Kop. Among these devices are; 1) the 
maximum indication trigger Tri, operating together with gradient trigger Tgr, which is controlled by keys 


and ke and 2) limitation detecting tigger T,,_,, operating jointly with trigger TQH, controlled by keys koy 


and ko);. 

The operation of T,,; amounts to the reaction to the sign of the increment of the working step as a func- 
tion of time, Since, with the passage through a particular extremum, the sign of AQ changes, this signal is used 
for sending commands to the gradient determination after several working steps, as this is required in operating 
by the method of steepest descent, To avoid false operations during the time of gradient determination, when any 
sign of the derivatives is possible, T,,,; is blocked during that time by means of trigger q)) in the controlling block, 
The operation of trigger Tq} amounts to the switching of the inputs to the AQ block if the boundary Hy = 0 is 
crossed in some direction or another, whereby there is simultaneously sent a signal to T,., which switches in the 
gradient determination, Trigger T;,,,, which gives the signals switching the instrument input from Q to Hy and 


‘conversely, is also blocked during the time of gradient determination to avoid false switchings during application 


of tentative steps. More details on circuit operation are described in [5]. 


To start the optimizer, it is necessary to establish initial conditions, i,e,, to establish the output voltages 
at the outputs of integrating links lj. For this, by special tumbler switches, the register is places in the initial 
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(null) state for which only trigger Tj) is excited, After the initial conditons are established, the blocking of the 
initial state is removed, and the excited state, under the action of the cycle pulses, generated by the cycle multi- 
vibrator, begins to circulate successively from one trigger to another in an order which depends on the sign of 

AQ at the time of the working step. 


Figure 2 shows a photograph of the optimizer with its supply block, The optimizer chassis has the dimen- 
sions 670 by 560 by 330 mm, its face being the mounting plate for the individual blocks, which are interconnected 
by plug-in connectors, When the optimizer is in operation, its state at any moment, and its operation as a whole, 
are easily inspected by means of the light signals on the mnemonic schematic, 


Description of Certain Blocks of the Scheme 





1, The AQ Determination Block, The AQ block determines the accuracy with which the extremum co- 
ordinates are found, Its circuit therefore calls for special attention. In the variant described here, the circuit 
shown in Fig. 3 was used, It is basically a cathode follower with a polystyrene capacitor, C = 1 microfarad, con- 
nected in series with the grid circuit, The cathode follower's grid can be shorted to ground via electronic key 
Ky. The capacitor's other end can be connected to one of the inputs (Q or H) by means of keys kop; and Kone 
operating like a transfer contact, The cathode follower'soutput is applied to the scale amplifier. When kg is 
closed, the capacitor charges exactly up to the input voltage and, after switching out of ko, the cathode follower's 
output voltage equals the difference between the voltage at that given moment and the voltage which existed at 
the moment when k, was switched out [6]. 




















Fig. 3. Circuit for the AQ determination block, 


Of fundamental influence on the accuracy 
of the AQ determination is the discharge of 
capacitor C by the grid current of the cathode 
follower and the back current of the diode in 
key k. If special measures are not taken, the 
drift of the output voltage is about 1 volt/ min, 
Better results are obtained by dropping the fil- 
ament voltage on the vacuum-tube diodes to 
5.5 to 5 volts, and by stabilizing the operating 
mode of the cathode follower by means of a 
series connected vacuum-tube, By means of 
these measures, it was easy to obtain drifts not 
in excess of 0,1 volt/ minute on vacuum-tube 
6NZP, 
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Fig. 4. Block schematic of the channel for variable x;. 
































2. The Block of the Channel for Variable 
x;. (Fig. 4), The memory link is a cathode fol- 
lower with a capacitance at its input. The ca- 
pacitance can be charged from the output voltage of the AQ block if electronic key k;, at the link's input, is 
open, In the closed position of the key, the voltage on the capacitor remains constant, 
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The integrator is a three-stage dc amplifier 
Qvolts Gvolts shunted by a capactitive feedback path. The input 
1 voltage is applied via electronic key k,. The ten- 
tative increment 6x is applied by means of a jump 
of the controlling voltage on key k; through a smail 
a 80 capacitor connected to the integrator's input, where 
the magnitude of the jump can be controlled by a 
limiting diode, The capacitors in the grid of the 
memory stage and in the feedback path of the integ- 
60 60 rator are polystyrene, With this, the drift rate of the 
integrator's output voltage,when the capacitance in 
the feedback path is 0.25 microfarads , does not ex- 
ceed 1 to 2 volts/ minute [6]. 
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- The initial conditions are established by switch- 
isos ing the integrators to the aperiodic link mode of 
Ae | operation, The entire system contains 12 channel 
1 z= search 20 blocks, 

} ae : Bridge-type electronic keys are used in the 
Ee l puff optimizer, since they allow the number of vacuum- 
tube diodes to be halved and also possess low forward 
a 0 b 0 impedance [6). 











Fig. 5. Diagrams of searches on specimen objects, Automatic Optimizer Test Results 





The preliminary tests of the optimizer were 
carried out on a reference object which formed an integral part of the instrument, This was.a finder of the sum 
of the moduli of twelve input voltages, i.e., its output was a function of 12 variables with one exact extremum 
for Vin, =... = Uing = 9. Typical search diagrams are shown on Fig. 5, a and b, where the initial steps of the 
search, in accordance with the method of steepest descent, are shown, 
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Fig. 6. Block schematic of the set-up for automatic approxima 
tion, I is start, II is stop, II is clear, IV is the synchronozing 
pulses, V is to the control EMU-6, VI is the function supplier, 
VII is control NP-1, 





When close to extremal points, the instrument transfers over to the gradient method, thus permitting one 
to obtain smaller steady-state oscillations of the function whose extrema are being found, The diagram 5,b 
shows a case when the optimizer channels were not tuned accurately, but where the transition to the gradient 
method permitted a smaller error in the determination of Qmin to be obtained than in the case when the method 
of steepest descent was used, 


Fig. 7. Photographs of successive steps in the search for an extremal. 


As a test problem for the automatic optimizer when operating with an electronic model, we chose the 
problem of the automatic approximation of a given function of time, In other words, the problem was to construct 
by means of the optimizer, a function of time which, in the sense of some criterion of nearness, would coincide 
most closely with a given function of time, As the criterion, the integral criterion 


g=\|U—U'|dt, 


. 


was chosen, where A is the interval of comparison of the two curves, U is the given function of time, generated 
in the model and U®* is the approximating curve constructed by means of the optimizer, i.e, U* = f(X,,.. «,%a2)- 





Fig. 7, c,d, Photographs of successive steps in the search for an extremal, 


Of course, the condition for the best coincidence of the curves in this sence is that Q(x,, . . ., X4g) = min, 


Figure 6 gives the block schematic of the set-up for the automatic approximation, 


For the generation of the approximating function by means of the optimizer we used a nonlinear function 
generator, developed by L. N, Fitsner, which played the role of a switch, switching the optimizer's output voltage 
X1, . . «» Xz to the input of an auxiliary operational amplifier. The auxiliary amplifier could be used either in 
the scaling or the integrating mode, depending on which one obtained either a stepwise or a linear approxima- 
tion, Both functions were compared with each other and their difference was applied to the integrator via the 
block for modulus determination, All the operational amplifiers are blocks of electronic model EMU-6 and are 
controlled by synchronizing the ring of single flip-flop oscillators "start", "stop" and “clear”, This ring gives 







































signals to the model, to the optimizer and to the memory stage between the model's output and the optimizer's 
input, The presence of the memory stage eliminates large oscillations of the optimizer's input voltage resulting 
from the periodization of the solution, It is convenient to observe the approximation process on an electronic 
oscillograph screen, particularly if the screen has afterglow, For this, both curves are shown on one screen, by 
means of an electronic commutator, and one can observe how the approximating curve approximates to the con- 
tinually repeated given curve, Figure 7 shows the successive stages in the approach of the two curves, the figure 
being obtained by photographing the oscillograph screen, Figure 7,a and b shows the case of a stepwise approxi- 
mation to exponents and parabolas, while Fig. 7,c and d shows the case of a linear approximation to these same 
curves, As in clear from the photographs of the curves, the maximum error provided by the approximation (in 
the case of the linear approximation) was about 4,5% in the case of the exponential curve and about 8% in the 
case of the parabola, 


Thus, the test results show that the process of convergence exists for any form of the function Q. 


Certain Questions Regarding Search Method 





The most important characteristic of the search process is the time taken to find the solution (certainly, 
with the necessary accuracy), particularly when there are limitations present and the search process is markedly 
retarded due to the zig-zag nature of the motion of the representative point about the boundaries [1]. In this 
connection, a variation of the search methods in the presence of limitations is worthy of attention [7]. Here, the 
extremum is sought, not for the function Q, but for the function Q + aH», where ao issome proportionality factor, 
Certainly, in domains free of limitations, Hy = 0, Let one boundary Hj be transgressed and immediately there 
is a counteracting signal, whereby only the component of Q normal to the boundary is compensated, If a is so 
chosen that the normal component is exactly compensated, then the condition for an extremum lying on the 
boundary will be that grad (Q} aH,) equal zero, It is easily seen that motion tangent to the boundary occurs 
both in entering a boundary zone and in leaving it, thanks to which the search process is significantly accelerated, 
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Fig. 8. Circuit for the automatic control of the coefficient a, 
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It should be stated that, since Q nor Hy is known beforehand, neither is the coefficient a, However, the 
choice of a can be made during the search process itself, Indeed, in a boundary region, the magnitude of the 
increment AH, during the time of a working step must always be less than zero (for the case of limitations of 
the type Hy <0) and, generally speaking, the magnitude is small, in order to avoid large oscillations at a nom- 
inal minimum point, Thus, control of a may be implemented by means of the simple circuit shown in Fig. 8. 


Here a summer of the quantities Q and Hy is located at the optimizer's input, where Hy is applied to the 
summer's input via a block of the controlled coefficient, for example, the product block, The input which con- 
trols the block of the coefficient to be controlled is applied from the integrator, at whose input the quantity 
4H + € is applied via the two series-connected keys k, and kjjm. These keys are controlled by the single flip- 
flop step oscillator and, respectively, the limitation indication trigger, 


The circuit operates only when kg is closed, and attempts to maintain 4H =— €, where ¢€ is a given size 
of step for H, which guarantees, on the one hand, that the boundary zone will be left and, on the other, that there 
will be small oscillations on the boundary at the nominal extremum, point, 
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Fig. 9. The cases of “false extrema” (a function of two variables, 
1) Excitation; 2) direction of gradient; 3) line of gradient; 4) tentative 
excitation; 5) direction of search, 
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Very important also is the question of the class of functions which can be optimized by the instrument, 
There has already been mentioned in the literature of cases of “false extrema" [8], Figure 9 shows these cases 
of “ridges” and saddles” for the simplest case of a function of two independent variables, The function is given 
by means of lines of equal function values, In both cases, the determination of the components of the gradient 
by means of providing tentative steps along the coordinate axes leads to zero partial derivatives at the “ridge” 
points and the “saddle “vertices [8]. However, an extremal value is not found in either case, 


These cases differ from one another by the character of the stability of the positions found, namely, in the 
of a “saddle”, the position of the working point will be unstable, and a disturbance of its position, which can be 
given in order to test the validity of the extremum found, leads to a departure from the “saddle"point towards 
one of the “humps”, 


In the case of a "ridge", providing a disturbance after a stable position on the “ridge” is reached, leads to 
the attainment of another point of this same “ridge”. This circumstance can be used for the indirect determina- 
tion of the direction and magnitude of the gradient, Indeed, since the first stable point may be defined by the 
radius vector ry, and the second by the radius vector rz, the direction of the ridge line is given by the vector 
Ar = r,—ry, where the direction of the gradient along this line will be indicated by the sign of the difference 
Q:- Q,. 


Thus, the automatic optimizer can carry out verification of the extrema found without making use of the 
Box-W ilson method [8] with its enormously cumbersome procedure of applying 2" combinations of tentative in- 
crements, To do this, it is only necessary, when the working point has reached a position where 0Q/0 x, = 
= = 0Q/dxp = 0, to remember the coordinates of the extremum thus obtained, and then to translate the 
working point in some arbitrary direction which does not coincide with the direction of approach to the extremal 
point, Thereafter, the regular search process is carried out, as a result of which a stable state will again be at~- 
tained, If the coordinates of the new stable state coincide with the remembered coordinates, the search should 
be discontinued, since the extremum will be a true one, If the coordinates do not coincide, it is necessary to 
Carry out a search along the line of a vector Ar = ryz—r, and in the direction given by the sign of AQ = Q,;~Q,, 
and when a new stable state is reached by means of the regular search procedure, the applying of a disturbance 
should again be repeated, The procedure just described should be followed until the coordinates of two or three 
successively found stable positions coincide, 


Figure 10 gives one variant of the circuit for realizing the algorithm just described, 


In addition to the ordinary block schematic of the optimizer, shown in the right side of Fig. 10 in simpli- 
fied form, there are n + 1 blocks, A, for determining the differences of the coordinates x; of two stable states 
and the values of the function Q which correspond to them, The outputs of the A blocks, via keys kj, are con- 
nected to a common input which, either directly or via an inverter, can be connected to the output of block AQ 
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Fig. 10, Simplified block schematic of the optimizer with 
a device for verifying the validity of the extrema found, 





through one of the verification keys Kye. Keys kye and k ye operate like a transfer contact, Keys kj at the blocks’ 
outputs are controlled in parallel with the keys k; in the channels, The inclusion, or noninclusion, of the inverter 
is determined by the sign of AQ = Q,—Q;. The keys which control the memories in the A blocks are controlled 
by single flip-flop oscillator O', which is added to the controlling block, 


In addition, a two-position commutator, shown in the left part of Fig. 10, is added to the controlling block, 
the purpose of this device being the commutation of pulses which pass through key k, sending these pulses either 
through the chain of single flip-flop oscillators 0", and O4 and thereafter to T},), or directly to T),) with a simul- 
taneous reset of T,.. Pulses can be applied to the commutator only in case a stable state has been reached, The 
attainment of an extremum (false of valid) is fixed by the. extremum indication trigger Tjpq which is connected 
to the output of the AQ block. 


Verification of an extremum's validity occurs in the following manner, 


After the stable state has been reached (AQ is very small during the working step), Tj,g switches keys k 
and & in such fashion that the pulse from the single flip-flop step oscillator is applied now, not to T},) (more 
exactly, not to the reset of the register to its zero state), but to output I of the commutator, After the pulse is 
applied, the commutator prepares for the following pulse from output I. 


The impulse from output I excites the single flip-flop oscillator for remembering the coordinates of the 
extremum reached, O',, and then the single flip-flop disturbance oscillator Og, which switches some previously 
established voltages to the inputs of the channel integrators, After the disturbances are applied, T},; and then 
Tg are excited, where T, resets the output voltage of the AQ block to zero and, since the magnitude of the most 
recent working step was small, the determination of the gradient occurs (Cf,, the description of the optimizer's 
block diagram), After the determination of all the partial derivatives, the working step follows and a large 
magnitude of AQ appears since, as a result of the disturbance, the working point moved away from the stable 
state, Therefore, Tind opens key k and closes k, after which the search process described at the beginning of 
this paper transpires, 


When the new stable state is reached, Tj,4 again operates and, and after the next working step, a pulse is 
applied to verification trigger Tye. This trigger, by means of keys kye and kye, switches the channel inputs from 
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the output of the AQ block to the outputs of the A blocks, either directly or via the inverter, depending on the 
sign of AQ = Q;—Q;. Now, when an excitation circulates through the elements of the T,— Tyg register, the Ax, 
from the outputs of the A blocks are applied to memory links Mj with the proper signs, since keys k, andk, are 
switched simultaneously, After Ax, is remembered, Tye is reset to its initial state, which corresponds to the 
ordinary search procedure. During the working step, Ting also switches key k to the ordinary ring. After this, 
the ordinary search process proceeds in the direction defined by the remembered Ax,. 


After the stable state is reached, disturbances will again be generated, and the processes of verifying an 
extremum are repeated, 


It is easily seen that the procedure described is also valid for the “saddle” case since, after reaching a 
vertex of one of the neighboring “humps”, the search is carried out in the direction of the line joining the “saddle” 
point with the “hump” point. If the “hump” point is a true extremum, the search process returns to the same 
point and Ax, = 0, {.e,, the search process terminates, 


SUMMARY 


1, The instrument described may be used for solving by direct methods a broad class of variation problems, 
both with limitations and without them, The simplest examples were given in the paper, 


2. The number of independent variables (12, in our case) can, if necessary, be greatly increased without 
undue complication of the instrument, 


The author wishes to thank A. A. Fel'dbaum for this scientific direction, and A, V. Kalinina, A. B. Shubin 
and I, N, Bocharov for taking part in the adjustment and testing of the instrument. 
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Descriptions are given of apparatus used in constructing pneumatic switch- 
ing circuits, Methods are demonstrated of realizing by pneumatic means the 
elementary logical functions and a one-cycle delay function, It is also shown 
that any finite automation canbe constructed by means of pneumoautomata, 




































The circuits are given for several pneumatic devices which implement 
analog~-to-digital conversions. Examples are given of the application of 
pneumatic switching circuits in industrial pneumoautomatic devices, 


INTRODUCTION 





Contemporary systems for the automatic control of productive processes contain, as their most important 
constituent parts, various computing devices and devices which implement certain logical functions (switching 
circuits), 


In industrial automation, a tremendously extensive use is made of pneumatic apparatus; this is related both 
to the superior features of pneumatic executive mechanisms and to the fact that this apparatus is simple, reliable 
and, by its nature, explosion- and fireproof. There are many factories where all the automation is effected by 
pneumatic means, If, in the control systems of such factories, it is necessary to use computers or switching cir- 
cuits and, if these latter need not, under actual production conditions, be very fast-acting, then it would certainly 
be eminently reasonable to use pneumatic variations of these devices and circuits, 


In gauging the domains where it would be advantageous to use pneumatic computing devices and switching 
circuits, one must also take into account that, in the first place, in chemistry, petroleum processing, thermo- 
technics, thermoenergetics, metallurgy and many other branches of industry one can use computing devices and 
switching circuits with speeds of operation which are fully attainable by pneumatic means and, in the second 
place, that pneumatic apparatus is in principle efficient for high-temperature operation, 


In the Laboratory for Pneumohydraulic Automation of the IAT AN SSSR, work has been carried out on 
the construction of pneumatic analog computers and pneumatic switching circuits, The present paper shows how 
switching circuitcs can be realized by pneumoautomatic means®* , 


At the base of the pneumatic switching circuits to be described here is the same relay-contact principle 
which is used in constructing electromechanical switching circuits [1], In the pneumatic variant this principle 
also leads to devices which consist of circuits with relays connected in them, but in this case both the circuit 
itself and the relays are pneumatic, The signals on which the switching circuits operate also have a pneumatic 
nature here, In the circuits under consideration, the signal levels will have the following correspondence with 
the pressure of the compressed air; a signal with the symbol 1 corresponds to a pressure of 1 atmosphere, a signal 
with symbol 0 to a pressure of 0 atmospheres, 


* The results of this work were present on Septembe: 26, 1958 at a conference on the application of digital (dis- 


crete) automatic systems, Also participating in this work were L, I, Rozonoer, A. A. Tagevskaya, N. V. Grishko, 
T. K, Morozova, P. A. Lebedev and M., M. Sharafetdinov. 
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Fig. 1, Apparatus used in constructing pneumatic switching circuits, 
a) A pneumatic relay; b) an uncontrolled pneumatic valve; c) a con- 
trolled pneumatic valve, 


With the method chosen for forming the signal levels in the pneumatic switching circuits being considered, 
the pneumatic circuits which form them will contain, in addition to pneumatic relay, also uncontrolled pneu- 
matic impedances, 


A pneumatic relay (Fig. 1,a), just as a relay used in electromechanical devices, consists of two basic 
parts: a reacting organ and contact (contacts), The pneumatic contact is realized in the form of a nozzle-flap- 
per pair, This element, established in the pneumatic circuit, constilates a local impedance whose magnitude 
depends on the distance between the flapper and the nozzle and may be varied from zero, when the flapper is 
removed a sufficient distance from the nozzle, to infinity, when the flapper is moved right next to the nozzle, 
If the nozzle, when controlled, can assume only these two positions, then the nozzle-flapper element will im- 
plement the function of a contact, i,e., will either completely open the pneumatic circuit or close it moreover, 
without the introduction of additional impedances in the circuit, in practice, In the future, when a nozzle-flap- 
per element is used as a contact, we shall say that the pneumatic contact is closed when the flapper is remote 
from the nozzle, and that it is opened when the flapper completely covers the nozzle, 


The reacting organ of the pneumatic relay consists of three plane membranes, They are connected axially 
to a rigid center whose movements are limited by two bearings, while the membrance peripheries are fixed in the 
frame, In the dead-end chambers which separate the membranes the command pressures are applied, these caus- 
ing the appearance of some total force which acts on the rigid center along its axis, This force shifts the rigid 
center towards the bearing located on the side to which the force is directed, 















The connection of the pneumatic relay's reacting organ with the pneu- 

matic contacts appertaining to it is implemented by using one of the bear- 

> Seemeepeee ings which limit the motions of the rigid center as the nozzle and by using 
the butt end of the rigid center which corresponds to this bearing as the flap- 


c , per of the pneumatic contact, Thus, the reacting organ of the pneumatic 

6 relay and its pneumatic contact turn out to be combined in the physical de- 
sign, 

e —4+—_ 


The uncontrolled pneumatic impedances are placed in the circuits 
either at those places where the pneumatic relay connects with the 1-atmos- 
g —D-— phere supply trunk, or where the pneumatic relay communicates with the 

atmosphere. If, with this, formation of a signal occurs in a portion located 

—() between a pneumatic impedance and the remainder of the pneumatic circuit, 
then closure (opening) of the pneumatic relay in the first case leads to the 
Fig. 2, a) A pneumatic appearance of the signal 0(1) and, in the second case, to the appearance of 
circuit path; b) the sup- the signal 1(0). 
plying end of a pneuma- 
tic circuit path; c) the 
“grounded” end of a pneu- 
matic circuit path; d) an 
uncontrolled pneumatic 
valve; e) a controlled 
pneumatic valve; f) a 
pneumatic contact; g)a 
pneumatic relay's react- 
ing organ; h) a pneu- 
matic capacitance. 


It turns out, in connection with such a method of using uncontrolled 
pneumatic impedances, that their magnitude has an effect on the speed of 
response of the circuit, on the air discharge in them and on the deviations 
of the levels of the formed signals from their nominal values (0 or 1 atmos- 
phere). However, if it is taken into account that the circuits to be realized 
can operate quite satisfactorily within a broad range of the characteristics 
just noted, and that only the magnitudes of the pneumatic impedances affect 
circuit operation, it then becomes understandable that this element is of 
modest importance, In practice, the uncontrolled pneumatic impedances 
are implemented in the form of capillary segments (Fig. 1,b). 


Pneumatic relays, consisting of reacting organs and pneumatic contacts 
and constant pneumatic impedances constitute the set of basic apparatus nec~ 
cessary for the construction of pneumatic switching circuits, In a device for pneumatic automation which con- 
tain switching circuits, it frequently happens that controlled pneumatic impedances are also used, The circuit 
for such an element is shown in Fig. 1,c. 


The connection of the individual pneumatic relays and pneumatic impedances in one device (the working 
mounting of the circuit) is done on special plates. The pneumatic relays and the controlled pneumatic impedances 
are mounted on the surfaces of these plates, and all communications (pneumatic circuits) occur inside them, 


Naturally, a pneumatic relay-contact apparatus, being analogous to the corresponding electrical apparatus, 
will possess the same functional capabilities, It is shown below in what ways pneumatic relay~contact apparatus 
should be used to create devices which bear the generic designation of “finite automata” (Cf., for example, [2]). 
We first give methods of constructing devices which realize the elementary logical operations, On the basis of 
these devices we can construct circuits in which operations are implemented which are defined by logical func- 
tions of any complexity, the so-called primitive circuits. We then give the circuit for the most important link 
in finite automata, namely, the delay element. It is then shown how the so-called nonprimitive circuits may be 
constructed on the basis of the devices in which the elementary logical operations necessary for the construction 
of primitive circuits are implemented and of the delay elements, 


In addition, we provide certain information on the construction of analog-to-digital devices which are based 
on the same apparatus, In concluding this work, we provide examples of the use of pneumatic relay-contact ap- 
paratus in the circuits of certain devices which are of important practical value. 


The basic content of the material given below is contained in the circuit diagrams which, for economy of 
space, are accompanied by short explanations only. 


Notation and Symbolism 





For the graphic representation of pneumatic switching circuits we use the notation enumerated on Fig. 2. 
In addition to these designations we also use a special symbolism which represents the method of connecting the 
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TABLE 1 pneumatic relay's reacting organ in the circuit, Essentially, this symbolism 
is a variant of the corresponding symbolism used in the theory of electrical 
Pa relay-contact circuits [1]. . ; 








x=P, The fundamental variable quantity in pneumatic circuits is the magni~- 
tude of the compressed air pressure, In those cases, when this variable can 
TABLE 2 assume only two values (0 and 1 atmosphere) it will be denoted by P with 

subscripts, If it can assume any value in the interval from 0 to 1 atmosphere, 

P, it will be denoted by P with subscripts, A second important variable is the 

x = B, quantity x which defines the state of a pneumatic contact; we shall say that 

x = 1 if the pneumatic contact is closed (the flapper is remote from the nozzle) 
and x = 0 if the pneumatic contact is opened (the flapper covers the nozzle), 
Variable x is the pneumatic relay's “output”, its value being determined by 
state of the pneumatic relay's "input", i.e., by the values of the command pressures P, and Ph applied to the 
pneumatic relay'‘s reacting organ. In the general case the relationship x = x(P ,; P},) is established. Its concrete- 
form is determined by the special characteristics of the command pressures and by the method of connecting the 
pneumatic relay, In the circuits which we shall be considering in the sequel, pneumatic relays, realizing five 
forms of this relationship, will be considered. 























Pneumo-con- Reacting ¢rgan ee matic 
tact notation representation realy 


a ~ hy. P P 


connection 














Fig. 3, The relationships realized in pneumatic relays and the cor- 
responding symbolic representations of the reacting organs and pneu- 
matic contacts; a) x = Sg(P,~Pp) = [P,-Pp}k b)x= Pa; c)x= 

= Py; d) x=P/ +P; e)x =P, Py. 


a) x = Sg(P,—Pp)=[Pa—Pp]} A pneumatic relay which realizes this relationship will be called a receiv~ 
ing relay; it is connected in accordance with circuit A, shown on Fig, 3, The command pressures P g and P), 
applied to it may vary within the interval of 0 to 1 atmosphere, 


b) x =P, (Table 1), To realize this function, the pneumatic relay is connected in accordance with the 
same circuit A, but has one input B which assumes values of either 0 or 1; the pressure supplied to the second 
input remains invariant (fixed) and equals approximately 0.5 atmospheres, 




























TABLE 3 c) x = P, (Table 2), In this case the pneumatic relay is 








again connected in accordance with circuit A and has one in- 
P, 0 | 1 0} 1 put, Ph; with this, the second input remains “fixed” at the value 
Pb 0 rte 1 , 1 of 0.5 atmospheres (P, = 0.5 atmospheres), 
— — -_ 








oa eye oe d) x= P, > P, (Table 3). Here, both inputs are "free", 
x= P,>Ph 0 1;0;]9 i.e., both P, and P), are variables which assume the values of 0 
and 1, To realize this function, a spring is provided to the re- 
































ele , ‘ n “a acting organ of the pneumatic relay, as shown on Fig, 3 in cir- 
TABLE 4 ech 
e) x =P, ~ Ph (Table 3). This case differes from the 
Py 0 1 previous one by the location where the spring is places (Fig. 3, 
SS circuit B,). 
P= P, 1 0 





The symbolic representations of the reacting organs and 
P=Py 0 1 the contact notations shown in Fig. 3 correspond to these five 
methods of employing a pneumatic relay, The method used 
to notate the contacts is exhaustive in relation to the informa- 
tiou on the methods of connecting the corresponding pneumatic relay, 




















Realization of the Elementary Logical Operations, Primitive Circuits, 





As follows from the propositional calculus [3], any primitive circuit may be constructed with the use of 
only some one elementary logical opegation, for example, the “Sheffer stoke" operation, However, it is more 
convenient, even theoretically, to construct logical functions and primitive circuits by the use of some two ele- 
mentary logical operations, for example, negation and conjunction, or negation and disjunction, Even more 
convenient is the widely used system constructed with the use of three elementary logical operations, namely, 
negation, conjunction and disjunction, Naturally, all other still broader systems, while providing no new theo- 
retical capabilities, may turn out to be more convenient still in practice, 


TABLE 5 
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P =P, V P2VPs\.. VP, 














We show below now, by pneumatic means, one realizes the operations of “negation”, "replication", "Shef- 
fer stroke", “conjunction”, “disjunction","implication", “equivalence” and “logical sum". Such a list of element- 
ary logical operations is known to be sufficient for the construction of logical functions of any complexity and, 
consequently, any primitive circuit may be constructed on the basis of devices which realize these operations, 


All these operations may be divided into three groups: 
1) One-input (Table 4) — negation (P = P,) and replication (P = P,). 


2) n-input (Table 5) — Sheffer stroke (P = P;/P)/ .../Pn), conjunction (P = PrAP2f/\, . ./\Pn) and dis- 
junction (P = P}\/P,\/.. VV P,). 
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TABLE 6 3) Two-input (Table 6) — implication (P = Py - Ps), equivalence 
| (P = P; ~ P,) and logical sum (P = PyVP3). 
Pt : | dh! 8 The circuits for the devices which realize the operations of the 
Ps 0 | 0;1)1 first group are shown on Fig, 4, All these circuits are constructed with 
the use of pneumatic relays of the type x = P, or x = Py, (Fig, 3), each 

P= P,; +P, 1 O;1)}1 type of relay having its corresponding variant of circuit construction, 
eA | ciel o =e The circuits for the devices,realizing the operations of the second group, 

ibe | are shown in Fig, 5, They are constructed with pneumatic relays of the 
P=P,V P, 0 1;1;0 same types. For each operation here, two variants of circuits construction 
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Fig. 4. Circuits of devices re- 
alizing the elementary logical 
functions of one independent 
variable (Table 4), a) “Nega- 
tion", "no", P = Py; b) “repli- 
cation", “yes, P = P;. 





are given: with series and with parallel connections of the pneumatic 
contacts, 


The circuits of the devices which realize the operations of the 
third group are shown in Fig. 6, These circuits are constructed with 


pneumatic relays of the types x = Pg — Pj, or x = Py, +P, (Fig. 3), Here 
also, each type of relay has its corresponding variant of circuit construc- 


On the last three figures, under the circuit which realizes a par- 
ticular operation, there is given the symbolism for this device, 


Realization of the Delay Operation, Delay Lines, 
Ring Counter Circuits, 








Figure 7 represents the circuit for an element which implements 
the operation of a one-cycle delay (the delay time is one nominal time 
unit), This element has two independent inputs P, and P; and one out- 
put P. Input P, has a special meaning: the successive appearance on it 
of ones and zeroes defines the division of physical time into cycles (into 
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Fig. 5. Circuits of devices realizing the elementary logical functions of n in- 
dependent variables (Table 5). a) The "Sheffer stroke", P = Py/P)/ .. ./P,; 
b) “conjunction”, "and", P= PsA P,A...A\Pp; ¢) "disjunction", "or", P = 


= P,VP,V.. V Po. 


units of nominal time), just as is shown on Fig, 7. For further convenience at each moment of time correspond- 
ing to some unit of nominal time we distinguish two successive moments; one, t,,, corresponds to the end of 


the n'th cycle, and the other, ta(n + 1), to the beginning of the (N + 1)'st cycle, In an idealized circuit these 
two moments coincide, We shall assume that tan = thin 4 1) + 0. The delay operation realized by the element 


under consideration relates its input and output by the relationship P = D'P, which, with the notation we have 
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Fig. 6, Circuits of devices realizing the elementary log- 
ical functions of two independent variables (Table 6), 

a) “Implication” P = P,P ,; b) "equivalence", P =P, ~ 
~ Py; c) "logical sum", P = PyVP3. 


assumed, has the following content: P(t) = P;(t— 0) if t = tan (with this, t — 0 = tn —,) ) and P(t) = const if 
tan <te< tan + 1): 

The graphs of Fig. 7 illustrate the operation of the delay element which, in fact, amounts to the "remem- 
bering” of the input quantity during the time corresponding to one cycle, Each cycle may be divided into two 
parts; in one of them, P, = 0 and, in the other, Py = 1, In correspondence with this, the delay element also con- 
sists of two parts: one of them “remembers” during the time when P; = 0, and the other when P; = 1. A series 


connecttion of n delay elements in one chain forms a delay line (Fig. 8) in which the following operation is 
implemented: 


P, = D'P,, Py = D* P,,..., P, +, = DP. 


of the 


By connecting the output, P 
last element of the delay line to the input, P;, 
of the first element, one obtains a ring counter 
circuit (Fig. 9). The graphs of Figs. 8 and 9 
show the operation of a delay line and a ring 
counter circuit, 


n+ 1° 





Nonprimitive Circuits 





The addition of delay elements to a system 

—~ Mmm fm of devices which realize elementary logical func- 
at thy :taz be tions qualitatively changes the capabilities of these 
batt ay i -- a Sele devices, Without the use of delay elements in 


4 : “ooo them, only primitive circuits can be constructed, 


3 :nal unit whereas when used with delay elements, they 
P ae || —ae give the possibility of constructing nonprimitive 


if ae — circuits of any degree of complexity, i.e., any 


A =. De finite automata, This is borne out by the follow- 
ing. It is known from the theory that any finite 
4 automaton can be constructed from universal 
elements, As such an element, one may take 


the model of a nerve cell (neuron) suggested by 
McCulloch and Pitts [2]. 








Fig. 7, A one-cycle delay element, P = D’P;, 
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Fig. 8. A delay line. Fig. 9. A ring counter circuit, P, = D'P,, Ps = 
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a 
Aes 
‘ Fig. 10, The neuron of McCulloch and Pitts, a) Simplified schematic 
of a neuron with n exciting and m inhibiting inputs and a firing (operat- 
ing) threshold of h; b) a circuit realizing a neuron with n = 3 and h =2, 
based on “no”, “and", “or and “delay"elements, 
m 
c- The fundamental property of a Pitts-McCullogh neuron (Fig. 10,a) is chat, at time t+ 1, there appears at 
ese its output P the symbol 1 (output excited) if and only if not less than h of its n exciting inputs was excited at time 
t and none of its m inhibiting inputs was excited at that time, 
, Figure 10,b shows a neuron with n = 3 and h = 2, its logical formula having the form: 
P = D{{(Py A Ps) V(P2 A Ps) V (Pi A Ps) A Ps4i A +++ A Ps4m): 
* This formula is constructed from the operations of “no” (“not”), “or”, “and” and “delay”, and the logical 


formula for any other neuron is also constructed with these same operations, Thus the capability of realizing 
these four operations is a condition for the possibility of constructing any neuron and, consequently, any finite 
automaton, It was shown previously how these four operations are implemented by pneumatic means, This means 
that it is possible to produce any neuron and any finite automaton by pneumatic means, 


As an example of the realization of a nonprimitive circuit, Fig, 11 gives the circuit for a pneumatic pulse 
divider (a pneumatic trigger). This simplest, but practically important, circuit consists of a “no” element whose 


output is connected back to its input via “delay” element, The 
logical formula corresponding to this circuit has the form: P = p’P, 


The graphs given on Fig, 11 illustrate the circuit's operation, 


A Device for Analog-to-Digital Conversion, 





In the previous sections we have considered the question of constructing 
finite automata by pneumatic means, HowevVer, digital (discrete) devices are 
widely used in automation which do not fall in the class of finite automata, 
Such are the devices,to be considered now, which pneumatically transform 
analog (continuoug quantities to digital ones. 


Fig. 11. A pulse divider, The basic element in these devices is the receiving pneumatic relay 


Py = D'P,, P, = Py. of the type x = Sg(P, ~ Py) (Fig. 3). 














Fig. 12. An interval indicator: P},», = Sg(Paz:—P), Phs = Sg(Pas — 
—P), .» +» Pon = Sg(Pan —P)» Pi = Phas P2 = Ph2 > Phs» Ps = Phs > 





























Fig. 13. A pneumatic isodrome regulator with automatic changeover, 
1) Pneumatic regulating block; 2) contains valves a and 6 with their vents; 
3) interval indicator, 









The receiving pneumatic relay, used as an - 
analog-to-digital converter, allows various receiving 
devices to be constructed when used in conjunction 
with the previously considered logical element, For 
‘example, Fig. 12 gives the circuit of a so-called 
interval indicator for continuous (analog) quantities, 
In it, the pressures P,; subdivide the entire range of 
variation of the input pressure P into n intervals; to 
each pressure P,; in the device there corresponds a 
definite section with its own output P; (we note that 


Pal = %- 

















Thus, the input quantity P is put into corres- 
pondence with some i’th interval by the output P;. 
In operation of the device, if the input pressure as- 
sumes a value in the i'th interval, then the P; which 
corresponds to this interval will have a zero value, 
and all the remaining sections will have output pres- 
sures of unity, 
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Fig. 15, Automatic optimizer of an object with several controlling organs, 
1) Object with n controlling organs; 2) the automatic optimizer's memory 
block and comparison block; 3) a pulse divider; 4) integrator; 5) block 
of. integrator welves with vents, 6) block of power amplifiers with vents; 
7) is the pulse counter;: 8) ring counter circuit, 


Some Applications of Digital Pneumatic Switching Devices 








A Pneumatic Isodrome Regulator with Automatic Changeover, The ordinary controlling block of a unified 
system aggregate (USA) realizes a proportional-integral control law, and has two adjusting organs, One organ 























































establishes the regulator gain, the other, the isodrome time constant, The.tunings are maintained invariant during 
regulator operation, In many cases, for example, when the object's properties depend to a large degree on the 
load, and the device operates over a wide range of loads, difficulties arise in choosing a regulator tuning which 
would be equally correct for all modes of object operation, As was shown in [4], it is meaningful here to use a 
regulating block with an attachment to provide automatic variation of the regulator tuning when some external 
condition or another changes, This attachment (Fig. 13) consists of an interval indicator and a block of valves 
and vents, When the regulator operates with such an attachement, there is automatic variation of the regulator 
tuning whenever the object's mode of operation changes (as measured by pressure P at the input of the interval 
indicator) from one interval of values to another, 


A Pneumatic Semiproportional Regulator, In industrial automata, the most widely used regulators are those 
which implement regulation laws (linear) which are constructed of some combination of proportional, integral 
and differential components, USA controlling blocks and lead blocks allow the realization of any of these laws, 
In [5], attention was turned to the advantages of using, in some cases, regulators which implement the (nonlinear) 
semiproportional control law suggested by V, Ferner, It is shown on Fig. 14, how by joining to a standard pneu> 
matic USA controlling block.an attachement containing two receiving pneumatic relays, one logical element 
(equivalence) and two vents, a proportional-integral block can be transformed to a semiproportional one. 





An Automatic Pneumatic Optimizer of an Object with Several Controlling Organs. The existing pneumatic 
automatic aptimizers for ordinary industrial use [6] are designed for operation on objects in which the quantity to 
be optimized depends on the position of only one controlling organ, However, by the use of a special attachment, 
this same optimizer can also be used in those cases when the quantity to be optimized depends on the position 
of several of the object's controlling organs (Fig. 15). In the case, optimization is carried out by successively 
moving the controlling organs to positions which provide the attainment of a nominal extremum of the quantity 
to be optimized. It is also assumed that three searching motions suffice for the attainment of such a position by 
each controlling organ, The circuit shown in Fig. 15 is constructed for such a method of optimizing. It includes 
the elements of an automatic optimizer of an object with one controlling organ: a sampled-data measurer of 
distance from an extremum (a memory block and a comparison block), a pulse divider (a stepping switch) and an 
integrator (a constant pressure drop block). In addition, the circuit contains a pulse counter which marks the ap- 
pearance of every third pulse, a ring counting circuit with a number of sections equal to the number of control- 
ling organs of the object, and a block of valves and repeaters with vents. 
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PUNCHED-CARD METHOD OF SYNTHESIZING SEQUENTIAL 


SYSTEMS OF MULTIPOSITION RELAYS 


V.i. Shestakov 


(Moscow) 


A punched-card method of synthesizing sequential systems of multiposition 
relays is presented, This method is valid for the synthesis both of autonomous 
and nonautonomous switching circuits, and is a generalization of the punched-card 
method of synthesizing sequential systems of two-position relays. 


INTRODUCTION 


A vector algebraic method [1] of synthesizing sequential systems of r-position relays has been presented, 
this method being a generalization of the analogous method of synthesizing sequential systems of two-position 
relays [2-4]. By using this method,one canobtain functions which formulate the necessary and sufficient condi- 
tions of operation of all the r-position relays in the system to be synthesized and the sequence of external stimuli 
to this system, In the general case, these functions are not Boolean functions, as in the case of systems of two- 
position relays, but are operations of a more general algebra which has been put forward by D, Webb *. 


As digital automation grows, an ever-increasing value inheres in automatic systems constructed of relays 
with more than two stable states, Of the greatest interest are the relay systems constructed of three-position 
relays, since such systems can contain the least number of elements, Work is currently going on in the designing 
and building of automatic computers which use the ternary number system, Th‘s number system is the most 
economical, as concerns the equipment necessary for computers, and therefore the work in designing computers, 
operating in the ternary system,is very pressing and very promising. 


Another pressing problem in computer technology is the development of computers using ten-position-relays 
for computing directly in the decimal system without translation of the numbers to some other number system, 
The solution of this problem is particularly important for the creation of high-speed computers designed for fin- 
ancial computations and for computations in various branches of the national economy, 


The tendency to make ever greater use of multiposition relays is observed, not only in computer technology, 
but also in other domains of digital automation. The development of discrete automata makes the problem of 
synthesizing sequential systems of multiposition relays an urgent one, 


To solve this problem, one must use a new algebraic apparatus, namely, the algebra of multivalued logic 
and in particular, the algebra of Webb. The use of this method entails a large number of computations in the 
synthesis of sequential systems of multiposition relays. As the number of independent signals and relays in the 
system to be synthesized increases, the number of these computations grows very rapidly, much more rapidly than 
the number of analogous computations necessary for synthesizing seqeuential systems of two-position relays, 
Therefore, the requirement of mechanizing these computations arises, 


This paper presents one of the simplest methods of mechanizing the vector algebraic method of synthesizing 
sequential systems ofr-position relays, our method being based on the use of special punched cards, This method 





* The Appendix to this paper contains a short sketch of Webb's algebra, 











is a generalization of the punched-card method [6] which was initially developed for the synthesis of sequential 
systems of two-position relays, It is assumed that the reader is acquainted with work [6]. 


The Vector Algebraic Method of Synthesizing Sequential Systems of r-Position 
Relays 





The vector-algebraic method of synthesizing sequential systems of r-position relays has already been de> 
scribed in [1], and will therefore not be presented in detail here, We shall present only such information regard- 
ing this method,which is necessary for the presentation of its implementation by means of special punched cards, 


For brevity, we shall henceforth refer to a sequential system of r-position relays simply as a system, If all 
relays of the given system have identical delays, and if the magnitude of the delay for relay operate equals the 
magnitude of the delay for relay release, and also if all controlled and controlling parameters of the system 
change virtually instantaneously and synchronously, then the change of state of this system can be described by 
the equation 


y(¢+ t) = f(x(?), y()), (1) 
where f is some single-valued vector function and r is the delay time, If we assume that r = 1 (which, without 
any loss of generality, may always be done), then the previous equation can be replaced by the following: 

y(i+1)=f(x(/), y(), (2) 


where j=0,1,2,... 


Equations (1) and (2) differ from the corresponding equations describing systems of two-position relays only 
in that the vectors x and y which enter into these equations are no longer binary, but are r-place vectors, where 
r > 2, and f denotes, not a Boolean function, but an arbitrary function of Webb's algebra, 


The set of sequences { x(j)} and{y(j)} is called complete if, among the terms of these sequences, there are 
contained all r™ yalues xq of vector x and all r" values yg of vector y, 


a=0,1,...,7r™—1; B=0, 1, ..., m—1. 


It follows, for the single-valuedness of function f , that only such sets of sequences {x(j)} and {y(j)} are 
valid for the synthesis of the system described by Eq, (2), which satisfy condition (B) in work [6] as applied to 
r-place vectors, 


The algorithm for synthesizing the system described by Eq. (2) is given by the formula 


N Ng-1 
1(x, y) = Di" Dd} Pajce () Pax (y) ¥ (7 + 4), (3) 
sx=1 j=0 


where is the symbol for the generalization of the Boolean sum, defined in Webb's algebra [5] by the formula 


1:011, N ss the number of pairs of sequences {x(j)} and {-y(j)} in the given complete set of them, N, is the number 
of states in the s"th sequence, and p,(x) and pg(y) are generalized constituent units, 


The generalized Boolean sum of two r-place variables u and v is defined by the formula 


u\/ v= max(u, v). (4) 
For binary variables, this operation coincides with the operation of Boolean addition 
The generalized constituent units po(x) are defined by the following conditions: 
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[, 7 eS 
Pa “10, eo w+ Re. (6) 


The generalized constituents Pg (y ) are defined analogously. 


The values of Xq and yg of the variable r-place vectors x andy will be numbered by integers a and 6 
by the same rules as in work [1], i.e., by the rules expressed by the formulas 


m n 
a= >) a ri—1, B = >) by rk—1 (6) 
i=] k=1 
where a; and by are the i'th and k*th components of the vectors 


Ka Mt ile, 5 «2, Oke reo Gene ¥p = [0,,..:, Op, .- +, Ouls (7) 
respectively, 


The generalized constituents p,(x) and pg (y) are expressed in terms of the components of the vectors 
Xa» % Yp and y by the following formulas: 


Pa (X) = lI bx, aj? Pp(y) = ll by, by (8) 
i=1 


k=1 


where 5 yj aj and 6 yk,bk #re the well-known Kronecker deltas, i,e,, functions defined by the conditions 


5 owes (9) 
potiaD Uo @ 926. 


For an autonomous system, i.e, a system described by the equation 
y(7+1)=f(y() 
all the formulas given above are simplified in an obvious manner, 


Punched Cards for Synthesizing Autonomous Systems 


J 





The vector algebraic method described above for the synthesis of autonomous systems of any r-position 
relays is a natural generalization of the similar method for the particular case of two-position relays, and there- 
fore, the punched cards intended to mechanize this method in the general case when = 2 can be made very little 
different from the punched cards for the synthesis of autonomous systems of two-position relays, In the general 
case, as in the previously described [6] particular case of two-position relays, for each value y of a variable n- 
dimensional vector y there is prepared an individual punched card, “punched card 6 * and, consequently, ther 
must always be prepared r" different “punched cards 6", 


Since it is almost impossible, with manual distribution of cards and visual reading of their numbers, to 
function with more than a thousand cards, the maximum number n, for a given r, may be defined by the inequality 
= 1000, from which it follows, in particular, than = 6 forr = 3 and n = 3 forr=10, In fact, 3° = 729 < 1000 < 

< a ie = 2187 and 10°< 1000 < 10*, i.e., the set of "8" cards for the synthesis of autonomous sytems must contain 
at most 729 punched cards in the case of three-position relays, and 1000 cards in the case of ten-position relays. 
Since both the number of cards in the set and the very construction of the cards must depend on r, we shall 
henceforth denote the punched card corresponding to the 6 ‘th value yg of the variable r-place vector y by the 
symbol 6 (r), and shall call this punched card “punched card 8(r)". We shall, in correspondence with this con- 
dition,call the punched cards described in work [6] for the synthesis of autonomous systems of two-position relays 
“punched cards 8 (2)". The sole difference, in the general case being considered here,is that each component by 
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of vector y = (by, ..., by, . . ., b,] will assume, not one of two values: 0 or 1, but one ofr values: 0,1,..., 
r —1, Thus, punched cards 6(r) can be prepared so that they differ from previously described [6] punched cards 
B (2) only for by = 2, 3,...,r-—1. 


By adhering to the * principle of least difference” of punched card 6 (r) from punched card 6 (2), we set 
aside the left half of punched card 8(r), as was previously done for cards 6 (2), for the perforation of the value 
of component b,, of vector yg , where the hole (perforation) corresponding to the value of the k*th component 
by is found in k*th row (counting down from above) of punched card 8(r). In the case when b, = 0, there will 
be no punch in the kth row. The value b; = 1, as in the case 6(2), corresponds to a rectangular punch at the 
edge of the k’th row, which, certainly, must not extend the entire width of the left half of punched card B(r), 
since it is necessary to leave space for the perforations corresponding to the values of bk = 2, 3,...,r—1. For 
these values of b,,, there is made a single punch in the k"th row of cards 6 (r) in the column corresponding to the 
given value of b;,, i.e., the punch is made at the intersection of the k*th row of punched card 6 (r) with the second 
column if by, = 2, or with the third column if b; = 3 or, finally, with the (r—1)'st column ifb, =r-1. The rect- 
angular hole, or punch, must have dimensions which allow one to see through the aperture all three digits of 
number 6 , of the punched card on which the given punched card is lying. By bearing in mind that the distance 
between punches equals the width of one digit, we obtain the fact that the width of the left half of a punched 
card must be not less than 4(r —1) times larger than the width of one digit. The right half of punched card 8 (r) 
must have exactly the same width, since on them must be placed r — 1 columns of three-place decimal numbers, 
each of which is the serial number of the given punched card 8(r). 





iS9 159 «159 159 159 159 159 159 159 
ro 159 «159 «159 159 ‘159 159 159 159 159 
(] 59 9 159 159 159 159 159 159 159 








a: 





123456789 
WI 23456789 
0 SI 23456789 
()s9 23456789 


b 








Fig. 1 


Figure 1, a shows punched card 6(10) = 159, prepared exactly in accordance with the method just described, 


The width of the punched cards can be decreased by a factor of one and a half to two, if their construction 
is changed somewhat, In order to obtain such a decrease in punched card width, it is necessary to do the follow- 
ing: 


1) To make the same rectangular punches as in the punched cards of the former construction not only for 
by, = 1, but also for by = 2,...,r—1. 


2) For by = 2,...,r—1, im the left half of the punched card, to punch an opening of one digit width at 
the intersection of the k'th row and the corresponding column, 


3) In the right portion of punched card, to leave the ordinal number 8 of this punched card in just the 
first (left-most) column, and fill the remaining columns, II, III and finally (r — 1), with the corresponding num- 
bers 2,3,.... (f—1). 


Figure 1,b shows punched card 6(10) = 159, prepared by this method. 


Punched cards 6(3) prepared by the method just described are approximately half the size of those prepared 
by the previous method, Figure 2 shows the punched cards for 8 (3): 


0 (3), 1 (3), 2 (3), 3(3), 4(3),..., 25 (3), 26 (3), 


prepared by the second of the methods described above, 


1460 


















































02 7 12 Lan 22 32 
02 12 22 S 32 
02 12 22 32 
02 12 22 $2 
02 12 22 32 
02 12 22 32 
a 42 25 2 0 62 
42 O 22 O #2 
42 isk D #2 O 22 
42 25 2 26 2 
42 25 2 26 2 
42) 252 26 2 
Fig, 2 


The 6 (3) punched cards omitted from Fig. 2 are easily reproduced by using a table which gives the first 
27 values yg of the six-place ternary vector y. 


Punched cards 6(r) are used for the synthesis of autonomous systems of r-position relays analogously to the 
way punched cards 8(2) are used for the synthesis of autonomous two-position relay systems, The punched ~card 
method of synthesizing the latter was described in detail in [6]. Therefore, it is not necessary to present the 
method of using punched cards 6(r) in the general form, We limit ourselves her to the consideration of one ex- 
ample of the use of punched cards 6(3) for the synthesis of an autonomous three~position relay system, 


Example, Let it be required to construct an autonomous system of three-position relays y;, y; and ys, such 
that they operate in the following sequence, (You Vir Yar + + +» Yas» Yas)» i.€., 80 that all the 27 possible states of 
the system of these three-position relays repeat periodically with period 27, and follow one another in the natural 
sequence given in the table, In other words, let it be required, for a given periodic sequence of values { yajcs) 
of the three-place ternary vector y which describes the state of the desired system, to find the function f, which 


satisfies Eq. (2a). For this, we take punched cards 6(3), represented on Fig. 2, and place them one after the other 
in the sequence 


0(3), 1(3), 2(3), 3(3),..-, By (3), Byy, (3), -- -» 25 (3), 26 (3), 03) 


such that the right half of each preceding card 6 ;(3) is covered by the left (perforated) half of the succeeding 
card 8; 4 ,(3). Figure 3 shows the initial and final fragments of this placement of punched cards 6,(3). The 
reader can easily regenerate the omitted middle portion by use of the table since, for the case considered, we 
have By =j for each value of j(j = 0,1,..., 26). 


a pal, a a Ba 





eeooco co 2° 
Ne NN NS & 























Fig, 3, 
Punched cards 8 (3) shown in Figs, 2 and 3 are useful. for the synthesis of autonomous system:which do not 


contain more than 6 relays Vj, Vz, ..., Vg. In particular, they are Valid for the synthesis of the autonomous 
system of our example, in which the system contained only three relays, In this case, the desired function f(y) 
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contains only three components, f ;(y), fy), and f(y). To determine these components, we use only the first 
three rows of punched cards 6)(3), placed as shown on Fig. 3, whereby each component f(y) is determined 

from the corresponding k'th row of these punched cards in the following manner, Instead of the serial number 

8, of punched card 6 ;(3), which is visible in the k'th row through the rectangular punch in the left half of punched 
card 8 j + 3(3) which covers the right half of punched card 6j(3), we write the designation pg j of the correspond- 
ing generalized constituent, where the symbol pgj is multiplied by 2 if, in punched card B; + 1(3), there is a 
punch in the second column of its right side through which the digit 2 is visible, By taking the genralized Boolean 
sum of all the generalized constituents thus obtained, both with and without the coefficient 2, we obtain the 
component fx(y ). 


By carrying out the procedure just described for all three rows of punched cards 8,(3), placed as shown on 
Fig. 3, we obtain the following definitive expression for the components of the function f(y) sought: 


11(¥) = (Po V Ps V Po \/ Po V Pi2 V Pts V Pas \V Par V Pra) VV 
V 2(P1 V Pa V Pz V Pio V Pas V Pris V Pio V Por V Pes): 
12(¥) =(P2 V Pa\V Pa V Pu V Pre V Pas V Pao VV Par V P22) VV 
V 2(Ps VV Po V Pr V Pras VV Pas VV Pie V Pas \/ Pos V Pos) 
Is(¥) = (Ps V Po V Pio V Pur V Piz V Pis V Pisa V Pas V Pris) V 
V 2(Piz V Pis V Pio V Poo V Par V P22 V Pos V Pos V P25): 
The components f;(y), f(y) and f,(y ) are the necessary and sufficient conditions for the operation of relays 


Yy, Yq, and Y3, respectively. It is easily verified that the expressions found for these components guarantee the 
operation of the system of these relays exactly as in the cyclical sequence given above. 












































i 
B by, | be | dy “is che bi, | be | bs brn 
NO ey 
0 () 0 0 eee 0 14 2 1 1 0 
1 1 0 0 owe 0 15 0 2 | 0 
2 2 0 0 oes 0 16 1 2 1 0 
3 0 4 0 ein 0 17 2 2 1 0 
4 1 1 0 vy 0 18 0 0 2 0 
5 2 1 0 ate 0) 19 4 0 2 > 
6 0 2 0 in 0 20 2 0 2 0 
7 1 2 0 ene 0 21 0 1 2 0 
8 2 2 0 oi 0 22 1 1 2 0 
9 0 0 1 = 0 23 2 1 2 0: 
10 4 0 1 0 24 0 2 2 0: 
11 2 0 1 0 25 1 2 a 0: 
42 0 1 1 0 26 2 2 2 & 
13 1 1 1 0 





In case the given complete set does not contain just one sequence, as in the example just considered, but 
several sequences, one should carry out the procedure described for each of these sequences, and then carry out 
Boolean addition of the components f;,(y) found for each of these sequences, The expressions thus obtained would 
be, in this case, the components f,(y) of the vector function f (y) sought, which described the autonomous system 
to be synthesized. 


The punched cards 6(r) are used for the synthesis of autonomous r-position relay systems analogously to 
the way punched cards 6(3) were used in the example considered, In the general case, instead of the serial 
number 6 ; of punched card 6;(r), which is visible through marginal rectangular punch in the k'th row of punched 
card Bj 4 s(t) covering the right half of punched card 8 j(r), one writes the symbol pg; of the corresponding gener- 
alized consituent Pg j¥)> where this symbol is multiplied by the number which is visible through the punch of 
the same punched card 6 j + 1(1) if, to be sure, such a punch exists, If no such punch exists, then only the symbol 
Pgy 1s written. If, in the k’th row of punched card 6 j + 1(1), there is'neither a marginal punch nor a hole, then, 
naturally, the symbol pg; will be lacking for this row. 
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By taking the generalized Boolean sum of all the generalized constituents thus obtained, both with and with- 
out the coefficients 2, 3, . . ., r, we obtain the k"th component of the vector function fl? ry) defined by the p "th 
given sequence. The generalized Boolean sum of the k'th components found for all the sequences entering into 


the given complete set will be the k'th component f ,(y)of the vector function f (y) sought, describing the auton- 
omous r-position relay system to be synthesized, 


Punched Cards for the Synthesis of Nonautonomous Sequential Systems 





As for the previously described particular case of the synthesis of two-position relays [6], for the synthesis 
of nonautonomous sequential r-position relay systems it is necessary to have cards of two forms— a cards and 
6(r) punched cards, For any value of r, the a cards are identical with the analogous cards as described in [6]. 
Cards 6(r) are different for different r, and are somewhat different from the corresponding punched cards intended 
for the synthesis of autonomous systems, In the sequel, we shall denote these latter cards by the symbol 6(r)},. 
Punched cards 8(r) differ from the corresponding punched cards 6(r), in the same way as punched cards 6(2) dif- 
fer from the corresponding punched cards 6(2),. The difference consists in this, that the perforated portion of 
punched cards 8(r) is extended by adding a column of punches which duplicate the rectangular apertures on the 
left margins of the punched cards, and the right portion of the punched cards are extended by an empty field 
along their right margins. The width of this empty field on the 8(r) punched cards must equal the width of the 


a cards, since it is intended precisely for the placement of the a cards in the implementation of the synthesis 
of nonautonomous sequential systems, 



































02 = Ls fe Sao 22 
02 12 z2 
02 12 22 
02 12 22 
02 12 22 
02 12 22 
32 CC] 22 OC) #2 
Lj 37 OC) a2 00) #2 
32 OC) az OC) #2 
32 25 2 26 2 
32 252 26 2 
32 25 2 26 2 
Fig, 4 
Figure 4 shows the punched cards 
B (3) == 0(3), 1(3), 2(3), 3(3),..., 25(3), 26 (3). 


Just as in the case of two-position relays, the numbers 6 (in the column to the right of the column of duplicate 


apertures) on cards 6(r) denote the natural constitutents pg (y), and the serial number borne by the card denotes 
the a-constituent p_(x). 


Punched cards A(r) and the a cards are used for the synthesis of nonautonomous r-position relay systems in 
exactly the same manner as in the previously described case [6] of such systems with r= 2, Therefore, we will 
not describe again here this method of using punched cards 8(r) and the a cards in the general case of an arbitrary 
I, We will limit ourselves to a discussion of Fig. 5, which shows the set of punched cards 


Bs (3) = 3, a = 5, Bypi(3) = 7, Oj41 = 12, Bj42(3) = 15, 





which is the fragment of the distribution of punched cards B,(3) and a cards in the sequence: 
Bo(3), %o, Pi(S), a1,---, Bi(3), %j, Bits, (3), X42, Bro (3),---, 
corresponding to the given (for synthesis) sequence of values of the ternary vectors y(j), X(J): 


YBor Mays Yrs Mayr +++ YB;) Kas YAj44) Masia YBs499 se. 
By reading, through the apertures of the punched cards shown on Fig, 5, the numbers 6 (3) and a), we ob~ 
tain the following fragment of the generalized Boolean sum: 


V Ps(X) Ps (y)¥ 2 V Piz (X) P2 (Y) Yas V - 


This expression is analogous to the fragment of the Boolean sum 


+ Ps (X) Ps (¥) ¥2 + Pre (X) Po (¥Y) is +» 


obtained from Fig. 9 of work [6], and differs from it only in the set of components, The first three components 
of the Boolean sum just given (i) are identical, and equal the Boolean sum of the products of constituent units 


Ps (X) Ps(¥) + Pra (X) Pr (Y)s 


while the first three components of the corresponding generalized Boolean sum equal the following generalized 
sum of products of generalized constituent units 


: [2] Ps (x) psy), 
J 2ps (x) Ps (y) V 2Pi2 (x) Pr (y); 





4 Cc QC 7ACZ P12 (X) P2(Y), 


(2) 


respectively, This difference is a consequence of the dif- 
ference in the make-up of the components of the single-place 
binary and ternary vectorsy, and yys5. If these vectors are 
binary, then y; = [1, 1, 1, 0) and yys = [1, 1, 1, 1]. If they are 
ternary, then y7 = [1, 1, 0, 0) and yy, = [0, 2, 1, 0]. 

















Fig. 5 
The difference just cited between the make-ups of the 

components of the binary and ternary vectors y; and y4, is easily observed by comparing the perforations in the 
punched cards 7(2), 15(2) and 7(3), 15(3), represented on Fig. 9 of work [6] and on Fig. 5 of the present paper. 


We note in conclusion that the 8(r) punched cards can also be used instead of the 6(r), punched cards in 
synthesizing autonomous systems, With this, ofcourse, the « cards are not necessary, The duplicating apertures 
and the empty fields on punched cards 6(r) will not then be used, and should be given no further notice, 


With manual positioning of punched cards 8(r) and of the ao cards, and with visual reading of the digits 
through the apertures in punched cards §(r), the punched card method described can be practicably used for the 
synthesis of relay systems with a total number of independent and dependent variables which does not exceed 
six, for r = 3, and three for r = 10. With mechanization of the punched-card method by means of punched card 
(EAM) machines, the total number of relays and independent variables can be doubled and, seemingly, even 
tripled, 


APPENDIX 
A Short Sketch of the Algebra of D. Webb 





Ina paper [5) published in 1936, D, Webb described an algebra, suggested by him, of an n-valued logic, 
this being a generalization of the two-valued Boolean algebra generally used as the classical propositional calculus. 































Using the author's somewhat simplified notation, we will here denote by the letters p, q and r variables 
which can assume the values 0, 1, ..., n, but we shall retain the numeration of the formulas used in the paper 
of Webb which was cited earlier. 


The fundamental operation in Webbian algebra is the operation defined by the equation: 


1.01 p/qg = 1-+ max (p, q), 





where the sign + (plus) denotes addition modulo n. In what follows, the signs + (plus) and — (minus) will always 
be used as the signs for adition and subtraction modulo n, so that we shall not bother writing the usual symbol 
(mod n) after each equation in which these signs appear. 


The operation p/q, defined by Equation 1.01, is a generalization of the well-known Sheffer stroke operation, 
All other operations of Webbian algebra are defined in terms of this operation: 





1.02 po-=p, p'ti—=pp' (i =0, 4,..., n—4) 
1.03 pV q=(piq)"" 

1.04 Np= p™—P-1 

1.05 p-q=N(Npy Nq) 

1.06 pdq=(p-q)™*” 

1.07 p=q=(p9)-(9>pP) 

1.08 pxq=p™'y qr". 


The operation p\/q and p- q are generalizations of the operations of Boolean addition and multiplication 
respectively, In Webbian algebra, as in Boolean algebra, one sometimes writes simply pq, instead of p-q. 


The operation Np is a generalization of the one's complement operation in Boolean algebra, or the opera- 
tion of negation in classical propositional calculus, 


The operation pq and p = q are generalizations, respectively, of the operations of implication and equiv? 
alence in classical propositional calculus, 


We now give some of the theorems proved in work [5]: 





p= p+h, 
PV q=max(p, q), 
Np=n—i—p 

p-¢ =min (p, 9), 


hig 0, if pa=0 
1.45 PX4=\max (p, 9), if pq0 


— 
iD whrs— 













Operation p x q is distributive with respect to the operation of generalized Boolean addition, t.e., the fol- 
lowing theorem 1,25 is valid 


PX(qVr)=(PX9)V (PX?r). 





By virtue of the commutativity and associativity (proven in [5) of the operations p\/q, p-q, we can agree 
to write the expressions 


PViaVr), p- (+r) mw px (¢xXr) 


without parentheses, {.e., to introduce the definitions: 






jpPVqaVr=PVQ@Vr), 
1.09 par p-(q-r), 
pxqxr=px(¢qxr). 






We introduce, by definitions, the following abbreviated notation: 





1.0100 Xy a poz, gi tgit! gel 
n—1 
1.041- > e,=mVnV...VqV...Ve__, 
i==0 
n—1 n—1 n—1 


and analogously for > >» pe 2 Xie j 
k=0 


i=0 j=0 


k» Where x; and x;, 4, , , ,, k are any elementary propositons of 


n-valued logic, 


The symbol oes introduced by definition 1,011 ts, for n = 2, the symbol for the Boolean sum (or disjunction, 
to use the terminology of the classical propositional calculus), In the general case, when n is arbitrary, this sym- 
bol will be called the symbol for the generalized Boolean sum, 


According to the definition 1,010, for n = 2, i.e., in Boolean algebra, there exist the two functions X, = 


= x' =X, X, =x" = x, where x is the one’s complement of x. 


n—] n—l n—l 
From theorems 1,27 — 1.33, proven in work [5] follows the formula >} >} jes >; (X,;xXVY;xXZ,)--1, 
i=0 j=0 k=0 


on the basis of which the functions X; x Y; x Z, can be called generalized constituent units. 


The following theorem was proven in [5]: 


FRR DD PG. eye ZL eM Ea), 
i=o j=0 k=0 


which asserts than any function F(x, y, . . ., 2) of n-valued logic can be expanded in terms of generalized con- 
stituent units or, as expressed by Webb, can be given in the normal form: 


F (i, fp eee Bh) X Xp X Vu_y X~-++ XK Ay_e- 
The normal form of a function F(x, y, . . ., 2) of n-valued logic is a generalization of the complete normal 


disjunctive form of a proposition in the classical propositional calculus, 


A number of purely logical questions were also considered in work [5], in particular, the relationship of 
Webbian n-valued logic to other systems of multivalued logic and to the logical scheme presented in the Prin- 
cipia Mathematica, 
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SHIFT REGISTERS WITH LOGICAL FEEDBACK AND THEIR 
USE AS COUNTING AND CODING DEVICES 


A. N, Radchenko and V, I, Filippov 


(Leningrad) 


Methods are considered for increasing the efficiency of counting and cod- 
ing devices based on shift registers, The distinguishing feature of such devices 
is that the shift register is closed into a ring, whereby the closing link is a logical 
circuit,whose state is a function of the code introduced into the register, Methods 
are presented for constructing logical feedback circuits in order to obtain count- 
ing devices or coding devices of arbitrary capacity, and methods of reading out 
information are also given, 


In various devices for remote control and computing, use is frequently made of pulse counting circuits 
built of binary flip-flop circuits in the role of code combination transducers, However, in many cases, serious 
difficulties arise in the use of such schemes, This springs from the inevitable delay when signals are transmit~- 
ted from low-order bits to higher-order ones, these delays being cumulative and, depending on the number of 
bits, of a greater or a smaller total value. As a result, false combinations may occur during the transient responses, 
For example, the transition of an ordinary binary counter from the state 0111 to the next state, 1000, is accom- 
panied by a series of false states (shown in parentheses): 0111 -(0110 - 0100 + 0000) +1000, The most notice- 
able are the appearances of false code combinations in comparatively slow devices, for example, in relay-con- 
tact ones, although it is obvious that the relative lag, expressed as the number of cycles of lag to the high-order 
bits as compared to the low-order bits, will be approximately identical for both slow and fast counters, if each 
of them operates near the upper threshold of its speed of action, The lag of the high-order bits with respect to 
the low-order bits can reach 5 to 10 cycles, If the counter circuit is being used as a delay line, as occurs, for 
example, in dynamic flip-flops, then the discrepancy can attain an even greater value, 


In order to eliminate the cause of false combination formation, it is necessary that the transition to each 
successive fixed state be occasioned by the operation of the flip-flop for one bit only, This route is a well- 
known one today, It is based on the use of special counters which operate in correspondence with the sequences 
of a Gray code, However, this is not the only method for eliminating false intermediate combinations, The 
nonsimultaneity of flip-flop operations, which is the cause of false code combination appearance, can be avoided 
if one makes the transition from a serial flip-flop circuit to a parallel one, i.e, one in which the pulses being 
counted are applied to all bits simultaneously, 


A shift register can serve as the basis for the construction of such a counter circuit, Simultaneous action of 
all its cells (bits) is provided by the shift bus in which the pulses to be counted are introduced, The elements 
entering into a shift register may be quite diverse: electronic, semiconducting, magnetic, etc, The method 
presented below for constructing computing and coding circuits is applicable to any case, However, the most 
technological and economic advantage is gained by the use of shift registers of magnetic elements with rect~ 
angular hysteresis loops. A condition for this is that the usual disadvantage of such elements — the necessity of 
special means for reading out ‘nformation — is completely excluded-in the present case, this necessity also be- 
ing lacking when dynamic triggers (flip-flops) are used, With this, particularly for low counting speeds, there 
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there is a sharp drop in power requirements, since magnetic registers ordinarily do not have internal power supplies 
and require power only during the transient responses, As a result, the circuits constructed in accordance with 
the method described below are some two to three times as simple and economical as flip-flop circuits of equi- 
valent capactiy, An additional superiority of shift-register circuits is the absence of errors in dynamic modes of 
operation, 
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To convince ourselves that, with a shift register closed in a ring, one can realize more than n discrete 
states (where n is the number of register cells), we consider the circuit of Fig. 1. It contains four ferrite cores, 
Ty, Tg, Ts and T,, joining to a common circuit by means of diodes D; — D, and capacitors Cy, Cs, Cs and C,, 
The distinguishing feature of the circuit is the use of the logical “not” element in the feedback path, whose 
presence permits a doubling of the capacity of any analogous ferrite circuit, In the given case, the capacity is 
increased from four to eight, in a decade counter the circuit would permit counting up to 20, etc, 


Let us follow through the operation of the circuit, At the initial moment, let all cores be so magnetized 
that a pulse applied to the input does not induce the reversal of their magnetic states, This pulse is divided on 
resistor R and passed through the “not circuit, charging capacitor C,. Until the termination of the pulse C, 
cannot be discharged, since diode D, is cut off by the voltage of the input pulse which acts on the diode's cathode 
via the input winding of T,. After the termination of the pulse, the cut-off voltage falls and condenser Cy, in 
discharging through D,, the input winding of T, and the internal impedance of the pulse source, reverses the mag- 
netic polarity of core T,;, A second pulse to be counted returns core T, to its previous state, whereby C, is charged, 
Simultaneously, the second pulse passes through the “not” circuit and again charges condenser C;, When this 
pulse terminates, condenser C, and C, discharge through the windings of T, and T,, and reverse the magnetic 
polarities of these cores, The character of circuit operation does not change until all the cores have had their 
magnetic polarities reversed, which occurs at the termination of the fourth pulse to be counted, The operation 
cycles so far considered can be conventionally expressed by the following table: 


0000 — initial state 

1000 — after application of the first pulse 
1100 — after application of the second pulse 
1110 — after application of the third pulse 
1111 — after application of the fourth pulse, 


Here, a zero denotes a nonexcited core state, and a one denotes a state of reversed magnetic polarity. 


The character of circuit operation now changes: the pulses divided by resistor R will not be sent to the 
feedback path, since there will simultaneously arise pulses at the output winding of T, which will cut off the 
“not” circuit, As a result, condenser C, will be discharged at the fifth cycle, and the circuit will go to the state 
0111, Thus, circuit operation continues in accordance with the following table: 


0111 — fifth cycle 

0011 — sixth cycle 

0001 — seventh cycle 

0000 — eight cycle (the initial state). 


Thereafter, circuit operation repeats the previous sequence, 
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This circuit will cary out counting up to 8. A circuit of n cells will count up to 2n, 


The logical function in the feedback circuit essentially affects the counting factor of the circuit, Figure 2 
shows the same counter as in Fig, 1, but with adifferent feedback circuit, Here, a core with a rectangular hysteresis 
loop is used in the feedback path, Pulses applied from core T; take it to the "1" state, Pulses applied from T, 
return Ts to the "0" state, If pulses from T; and T, appear simultaneously, core T, does not change its state, since 
the corresponding magnetization polarity reversing windings act in opposition, Reversal of magnetization polarity 
is accompanied by a voltage on its output winding. This winding, in conjunction with diodes Dy, and D,,, forms 
a two-half-period rectifier circuit at whose output appear pulses of positive polarity, independently of the direc- 
tion in which the magnetic polarity of T,; was reversed, 


Let us consider circuit operation, with the assumption that, in the initial state, the first core isin the state 
"1" and the others are in state "0". The following table gives the sequence of circuit states: 


1) 1000 6) 0110 41) 1110 
2) 0100 7) 1044 42) 1414 
3) 0010 8) 0104 13) O111 
4) 1001 9) 1010 14) OO11 
5) 1100 410) 1101 15) 0001 
1) 100C etc, 


The first three states require no explanation. When the fourth pulse appears, core T; is returned to zero, 
thus occasioning the reversal of magnetic polarity of T, and charging of condensers C, and C, (via T,). At the 
termination of the fourth cycle, the cores are excited in the combination 1001, The fifth pulse returns cores T, 


and T, to the "zero" state, with which core T, in the feedback path is also returned to the zero state, This leads 
to the charging of C, and C3, which provides the appearance of the combination 1100, etc, 


Transmission of a one to the first cell occurs if one of the cores T; and T, is in state "1". No transmission 
of ones occurs if both these cores are in state "0" or in state "1". A necessary condition for the transmission of 
units is the alternate action on T; by Ts and T,, It can be verified from the table that this condition is satisfied. 
In fact, by writing the individual states of T, and T,, we get 


Ts 001001101011110 
T, 000200110101111 


The underlined ones transmit information to the feedback path, It is easily seen that they follow alternately 
in the upper and lower lines, Otherwise, creation of a feedback path with a core with a rectangular hysteresis 
loop would be impossible, 


The feedback circuit shown in Fig, 2 is sufficiently universal, and can be used to construct code counters of 
arbitrary capacity, For example, a circuit of seven binary cells with the same feedback circuit provides counting 
up to 127, 








The method considered for constructing counting circuits is general, t.e., for 





~~ any given number n of relay elements, it is always possible to construct a feedback 
path which will provide counting to a number m =< 2". The theoretical basis of 
we ~ this assertion is the proof of the existence of code rings [1]. 
100 1001 = 001 
ym, on “0 Code ring is the name given to a particular closed sequence of symbols, in the 
n-term segments of which are contained code combinations, where each of these 
Fa | S 3 1s combinations is contained once and only once, For example, the closed sequence 
7 yo ty 0011101, the code ring for a three-term binary code, contains in its segments all 
ww) Mi = possible three-place combinations; 000, 001, 911, 111, 110, 101, 010 and 100, The 








number of terms in a code ring, independently of the number of terms in each code, 


07 
“4 Mt \\\ 
. always equals the number of elements in the code ring. 


mi In a code ring, neighboring (intersecting) code combinations have n-1 common 
terms, Therefore, to obtain each successive code combination, it suffices to replace 
only one term in the previous combination, Precisely this operation is implemented 
by a register with a shift of combinations by one step: n-—1 terms are retained (shifted 
to more-significant orders) and, in place of the disappearing right-most term, the new term of the next combina- 
tion of the code ring appears on the left by means of the feedback path, Since the sequence of symbols in a code 
ring is always known in advance, one can always cetermine the previous or following terms of the ring for any n- 
term combination, It follows from this that the logical circuit of the feedback path can be constructed by realiz- 
ing a Boolean function which is given in the following form: the function's arguments are all the possible n-term 
code ring segments, and the values of the function are the code ring symbols which are preceded by the segment 
which is the current argument, 


Fig. 3 


Therefore, the method of setting up the logical feedback path for a counting circuit reduces to the following, 


1, Choice of the code ring corresponding to the conditions of the problem, The length of the code ring (the 
number of its terms) must equal the counting coefficient of the circuit to be developed (i.e., the number up to 
which it is supposed to count), A method for constructing such a ring is presented in [1]. Synthesizing code rings 
from the graphic representation of the set of code rings, as suggested in [2], can be recommended, Figure 3 shows 
such a set for four-term rings. To construct a code ring, one starts from some node and completes a circuit in the 
direction given by the arrows, returning to the original node, With this, a two-fold passage over any one segment 
is prohibited, By writing the sequence of symbols appearing during the passage, and eliminating the terms com- 
mon to each two neighboring scgments, one obtains a code ring. 


2. Setting up the operating conditions of the feedback circuit, By “operating conditions", we understand 
those states of the individual counter cells for which the logical feedback path is closed and for which a "bne” is 
written at the circuit's input. These conditions are written in the symbolism of logical algebra in a form con- 
venient for future transformation, 


3, Enumerating the unused states, i.e., code combinations which do not enter into the ring. 


4, Setting up the function of the feedback path, The function of the feedback path is defined as the logical 
sum of the operating conditions of this circuit, Moreover, we have the right, with the aim of simplifying the result, 
to add to this sum any expression from the list of unused states, if by so doing we ease further operations [3], 


5. Setting up the functional schematic of the shift register with logical feedback, On the basis of the feed- 
back path function just obtained, and using the rules of logical (Boolean) algebra [4], we set up the functional 
schematic of the counting device, 


6. Verification of the correctness of circuit operation, For this, we write the sequence of operating states 
and compare them with the sequence of codes of the selected code ring. 


We now elucidate what has been said with examples, 
Example 1, To develop a feedback circuit for a counter with counting coefficient of 10, 

1, We choose any code ring of the ——— possible four-term binary rings, for example 
—1001111000- 
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2. We set up the operating conditions: 
1) 1000 


2) 1100 
or 
3) 1110 


4) 0011 
5) 0001 
3. We enumerate the unused states: 
1) 1011 
2) 0101 
3) 1010 
4) 1101 or 
5) 0110 


6) 0000 


abcd 


abcd 
abcd 
abcd 


abcd 


abcd 
abcd 
abed 
abcd 
abcd 


abcd 


(2) 


(3) 


4. By writing conditions (2) and adding to their sum the expression abcd from (3), we obtain the following 


feedback-circuit function: 


) = abed + abed + abed +- abed +-abed +- abed = ad + abd, 
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(4) 


5. We now set up the feedback circuit from the wind- 
ing of the shift register elements, The first term of expres- 
sion (4), ad, can be reproduced by a contrary connection of the 
windings of cells A and D in series with a diode oriented away 
from A, It is obvious that there will be a pulse at the output 
'¢ of such a circuit only when there is a pulse in the winding of A 

, and no pulse in the winding of D, The second term, abd, can 
be reproduced analogously by connecting D counter to A andB, 
°° where the diode now must be oriented away from D, Both 
ae circuits are connected in parallel, The resulting circuit is 


shown in Fig. 4, a, The input, output and cycling windings 


Fig. 4 of the cores are not shown on the diagram, 


6, We verify the correctness of circuit operation, For this, we give the sequence of working states; 


1) 1000 
2) 1100 
3) 1110 
4) 1111 


5) 0111 


6) 
1) 
8) 
9) 
10) 


0011 
1001 
0100 
0010 
0001 


(5) 


We thus convince ourselves that the circuit has ten states, and that their sequence corresponds with the se- 


quence of codes in the selected code ring, 


Example 2, To develop a feedback circuit for a 20-cycle counter, Since 2 < 20< 2°, the first circuit 


must be a five-place one. 


1, We give ourselves a code ring 


—11110101001100010000 — 








We set up the conditions for transmission of a "1" by the feedback path: 
a) mandatory 1) 00001 6) 01010 

2) 00010 1) 10101 

3) 10001 ~=8) : 11010 (7) 
4) 00110 9) 11101 

5) 01001 


b tional 
a have 1) 01110 7) 11100 


2) 10111 += 8) - 11011 


3) 01011 9) 01101 

4) 00101 10) 10110 (8) 
5) 10010 11) 10000 
6) 11001 12) 11111 


We set up the feedback circuit function, writting the conditions, for convenience, in the following order 
(the optional conditions are given in parentheses): 4, (10), 6, 8, (5), 2, (4), 7, 9, (9), 5, (6), 3, 1. 


p=abcde+tabcde+abcede+abcede+abede+ 
4+abcde+abcde+abcde+abcecde--abcde-+ 
+abede+abcde+abede+abcde. 


After some uncomplicated transformations we get 
) = (b+ c) de + de (9) 


4. We set up the feedback circuit (shown on Fig. 4,b). 
5. We verify the correctness of circuit operation, 


The number of possible variants of the feedback circuit of a device of given capacity is large, since the 
number of nonidentical rings of a given length can be quite significant, For example, the number of complete 
rings (of length 2") is determined by the formula 

n-1_ 
p=2° 4 (10) 

By computing the number of variants of full-length code rings from formula (10), for example, for n = 5, 

we obtain p = 2048, 


Thus, the number of possible feedback circuits for a counter of fixed capacity is very large. Characteris- 
tically, among these variants are both very simple circuits and very complicated ones, No one has yet found 
a universal criterion for choosing the simplest feedback circuit, although it is obvious that the regularity of the 
code-ring construction immediately affects the complexity of the circuit, There is,as yet.no literature on the 
synthesis of code rings, Only two sources [2, 5] are known where attempts were made to synthesize code rings, 
where the method presented in the second source is far from exhaustive, 


The presentation, or read-out, of information deserves a brief mention. 


Parallel read-out of the codes in all orders of the register is automatically provided at each cycle of cir- 
cuit operation; each pulse applied to the register’s input simultaneously acts on all cells, transferring them to 
the “zero” state, With this, current pulses arise in the windings of those cores which contained “ones", these 
pulses supplying information to other blocks of the device, Parallel read-out of information {is characteristic of 
only some computing devices; in the majority of telemetering and remote control devices, the code must be 
transmitted along a communications line and must therefore be reproduced in the form of a temporal sequence. 
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of staggered pulses, In the given case, a serial code can be obtained by the use of a delay line, For this, the 
pulses from only one order are applied to the communications channel without delay; the pulses from the other 
orders are delayed by different time intervals, as a result of which the parallel code is transformed into a serial 
one, 


This method is difficult, however, in the case of large time intervals between code pulses, Therefore, it 
is more advantageous to carry out reading by introducing a series of pulses into the counter, 


For repeated read-out of the same information without disturbing the cores" magnetic states, one can use 
the circuit shown in Fig. 5. Here one uses a second feedback path containing an “and® circuit, while a “not” 
circuit is connected in the first feedback path, 


When there are no read-out pulses present (these pulses being applied to the “read-out” terminal), the 
"not" cell does not affect the character of circuit operation, i.e., the feedback signals are passed, In the read- 
out mode, the "not" circuit interrupts the basic feedback circuit, but the code, which we read out, is rewritten 
in the counter via the “and” circuit, 


LITERATURE CITED 


[1] A. N, Radchenko, Code Rings and Their Use in Contactless Coding Devicés [in Russian) Dissertation 
(1958). 


[2] N. G. De Bruijn,"A combinatorial problem,* Koninklijke Nederlandsche Akademie van Wetenschappen 
t, XLIX, 1946, 


[3] V. N. Roginskii, "Taking unused states into account in the synthesis of switching circuits,” [in Russian] 
Automation and Remote Control (USSR) 15, 3 (1954). 


[4] M, A, Gavrilov, Theory of Switching Circuits [in Russian} Izdatelstvo AN SSSR (1950). 


[5] Lippel and Epstein, "A method for obtaining complete digital coding chains,” Transactions IRE, vol. 
EC-6, No, 2, p. 121, 1957, 


Received January 2, 1959 








ON CALCULATING FUNCTION-GENERATING SHUNTED 
POTENTIOMETERS FOR CHEBYSHEV APPROXIMATIONS 


S. G, Kislitsyn and F, L, Litvin 
(Leningrad) 


A method is presented for the calculation of a shunted potentiometer 
intended to reproduce a nonlinear function, The potentiometer is chosen so as 
to provide the best (Chebyshev) approximation to the function to be reproduced, 
The design method presented is extended to potentiometers with finite and with 
infinite load impedances, The method presented in the paper is illustrated by 
a design example, 


Linear potentiometers with homogeneous rheostats and shunting impedances (Fig. 1) are widely used fcr the 
approximate reproduction of functions of one independent variable, By controlling the magnitude of the shunt 
impedances and the places at which they are connected, one can so arrange matters that the error has a Cheby- 
shev character, i.e,., so that the following expression is minimized: 





| A (x) |max = max | f(x) — P (x, rx) |, 


where f(x) is the function to be reproduced, P(x, r,) is the function actually reproduced by the potentiometer 
(the approximating function) and r, (k = 1, ...,m) are the parameters of the function-generating potentiometer, 
i.e,, the magnitudes of its shunt impedances, 


- Two cases of potentiometer design should be distinguished: 





a) for an infinite load impedance, the approximating function 


7, it vi Srl 4 =| | is piecewise-linear; 
b) for a finite load impedance, the problem of a piecewise- 


nonlinear approximating function arises. 


























1. The Piecewise-Linear Approximation 





Fig. 1 





The necessity for piecewise-linear approximation arises not 
only in the design of shunted potentiometers, but also in the case of 
designing devices with diodes [1]. To solve this problem, it is necessary to determine the parameters of the piece- 
wise-linear function and the boundaries of the pieces, For the minimum number of segments, the approximation 
must be a chebyshev one and the error is | A(x)| < €, where € is a previously given positive quantity, The solu- 
tions must encompass two types of functions: with a constant sign of the second derivative (no points of inflection), 
and with a variable sign of the second derivative in the interval of reproduction. 


Let there be given the function y = f(x) which must be approximated by the piecewise-linear function P(x) 
such that, in the interval (xp; xx], the error | f(x) — P(x) | < €, and the segments of approximation will be 

























as large as possible. With the error symetrically distributed on the positive and negative sides, the approximating 
function must lie within the band defined by the inequality (Fig. 2) 


{ (x) 4-€>P(2) > (z)—e. 


The best (Chebyshev) approximation in the first segment will be attained if the graph of the linear function 
P(x) is drawn through the point [x», f(x) + €] and is tangent to the curve f(x)— € ; the boundary, X;, of the 
first segment is defined as the abscissa of the point of intersection of P(x) and f(x)—e€ (Fig. 2). 


Inside the interval thus defined, x» <x < X,, the error is | A(x)| < €, and max |A(x)| = € occurs on the 
boundaries of the segment, i.e., for x = x9 and x = x;. 


For the constant sign of f"(x), the parameters of the plecewise-linear function are defined by the equations 


1 (2n) — J" (2n) (Zn — Tm—1) — J (Zm—1) F Ze = 0, (1) 
} (Zm) + 2e — f (Xn) — f’ (&n) (Lm — Tn) = O. (2) 


In these equations, x, and f(x,) + € are the coordinates of 
ia) the points of tangency of the graphs of the linear functions with the 
boundaries of the error band, f*(xn) is the slope of the linear func- 
tion on the unknown segment, X,,_, and f(Xpp.7) are the known 
coordinates of the beginning of the segment of xy, and f(x;p) are 
the unknown coordinates of the end of the segment, 


ffz) 
fz}re t)-€ 


P The subscript m is the ordinal number of the segment (m = 
Y Z, aot = I, Il, Il, . . .) and the subscript n is the ordinal number of the 
tangency points inside the segments (n = 1, 2, 3,...). When the 
doubled sign (+) is used, the upper and lower signs relate, respec- 
tively, to the cases when f"(x) < O and f(x) > 0 











> Fig. 2 


In the general case, Eqs. (1) and (2) are transcendental, and it is necessary to use the method of successive 
approximations to solve them, 


For the solution of Eq. (1), let the rough value x, ~ a by chosen, and 


f(a) — f’ (@) (@ —%m—1) — f (Sm—1) F 2e = Ac, 


where Ac is the discrepancy. 


The correction Ax for decreasing the discrepancy Ac is found from the equation 











on 
Ac 
Az = ; 3 
=F @) Gt, ©” 
Analogously, one can find the relationship between the discrepancy Ad and the correction Ax for the solu- 
tion of Eq, (2). Let 8 © Xp, be an initially chosen value of the independent variable; then the correction is 
Ad 
= — : 4 
7 —T @,) (*) 
of 
— In determining the segments of the approximation, it usually transpires that the end of the last one of them 
ne does not coincide with the end of the interval of reproduction, i.e., x, * xk. One can achieve coincidence of 
u- 


the boundaries of the last segment and of the interval of reproduction by an insignificant change of the error € 


tion), over all the segments of the approximation, 


Let the value at the end of the last segment be x = y,; < x, for € = €4, By choosing a new value, € = €3, 
we find that x = y_>x,. Then, by linear interpolation, we determine the value €, for which x = ys © X,: 


a. = = (&: — 8). (5) 


(x) 


&3 = & + 








In approximating a function with a variable 
sign of the second derivative (Fig. 3), one must 
choose one of the segments in such manner that the 
graph of the linear function is tangent to the upper 
and to the lower boundaries of the error band. The 
coordinates of the points of tangency of the linear 
function to the function to be reproduced are thus; 


fa) 


IZn, / (Zn) + €]; 
[z1, } (i) — s}. 


The boundaries of the approximation segment 
are defined by the abscissas x,,_, and x;,. On such 
a segment, the approximation satisfies the conditions 
of Chebyshev approximations, and the length of the 
segment (for a given size of €) is maximal, The 
error function, A(x) = f (x) — P(x), has the form shown 
in Fig. 4. 

















Fig. 3 


To determine the slope of the linear function and the boundaries of the segment, one must use the following 
system of equations; 


I’ (fn) = 7’ (x), (6) 
(%n — 2) f’ (Zn) — f (Zn) + f (21) — 2e = O. (7) 


From these equations, x, f (X,,), X; and f(x;) are determined. In the general case, Eq. (7) is transcendental, 
and the method of successive approximations must be used for its solution, Let the initial value x, = a be chosen, 
By using (6), we find an intial value of x; » 8. For x, = a and x; = 4, the following discrepancy arises 


Ac = (a —B)/"(«)—7(«) 4 
+ }(B) —2¢. 
The corrections Ax which allows the discrepancy Ac to be decreased is expressed by the equation 


Ac 


At = — ap) 





(8) 


After x, and x; are determined, we can find the coordinates of the left end of the approximation segment 
by using the equation 


J (@n) —f' (&n) (fn —- Fm--1) — J (Xm—1) + 2e = O. (9) 


If, for Xyp-3 = y, the discrepancy Ac arises, then the correction 








Ac 
Az= — 10 
7 @)—T @ 
gives the more accurate value of Xm-1 = y + Ax. 

The coordinate x,,, of the right side of the approximation segment is found from the equation 

} (Xm) — f (1) — 7° (#1) (m — 21) + 2e = 0. (11) 
If, for Xm, = 9, the discrepancy Ac arises, then the correction 

Az = he (12) 


~ FOr) ° 
allows one to determine the more accurate value, x,, = $+ Ax, 
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pP After one determines the approximation segment in 
ae ae nal BOYS aX 4 which f (x) has a point of inflection, one determines the 
J ‘ n approximation segments to the left and to the right by using 
ne 5 fe Eqs, (1) and (2), With this it must be kept in mind that the 
a { boundaries of the extreme approximation segments close to 
the end and the beginning of the given interval [X9; x,] 
Fig. 4 do not coincide with the values of x, and x,. By a small 
change in the magnitude of the error € one can arrange 
for coincidence of the boundary of one of the extreme approximation segments with the corresponding end of the 
interval, To obtain coincidence of the boundary of the other extreme segment with the corresponding end point 
of the given interval, one can decrease slightly the size of the approximating segment which contains the point 
of inflection of the function being approximated, Such a possibility is based on the circumstance that, in the 
segment with the point of inflection, the number of greatest deviations equals four (Fig. 4), while for a Chebyshev 
approximation by linear functions, the minimum number of such deviations equals three, 





é 


The material already given allows us to carry out the piecewise-linear approximation of a function to be 
reproduced by a potentiometer, The following step in the design is the determination of the potentiometer’s 
shunt impedances, for which one must use the equation 


P(2) a & me (13) 


where Rx is the impedance of the lead-in portion of the potentiometer, R is the total impedance of the potentiom- 
eter together with the shunts, ¢ and E are, respectively, the tapped-off and the applied voltages, P(x) is the piece- 
wise-linear function. 


2. Designing a Function-Generating Device with a Finite Number of Parameters 





For the approximation of f(x) by a nonlinear function (x, %)(i=1,..., 1), where r; are the device’s 
parameters, the initial systern of parameters id defined by the condition that, at n different poins, ¢= f. With 
this, there arises the errorAx = f — ¢. To decrease the error of approximation, one must so vary the initial sys- 


tem of parameters that the value of Py an dr; — A(x) differs as little as possible from zero, The number n of 
i=, * P 

parameters is finite; from the set of values of A(x) we choose those which Correspond to a sufficiently dense series 

of values of x, contained in (m—1) intervals, where m > n, The problem of a Chebyshev approximation reduces 

to the search for the best solution of the inconsistent system of linear equations 


n 


Di Fe ars A (z-) (t = i, m: m; m>n),. (14) 


In this paper (Cf., the design example), we use as a procedure the so-called gaugeless equalized descent 
method, developed by E, Ya, Remez [2]. 


We note that, if the approximating function ¢(x, 1) is linear in its parameters, then the initial system of 
parameters can be chosen arbitrarily, With a nonlinear dependence of ¢ on the r,, the procedure mentioned only 
guarantees such a best approximation, which can be obtained by small variations of the initial system of param- 
eters, Generally speaking, it is necessary to try many initial systems of parameters of the approximating function, 
each time finding the best approximation by small shifts of the interpolation nodes, in order to guarantee the at- 
tainment of the best approximation, which is characterized by the impossibility of decreasing |A(x)|n,, by any 
variations of the parameters of the approximation function, 


8. Determination of the Approximation Segments and the Choice of the Inittal 
Magnitudes of the Impedances for a Finite Load Impedance 








With a finite-load impedance, the output voltage loses its piecewise-linear characters, The load impedance 
is ordinarily quite large in comparison with the poteniometer impedance, With this, the output voltage actually 
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Fig. 5 Fig. 6 
differs little from a piecewise~linear one, and the method presented in Section 1 can be used for the determina- 
tion of the approximation segments, 


After the approximation segments are determined, the initial values of impedance of the individual rheo- 
stats are chosen, In order that the possible load oscillations notshow noticeable influence on potentiometer opera- 
tion, it is necessary that the total potentiometer resistance z be several times smaller than the given load resis- 
tance H, 


For a fictitious nonhomogeneous rheostat (Fig. 5), the voltage on the load is 


EHZ 





(%) = WZEE—HUa+D)’ (15) 
and Z is the positive root of the quadratic equation 
HE ’ 
2 ___ -- —_—— — 
Z +(= z)Z—Hz = 0. (16) 


By setting yx) = f(x), and then determining Z(x) from (16), we obtain an exact reproduction of the re- 
quired functional dependence. Since the preparation of rheostats with continuously varying specific resistance 
(nonlinear potentiometers) is technologically quite difficult, one uses instead rheostats, consisting of several 
homogeneous portions, With this, the equation y(x) = f(x) will hold only at individual points, called interpola- 
tion nodes, 


The piecewise~homogeneous potentiometer represented in Fig. 6, with different specific impedances of the 
individual portions, is usually replaced by a shunted potentiometer (Fig. 7), in which the sliding contact moves 
along rheostats R,, .. ., Ry, with identical specific impedances, Impedances ry, Ry, tz, Rg, . . . are so chosen that 
the parallel connection of impedances r, and R; will be equivalent to resistance z,, of impedance rz and R, will 
be equivalent to resistance Zz», etc, 
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Fig. 7 Fig. 8 


4. Computation of the Actual Values of Voltage on the Load Impedance and Set- 
ting up the Correction Equations 


The potentiometer shown in Fig. 6, with the method cited above for choosing the impedances entering into 
it, will give an exact value of the function to be reproduced when the sliding contact is at any of the segment 
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ne 


into 











boundaries, For the other contact positions, the values of the func- 
tion to be reproduced, f(x), and the approximating function ¢(x) 
will not, in general, coincide, and the error A = f — ¢ will arise, 


The approximating function ¢ depends, as will be shown 
below, on the impedances of the individual potentiometer sections, 
The parameters of this function which are subject to variation as 
one attempts to obtain a Chebyshev approximation are assumed to 
be the magnitudes of the shunt impedances and 





n » 

0 

de ’ ar dr;. (17) 
t=} 








To set up system (14) of equations, it is necessary to find an 
expression for the approximating function ¢ which gives the values 
of the actual voltages on load impedance H, 


Fig. 9 


For any position of the sliding contact, the potentiometer can be replaced by the equivalent circuit shown 
in Fig. 8, where the letters H, K,. . . denote the magnitudes of the corresponding impedances, If the contact 
is placed as shown on Fig. 7, then p = r’, and the sense of the other impedances is clear from a comparison of 
Figs, 7 and 8, 


The computation of the quantity ¢ can be carried out by the standard electrical engineering method of 
loop currents, Thus, for example, if one chooses the loopes as shown in Fig. 9, one obtains 


HE(|Roe+(R4+p9+W) KY 


: (18) 





where 


R —H kK 


A=-|R+p+W Ww ons 
Ww H4+w44u U 


Expression (17) can, in general, be obtained by differentiating formula (18), but this way leads to a tremen- 
dously complicated expression. A significantly more convenient formula is obtained as a result of giving a physical 
interpretation to the algebraic relationships given below, 


Let there be given a system of linear equations 


>} a,;¢;=b, (cari, ..., m)s (19) 
j=1 


with a definite solution, If the coefficients of the unknowns are given infinitesimal increments da_,, then the 
unknowns will also receive infinitesimal increments dx;, The new values of the unknowns will be found from the 
system of equations 


n, 
> (@.5 + da-;) (x; + dz;) == be. 


j=l 


If we discard second-order quantities form these expressions, and take into consideration that the x; satisfy 
Eqs. (19), we obtain the following equation for finding the dx; 


>) ajax; = — D} a;da,}. (20) 


j=1 j=1 


To determine the currents in the loops of an electrical circuit, one sets up the equations which express 
Kirchhoff's first and second laws, When the loop impedances change, the equation which expresses Kirchhoff's 
first law does not change. The equations for the second law are transformed to the form 


>) 2,1; = — Dj1,aZ;. (21) 


Equations (21) allow the following assertion to be made, 


With small variations of the loop impedances in an electrical circuit, the variations in the currents will 
be approximately equal to the currents, which would be induced in the circuit by the compensating* electromot- 
ive forces connected in series with the changed impedances. 


If it is assumed that only impedances Zio Tis + - +» M and Zj are subject to correction, then only K, p, U 
will be varied in the circuit of Fig. 8, and we obtain the scheme of compensating emf‘s shown in Fig. 9 and cor- 
responding to increases in K, p and U, 


Now, for the computation of ¢, it will only be necessary to compute the current dI,, through impedance 
H under the influence of the compensating emf's, After this, dg is found from the formula dg = Hdlj;. 


By carrying out the computations, we get 


dp = = { 


R+p+w — | R 
Ww U||R+'p+W 


R K\|\|jR+K —H R —'K —R H+K 
+l o\ |e aw wi%+|rap+w ol lete+w w +p |Z}: 
Formula (22) serves as the basis for setting up correcting Eqs, (14). Computation of formula (22) is greatly 


eased by the fact that A and three of the determinants within the curly brackets were previously computed in the 
process of determining ¢. 


yw |ak + 





APPENDIX 


A Design Example 





Let it be required to reproduce the function f(x) = sin x in the interval 0 <x < 1.4486 with an error of 
about 0,01 for a bounded magnitude of load impedance, 


On the basis given inSection 1, we determine the four segments: 


[0; 0.5323},  [0.5323; 0.8814]; (0.8814; 1.1868], [4.1868; 1.4486]. 


Since (0) = 0, we shall assume, to simplify the potentiom- 
eter, that the sliding contact also moves along impedance Z,. 
Then, for the four segments of the approximating potentiom- 
eter, we shall have the scheme shown in Fig, 10, 


is 

We take H = 50 ohms, the total potentiometer resistance 
(unloaded) as z = 14 ohms and the voltage applied to the poten~ 
tiometer as E = 1,04 volts, As the initial interpolation nodes, 
we take the ends of the segments, By carrying out, under the 
assumption made, the computations described at the end of section 3, we obtain the following values for the 
individual impedances: 














Fig. 10 


Z,. = 7,311, Le 14.1302, re=> 3.3808, rs => 0.9160, 
Z,= 0.507, Ry = 4.7948, R,= 4.1946, Rs = 3.5957, 


expressed in ohms, 





* The compensating emf in a loop equals the product of the initial current in the loop by the absolute magnitude 
of the increment of impedance, With an increase in impedance, the compensating emf is directed counter to the 
initial current, and for a decrease parallel to it. 
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TABLE 1 












































x | f (x) = sin =x @ (x) 100 (f — 9) x { (x) = sin x @ (x) 100 (f — @) 
0.1 0.09983 0.09956 0.03 0.8814 0.77163 0.77172 —0.01 
0.2 0.19867 0.19543 0.32 0.9 0.78333 0.77878 0.46 
0.3 0.29552 | 0.28927 0.62 1.0 0.84147 0.82150 2.00 
0.4 0.38942 0.38258 0.68 1.4 0.89121 0.87333 41.79 
0.5 0.47943 0.47679 | . 0.26 1.1868 0.92717 0.92719 0 
0.5323 0.50752 0.50764 | —0.01 1.2 0.93204 | 0.92792 0.41 
0.6 0.56464 0.55456 1.01 1.3 0.96356 0.94209 2.15 
0.7 0.64422 0.62686 1.74 1.4 0.98545 0.97192 1.35 
0.8 0.71736 | 0.70350 1.39 1.4486 0.99254 0.99265 —0.01 

TABLE 2 
t x any are ae ax4 by 
4 0.4 0.458 2.186 4.529 7.130 —0.3 
2 0.2 0.905 4.324 8.958 14.103 —3.2 
3 0.3 1.357 6.480 13.425 21.135 —6.2 
4 0.4 1.825 8.717 18.059 28.431 —6.8 
5 0.5 2.324 11.099 22.995 36.201 —2.6 
6 0.5323 2.494 41.911 24.676 38.848 0.1 
7 0.6 1.776 13.181 27.308 42.994 —10.1 
8 0.7 0.754 15.223 31.538 49.651 —17.4 
y 0.8 —0.216 17.490 36.235 57.045 —13.9 
10 0.8814 | —0.963 19.633 40.674 64.033 0.1 
11 0.9 —().928 18.127 41.158 64.796 —4.6 
12 1.0 —0.743 10.506 44.114 69.449 —20.0 
43 4.4 —0.562 3.584 47.765 75.197 —17.9 
14 1.1868 | —0.403 | —1.925 51.642 81.301 0 
45 1.2 —0.390 | —1.865 48 . 361 81.418 —4.1 
16 1.3 —0.299 | —1.429 25.882 83.236 —21.5 
17 1.4 —0.212 | —1.041 6.695 86.756 —13.5 
18 1.4486 | —0.169 | —0.809 —1.676 89.167 0.1 




















In using the initial values just found for the impedances to compute ¢, we obtain the results shown in Table1, 
Then, by means of formulas (22), we obtain the correcting equations, given in Table 2, These equations were set 
up under the assumption that impedances ry, rp, ts, and Zj were subject to variation, 


We assume the equations are written in the form 


a,,4r; + a,,4rs + a,,Ars + 4,42; = b_. 
where 


b, =A (z,) = tf (=,) acd. (z,). 


In setting up the correcting equations, we used the following expressions for K, p and U. 


1, The sliding contact on impedance Z,: 


ao Riri Rars Roars 
K=O pmo U= tn + Rtn + Retr te 





2. The sliding contact on impedance R;: 


Rare Roars 
K=Z,, p=n, U= Ritr + Ber, + Z% 
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3. The sliding contact on impedance Ry: 


Riri Roars 
Kat Rp P= T= Beery + 


4, The sliding contact on impedance R;: 


Riry Rare 


K=2.+ Rr, + Retr? Cue U HX, 





In correspondence with these, we obtained the following values of dK, dp and dU on the segments cited, 
On the first segment 


R, \2 R, \2 Rs 2 
dK=0, dp=0, dU= (45) dry + (asx) dra + Regan) dr, + dZ,, 
On the second segment 


R 2 R 2 
dK =0, dp=dr,, dU = (as) drs + Ges) drs + dZ,, 


On the third segment 


Ry 2 Rs 7 
x= (ate! ee ae HS (Retr) ore + oy. 


The Chebyshev solution of the equations of Table 2 will be 


Ar, = 1.2171, Arg =— 0.0882, Ars = —0.0162, AZ; = —0,1194, 


in accordance with which we obtain the following definitive values of the corrected impedances; 


ry = 15,34739hm 2 = 3.29260hm rs = 0.89980hM = Z, = 0.3876 ohm 























TABLE 3 
x 1008 1008* x 1008 1008* 
0.1 0.02 0.06 | 0.8814 1.09 1,13 
0.2 —0.22 | —0.21 | 0.9 0.61 0.66 
0.3 —0.47 | —0.36 | 1.0 —0.94 | —0.91 
0.4 —0.48 | —0.45 } 1.1 —0.72 | —0.71 
0-5 —0.00 0.04 | 1.1868 1.09 1.10 
0.5323 | —0.29 0.34 | 1.2 0.67 0.68 
0.6 —0.58 | —0.53 | 1.3 —1.09 | —1.08 
0.7 —1.09 | —0.88 | 1.4 —0.28 | —0.28 
0-8 —0.51 | —0.48 | 1.4486 1.09 1,09 


With these impedance values, we might, on the basis of the correcting equations, expect the error 4 given 
in Table 3, However, it should be kept in mind that the correction equations give approximate results (since 
finite, and not infinitesimal, variations of the impedances are made), For this reason, Table 3 also includes the 
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actual values of the errors 5* , computed by means of formula (18) on the basis of the definitive values of the 
impedances, A comparison of the values of 5 and 6* show that they accord satisfactorily with one another. 
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TEMPERATURE STABILITY AND ADMISSIBLE VALUES OF 
POWER AND VOLTAGE ON TRANSISTOR COLLECTORS 


B. S. Sotskov 
(Moscow) 


We use the notation: 0; is the n—p junction temperature, P; is the power dissipated at the n—p junction, 


equal to Py = VI, Rey is the thermal impedance with heat transfer from the n—p junction, The value of Rep 
may be presented (Fig. 1) as 


Rix = Ro + Rim, 
where Rr» is the thermal impedance for the internal portion (from the n—p junction to the body of the transistor) 


and R,,,, is the thermal impedance for transfer of heat from the transistor body to the ambient medium, The 
increase of temperature 6; for the n—p junction layer is defined as 





d0; = RirdP. (1) 
On the other hand, the power dissipated in the n—p junction layer is defined 
4% by the relationship 
$ 
3c dP, = Vidl.. (**) 
P 
t ’ Since I, = Icg + Igy, where Ico is the portion of colector current which does 
h,, f; not depend on the cotnrolling voltage and Icy is the portion of collector current 
: - which does depend on the controlling voltage but does not change with tempera- 
— Sy) oa ture, then 
6 . 
Fig. 1 dP; = Ved (Leg + Lev) = Vedleg. 


The dependence of current Icg on ternperature is expressed as 
Teg = Tenge? 4 = Teq910*°5, 
where the value of a equals: for germanium, « * 0,08 to 0.85, or a's 1/ 30, and for silicon a 1/40, whence 


it follows that 


dI co = Ieqe**i ad. 
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Temperature stability for transistors is defined by the condition 
that the increment of power dP, dissipated in the n—p junction layer 
not exceed the increment of power dP transferred through thermal 
impedance Rey l.€., 


dP, <P. 





By substituting the values for dP, and dP, we find that 











Risa legge ™ i < 1. 
Fig. 2. 1) With radiator for nor- As 6, increases, the left member will be increased, For limit- 
mal air pressure; 2) without radi- ing admissible value, j;, we have that the limiting admissible value 
ator for normal air pressure; 3) of collector voltage, Vc = Vo. cri, equals 
without radiator for evacuated air; 
nxtany / ny = Ry, where n, is the V < mE 34 
scale of the axis of abscissas in crit ““ “align Rex ° 
grams /cm and ny is the scale of 
the axis of ordinates, in watts/ cm. With this, Rj, must be chosen with the largest value in accord 


with operating conditions, If there is the possibility of the transistor 

operating in an evacuated atmosphere, for which Rim increases due 
to a decrease in the heat transfer coefficient, then it is necessary to choose an Rj, which corresponds with the 
worst heat transfer conditions, The value of 6); is given on the rating sheet for the given transistor type. 


When the load impedance R, # 0, the value of V,, is 
V. = er T.R,, 


where Vc, is the voltage of the collector circuit's power supply. With this, the power dissipated in the collector's 
n—p junction layer is defined as 


Pt = (Veg —1 Pi) Ie, 


and, consequently, 





dP, = (Veg — 21 ¢R,) dle = (Ee — 21-R,) deg 


so it will follow, form the condition that dP, < dP, that 


Riza (Veo — 2 eR) Leae** i <1. (m) 


It is clear, from a comparison with (I), that the presence of a load impedance in the collector circuit makes 
it easier to attain the conditions of thermal stability of the transistor, 


By knowing the limiting admissible n—p junction layer temperature, 6 ;};, and the normal temperature of 
the ambient medium, 69, we have the following relationship for the admissible power dissipated in the transistor 


PR = Vl Rix = Oi — 4,. 


In case the temperature of the ambient medium is 6, = 6, + 46x, then the admissible transistor power is 
defined by 


Pu = i (9515 — (+ A0,)}. 
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Figure 2 shows P,, as a function of 6, for various values of 


a a= Rip. Mounting a transistor on a radiator is analogous to connect- 
YG Vo 


ing a secodn thermal impedance 'R,, in parallel with R;,, and, con- 
sequently increasing the total thermal resistance 


R 


Ry = Rig + Rim (x-rx_) 





and increasing the admissible power P,, for the same values of 














Let the characteristic |, = f(Vc, I) be given for a transistor. 
Then, by knowing o, Ic 99, 9j)j. Ox = 65+ 40, and Rese we find 
the admissible values for: 








Fig. 3 a) The power dissipated in the transistor 
Py. = — (8; — 9]; 
Ris 


b) The critical value of collector voltage 


on ai 
Vecrit SSR” 


For this it is necessary to choose the greatest values of I, 99 and R,; which arise in the given transistor, On 
the characteristic i. ® f(Ve. l,) (Fig. 3), we draw the functions Ic = Pry / Vo and Vo crit = Const (dotted lines), 
The admissible I,, Ve region for extended transistor operation is bounded ~ the coordinate axes and the dotted 
curves just cited, 
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A DIGITAL SERVO SYSTEM FOR THE AUTOMATIC PROGRAMMED 
CONTROL OF A REVERSIBLE ROLLING MILL SCREWDOWN 
MECHANISM 


V. M. Ozernoi and G, S, Rozenberg 
(Moscow) 


A short description is given of a digital servo system for the automatic 
control of the screwdown device of a sheet steel hot-rolling mill, 


The circuits of some of the blocks, as well the formulas of the basic 
logical cells, are given, 


The development of the complex automation of productive processes creates a host of problems related to 
the programmed control of mechanisms, The programmed control problem arises, in particular, in the automat- 
ing of the screwdown mechanisms of reversible rolling mills, 


A system of automatic programmed control of a revers- 
ible rolling mill's screwdown mechanism must provide an ac- 
curate placement of the mill's upper roller for each pass of th 
metal in correspondence with the given program of width re- 
duction, There may exist several programs for each rolling 
mill, depending on the assortment of ingots or sheets to be 
rolled, 





For the automatic displacement of the upper roller, elec- 
tric servo systems are used, whereby these systems may be 
either continuous or discrete (digital), 





In continuous (analog) servo systems, the actual position 
Fig. 1 of the upper roller is transformed to a continuous electric quan- 
tity, for example, to the amplitude of a dc or ac voltage, to 
phase, etc, The transducers which measure the upper roller's position and transform it to a continuous electrical 
quantity are ordinarily selsyns, phase inverters, rheostats and inductive potentiometers, Preparation of these 
transducers with a high degree of accuracy of transformation is very complicated, 


For certain types of mills, wehere there are particularly high requirements on the accuracy of upper-roller 
placement, the use of transducers of the continuous type is only possible when a system of rough and accurate 
(vernier) computation is used, which significantly complicates the circuitry, The set of programs for screwdown 
mechanism operation, and also the devices for remembering these programs in the analog systems are very cumber- 
some and inconvenient, This is particularly essential for mills which operate with a large number of programs, 


In digital servo systems, the actual position of the upper roller is transformed to a number, usually expressed 
in the binary number system, i.e., in a binary code, The program for screwdown mechanism opetation ts also 
given in the form of a binary code, The use of digital servo systems allows one to obtain a high accuracy of 








screwdown mechanism operation, and to implement the control of screwdown mechanism operation with a min- 
imum of apparatus, as well as a rapid and easy switch-over of programs by means of punched cards, 


We give below the description of a servo system, developed by the Central Automation Laboratory of the 
Construction Ministry of the RSFSR, for the automatic control of the screwdown device of a sheet steel hot-rol- 
ling mill, The basic technical requirements to + 0.1 mm for an operating travel of 400 mm, and the capability 
of fast and easy switch-over programs for screwdown mechanism operation with a large volume of programs, These 
requirements conditioned the selection of a digital servo system, 


The block schematic of the servo system is shown on Fig. 1, where the following notation was introduced; 
1) is the giving of the program by means of a punched card; 2) is the program-furnishing device; 3) is the 
arithmetic device; 4) transforms the Gray code to an arithmetic code; 5) transforms the difference code to a 
voltage; 6) is an electronic amplifier; 17) is an electromagnetic amplifier; 8) is the generator; 9 is a motor; 
10) is a reducing apparatus; 11) is the screwdown mechanism; 12) is a selsyn-transducer; 13) is a selsyn-dif- 
ferential; 14) is the null correction system; 15) is a selsyn-receiver; 16) is a reducing apparatus; 17) encodes 
the position of the screwdown mechanism; 

















Fig. 2 


The program of screwdown mechanism operation, i.,e., the given position of the mill's upper roller at each 
pass, is written in the binary code in the form of round punches on a metal-disk punched card, There is just one 
program written on each disk, To switch programs, it is necessary to replace the disk in the program-furnishing 
device (Fig, 2), The disks are divided into sectors by radii, The number of sectors equals the maximum number 
of passes for rolling mill operation, The holes (punches) placed on one radius determine the given position of the 
upper roller in the corresponding pass. 


For reading the program, germanium photodiodes are used, incandescent lamps serving as the light sources, 
The photodiodes and light sources are distributed on both sides of the disk, and are maintained immovable, The 
punched disk, by means of a special stepping mechanism, is rotated one step at each pass, so that a definite radius 
with punches is positioned between the photodiodes and the light sources, 


When a photodiode is illuminated, the magnitude of its back impedance is decreased by several ten-fold. 
The signals form the photodiodes are amplified and are presented in the form of high and low potentials (with 
respect to ground) in parallel for all bits of the code, An illuminated photodiode provides a low potential, i.e., 
0, to the corresponding bit, while a darkened photodiode provides a high potential, i.e,, 1. The number of 
photodiodes in the program-furnishing device equals the number of bits in the code, 


The measurement of the actual position of the rolling mill's upper roller is done by means of a coding 
transformer, The number of bits in the transformer is determined by the formula 


n = logs —, (1) 


where A is the working travel of the upper roller and q is the measurement accuracy, 
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With the given working travel of the upper roller, it is necessary to use a 12-place (12-bit) transformer to 
obtain the desired accuracy, corresponding to an accuracy of upper roller measurement of + 0.025% of the maxi- 
mum value of the measured quantity, 


The use of the binary arithmetic code in the coding transformer could lead to ambiguous transformations 
due to indeterminacies arising at the boundaries between two discrete neighboring values, Special measures 
were taken to eliminate ambiguities, One of them is the use of a cyclical progressive code (a Gray code), The 
distinguishing feature of a Gray code is that only one bit changes when a transition is made from one number to 
an adjacent one [1] This makes it possible to decrease the errors due to ambiguity in a transformation, involving 
a low-order bit unit change, 





Fig. 3 


Figure 3, a shows a disk with a Gray code, containing the six low-order code bits, and Figure 3,b shows 
a disk containing the higher-order six bits, The disks are prepared from organic glass; the code patterns are placed 
on the disks by photochemical means, The code disks are connected togher by means of a discrete reducer with 
a 1:32 transfer ratio, and are connected to the screwdown mechanism by means of a selsyn transmitter. 


The selsyn-transducer is mounted on the rolling mill housing and, via a reducing apparatus, is connected 
directly to the motor shaft of the screwdown device (Fig. 1). The selsyn-receiver is located with the code disks, 
and is joined to the first disk via a step-down reducer, A selsyn-differential is connected between the selsyn- 
transducer and the selsyn-receivei for the null correction of the servo system, 


The use of intermediate selsyn transmission allows the coding transformer to be placed, not on the rolling- 
mill housing, where it would be subject to vibration, high temperatures and emulsion, but in more favorable con- 
ditions, The transmission factor of the step-down reducer between the selsyn-receiver and the coding transformer 
was chosen sufficiently high so that the errors of selsyn transmission would not lower the accuracy of operation of 
the digital servo system, 


The reading of inform~tion on the actual position of the upper roller from the code disks is implemented by 
germanium photodiodes and the light sources, The signals from the photodiodes are amplified and, for each dis- 
crete position of the upper roller, there is presented, in parallel, a combination of high and low electrical poten- 
tials in 12 orders (bits), Just as with the furnishing of the program, a low potential corresponding to a 0, and a 
high potential to a 1, 


The actual position of the upper roller, thus coded, cannot be directly compared with the given code, since 
the Gray code is not arithmetical, Therefore, it is necessary to transform it first to the binary arithmetic code. 


The transformer of the Gray code to the arithmetic code, and the arithmetic device itself, consist of the 
simplest circuits, which execute the logical operations of addition, multiplication and negation, The controlling 
signals at the input and output of the circuits are combinations of high and low electrical potentials, 


The use of potential elements allows us to obtain simple and reliably operating circuits distinguished by 
high stability and noise stability, The synthesis of the basic logical circuits was done on the basis of the apparatus 
of the algebra of logic (Boolean algebra), The basic criterion in choosing the circuits was to minimize the number 
of vacuum tubes, which simultaneously implement the functions of amplification and inversion (logical negation) 











of their input signals, In the synthesis a matching condition was adhered to, whereby it was inadmissible to con- 
nect an “or” circuit in series with an “and” circuit, since such a connection cannot be carried out with diode 
elements without the use of intermediate amplification. 


We consider now the circuit for transforming the Gray code into the binary arithmetic code, We denote the 
digits of the successive orders (bits) of the Gray code, in the order of increasing order weight (increasing signif- 
icance) by a4, a9, a3, . . ., 4, and the digits of the corresponding order of the binary arithmetic code by yy, y», 


Yar ++ ++ Yn 


The transformation of the Gray code to the arithmetic code is done serially, starting with the most signif- 
icant bit, in accordance with the following logical rules: 


a) for the most significant bit 


Yn = 4ni 


b) for the remaining bits 


f0, ecif Yrtr=@e 
Ye = , 
\1, ec if yess eae. 


The transformer of Gray code to arithmetic code consists of 12 identical cells, The logical formula of each 
cell is analogous to the formula for a one-bit half-adder [2}: 


Yr = Yay + Deas (2) 


Use in the transformer of cells constructed in accordance with formula (2) is not possible, since the connec- 
tion of such a cell with another transformer cell or with the arithmetic device would give rise to the inadmissible 
connection of an “or" circuit to an “and” circuit, 


To obtain the optimal circuit, we transform formula (2) to the form 











Yu = Ax + Year t GeYrir (3) 
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Fig. 4 


Figure 4 shows the block schematic (a) and the simplified functional schematic (b) of one transformer cell 
corresponding to formula (3), as well as a table of the logical operation executed by the cell. 


A comparison of the given (desired) and actual positions of the upper roller, coded in binary, is implemented 
in the arithmetic device, The codes to be compared are applied to the arithmetic device's input; at its output, 
the difference of these codes is obtained. To obtain this difference, subtraction is replaced by addition of the 
complemented code (of the subtrahend), the complemented code being obtained from the uncomplemented (direct) 
code by replacing all ones by zeroes and all zeroes by ones, 


To avoid additional transformations, the given program is implemented in the complemented code, 
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Replacement of subtraction by addition of the complemented code makes it possible to use a parallel 12- 
bit adder in the arithmetic device, Each one-bit cell has three inputs and two outputs, The three inputs are the 
complemented code from the corresponding bit of the program-furnishing devite, the direct code of the actual 
position from the transformer of Gray code to arithmetic code, and the carry signal from the next lower-signif- 
icant bit, The two outputs are the sum bit, and the carry bit to the next higher significant bit, The carry signal 
from the output of the most significant (twelfth) bit cell is applied to the input of the least significant bit cell, 


The logical formulas for a one-bit adder has the following form [2): 


S, = [(z, + Y,) (x, + Y_)+P,_1] (r,¥, + Vp + P,_.): (4) 
Py = [(2, + Yy) (2, + y,)I Py_y + yy. (5) 


Here, S) is the sum bit at the cell's output, x, is the k’th bit of the complemented given (program) input 
y, is the kth bit of the direct code of the actual position, P;,., is the carry from the next lower-order cell and 
Py is the carry signal to the next higher-order cell, 


Formulas (4) and (5) are not in accordance with the matching condition, We transform them to the form: 














Se = (Fq YQ) Pay HEV Py F TY HF (Fe A Vg) + Ppp (6) 
(7) 





Pe = Bye + Yq + My + Py_y + Tee: 


It follows from expression (6) that the inverted (complemented) carry signal from the next lower-order bit 
must be applied to the input of the current bit, The carry signal to the next higher-order bit must also be in- 
verted 






































































































































Py = 2% + yr 4 Teyx + Pry + Lrye- (8) 
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Fig, 5 


Figure 5 shows the block schematic (a) and the simplified functional schematic (b) of a one-bit adder, con- 
structed in accordance with formulas (6) and (8). Figure 5 also gives the table of the logical operations imple- 
mented by the one-bit adder cell, 





-With subtraction by the method of complement addition, the difference at the adder’s output may be ob- 
tained either in the direct or in the complemented code, depending on its sign. 


If the difference is positive, then it is given in the direct code, and there is a carry signal from the output 
of the most significant (twelfth) bit of the arithmetic device, If the difference is negative, it is given in the 
complemented code, and there is no carry signal, Thus, the presence or absence of a carry from the most sig- 
nificant bit indicates both the sign of the difference and the code in which it is expressed. If the difference ts 
obtained in the complemented code, it must be transformed to the direct code. 
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Fig. 7 


The block and functional schematics of a one-bit cell for choosing either the direct or the complemented 
code are shown in Fig. 6,a and b, The table of logical operation is also shown there, 


The logical formula of the cell has the form: 
Fy = SpP y+ Se t+ Pr. (9) 


Here, F,, is the output bit of the code choice cell, Sy is the sum bit at the output of the adder's k’th bit 
and Py, is the carry bit from the highest-order (twelfth) bit of the adder, 
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The block schematic of the logical circuit of the digital servo system is given in Fig, 7, where a shows 
the transformer of Gray code to arithmetic code, b shows the adder and c is the circuit for choosing between the 
direct and the complemented code, ' 


For the control of the electric drive system, the difference code obtained is transformed to a sinusoidal 
voltage whose amplitude depends on the discrepany (i,e., the error), 


Stability, and the required transient response quality in the system, are attained by introducing flexible 
feedback paths in continuous (analog) form, For this, it suffices that the controlling voltage depend linearly only 
on those errors which are close to zero, Because of this, only the six lowest-order bits are transformed to a pro- 
portional voltage, this corresponding to an error somewhat in excess of the braking travel of the screwdown device 
motor when oen brakes from the maximum velocity to zero in the case when the maximum braking torque is ap- 
plied. 
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The transformation is effected by the method of summing currents, The circuit for transforming difference 
code to voltage is shown in Fig. 8, It contains seven pentodes, The resistors in the vacuum tubes" cathode cir- 
cuits provide the flexible feedback and, consequently, a high stability of the ac components through the pentodes, 


Tubes L1-L6 being cut off corresponds to a 0 in the corresponding bit of the difference code, the conduct- 
ing of any of them corresponding to a 1 in the corresponding bit, Tube L7 will conduct if there is a 1 in at least 
one of the six higher-order bits, In addition to the signals of the difference code, sinusoidal voltages of identical 
amplitude are applied to the control grids of the tubes, If a tube is conducting, this voltage determines an ac 
component in the anode current, In the tubes’ anode circuits are connected resistors whose magnitudes are mul- 
tiples of the numbers 1, 2, 4, 8, 16 and 32, The ac component of the resulting voltage, taken off the resistors 
by means of the transformer, equals 


u=i(Ro+Ry+...+Rr), (10) 


where { is the ac component of the tube's anode current and Rg, Rp . . ., Ry, are the total impedances in the 
anode circuits of the tubes which are conducting, 


The phase of the error voltage U, thus obtained is chosen as a function of the sign of the error, determined 
by the presence or absence of a carry bit from the adder's twelfth bit, . 


A constant potential, equal to half the voltage corresponding to a 1 at the adder output, is applied to the 
middle of the transformer's secondary winding, The carry signal from the twelfth bit equals 0 or 1, correspond- 
ing to a high or a low potential, and permits a voltage of the necessary phase to be obtained in the load on the 
transformers output, 





The error voltage thus obtained is applied to a phase-sensitive electronic amplifier which forms the voltage 
which controls the electric drive system. 












The operation of all the system blocks was verified on mock-ups, They are simple to adjust, very non- 
critical to variations in circuit parameters, stable under external noise and operate stably when the supply voltage 


varies by + 15%, 
The digital servo system herein described can be used, not only for the automatic control of rolling mill 


screwdown mechanisms, but also for other program-controlled mechanisms, 
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A COMPENSATION INSTRUMENT FOR TRANSFORMING 
POWER TO DIRECT CURRENT 


M. Duma and M, Syrbu 


(Bucharest) 


A compensation instrument, developed by the authors at the Research 
Institute of Electrical Engineering of the Ministry of Heavy Industry of the RPR, 
for transforming power to direct current,is described, The instrument's con- 
struction is based on a standard three-phase, three-stage meter, 


In power telemetry technology, there is frequently the question of tranforming power to direct current, 
both as a parameter which is convenient for short-range transmission, and also for the summation, and transmis- 
sion of the summation indicators, of electric substation power, 


In the transformers of power to direct current, which have hitherto been developed, the following properties 
are inherent, 





T J ¢ a) The existing static transformers [1-4] can provide ac- 
AT MS TA —+{/) curacies of from 1 to 1,5%, but possess a number of disadvantages, 
For example, the thermoelectric transformers [2; 3], provide low 
v, Ic dc power levels at their outputs, so that special amplifiers are 
| CS required if they are to be used in short-range telemetry devices, 











The magnetic transformers [1, 4] require that complex operations 
of sorting, calculating and regulation be implemented, which is 


his not readily applicable to mass production, 


b) The existing transformers which contain moving parts, operating in closed control loops, for example, 
in accordance with the principle of equalizing the active torque induced by a three-phase wattmeter with the 
counteracting torque produced by a magnetoelectric milliammeter, can provide the necessary output power for 
long-range transmission of indications, but require specially developed accurate mechanical construction, which 
is undesirable. The accuracy of these devices depends on the accuracy with which the mechanical construction 
is implerr ented, 


In 1957, the authors, working in the Research Institute of Electrical Engineering of the Ministry of Heavy 
Industry of the RPR, developed a compensation transformer of power to direct current, with class 1 accuracy and 
allowing telemetry of indications to be made over distances of tens of kilometers over wire communications 
channels, The transformer is easy to reproduce and is of simple construction, since it is constructed of standard 
elements of large-scale mass production, 


The principle of operation of the developed transformer of power to direct current (Fig, 1) consists in the 
equalization (balancing) of the active measured torque T,,,, generated in the two measuring systems MS of an 
induction meter, by a counteracting torque Tc, generated by means of transformer amplifier TA in the counter- 
acting system CS of the same meter, The counteracting torque T,, follows the changes in Typ and tends to equal- 
ize (balance) them, As is generally known, the torque developed by the moving system of an induction instrument 














with moving poles is proportional to the power measured: 


T m= Klin where a = Ut £08 ?m. 


*, = K. P. , where P., = U. I £08 9. 


Here the subscripts "m" and “c* denote the quantity relative, respectively, to the measured and counter- 
acting systems, 


For a balanced system, the counteracting torque equals the active measured torque, i.e., T,,, = Te. If the 
voltage U, of the counteracting system is maintained constant, and the angle ¢, equals zero, then the measured 
power is proportional to the actual value of the counteracting current: 


K, Ua 


Pit Km 


a 


After rectification , the average current value I, at the output of the transformer amplifier will be pro- 
portional to the effective value of the ac current I, of the current winding of the counteracting system, In the 
case of two-half-period rectification, if we consider K, as the transformation coefficient of current transformer T, 
we obtain 


Ta= aK, a 
Then, 
uK,K,U 
ns nt Sh. 
m 4K ny 


As is seen from the last expression, the magnitude of the direct current at the rectifier's output is proportional 
to the value of the measured power. 
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Fig. 2 


To derive these formulas, we made the following assumptions, 


1, The expression U,, cos 9, is constant even with oscillations of the line voltage over the entire range of 
transformer operation: this can be accomplished by stabilization of the value of Uc cos ¢, or by stabilizing only 
the effective value of voltage Uc, of the counteracting system for an angle ¢c between U, and Ic approximately 











tonal 











equal to zero, such that cos ¢, is alwyas close to unity, Then, small variations in the angular shift between the 
two vectors during system operation do not have significant effect, 


2. Over the whole range of transformer operation, the coefficient of proportionality, Ky, between the average 
value of output current I, and the effective value of counteracting current I, remains invariant, This means that, 
over the entire range of transformer operation, the transformation coefficient K; of current transformer T; (Fig. 2) 
remains invariant, 


3, There are no distortions of the form of the counteracting current curve which would disturb the proprotion- 
ality between the average value of the current and the effective value of its first harmonic component, The 
counteracting system's voltage winding is supplied by a virtually undistorted line voltage,so that only the first 
harmonic induces a counteracting torque. 


4, Over the entire transformation range, the ratio K,/K,, remains constant, In actuality, the error curve 
of an induction meter for variations in load is analogous for each system, Consequently, the changes of coef- 
ficients K, and K,, with the load follow the same law, so that the ratio K-/K,,, can be considered constant, 


Figure 2 shows the simplified schematic of the compensation transformer of power to direct current, As the 
comparison element, we used a standard three-phase, three-stage electric energy meter in which the necessary 
changes were made; the first of these was the removal of the metering mechanism, 


Measurement of three-phase power is implemented by the method of two wattmeters in measuring systems 
MS1 and MS2, On the meter's moving system there acts a rotational torque proportional to the measured power, 
Windings W, and We, together with capacitors C, and Cz, comprise the tuned circuits of a high-frequency gener- 
ator, A notch was made in the lower disk, providing an increase in the coupling factor between the two windings 
when the disks are rotated to the side corresponding to the effective active torque, induced by systems MS1 and 
MS2, by which the amplitude of the auto-generator's high-frequency oscillations is correspondingly varied, 


When the measured power equals zero, there are no high-frequency oscillations, and output tube 6N7S is cut 
off by the voltage from transformer T, which is in phase opposition to the anode voltage, When an active torque 
from the measured power appears, thus generating high-frequency oscillations, tube GN7S conducts to the extent 
required to induce, in the counteracting system CS, the ac current I, necessary to equalize the torques, It is 
clear from the simplified schematic that the dc current of the counteracting system is produced in current trans- 
former T; by the addition in the secondary winding x of the two half-wave currents in the primary windings, The 
primary windings are the anode loads on the two halves of tube 6N7S, operating in a rectifier-controlled two-cycle 
circuit, The rectified current thus obtained is filtered and directed to the communications channel for power 
telemetry, or to the control winding of summing device SD if summing is carried out at the sending points, 


It is clear from the foregoing exposition that the transformer developed is a balanced-current instrument, 
In such devices, one eliminates almost completely the errors stemming from variations in circuit parameters, 
line parameters, voltage supplies, aging of the electronic tubes of the transformers electronic amplifier, etc, 
Since, the instrument is astatic, the error rapidly tends to zero if any such external excitation appears, 


The physical construction of the transformer is comprised of two blocks (Fig, 3); the compensator block, 
with dimensions of a three-phase meter, and the electronic block, 


Testing of the compensation power transformer showed that basic instrument error did not exceed 41% 
this error being the deviation of the characteristic from linearity with changes of the input quantity from zero to 
120% of its nominal value, when the power coefficient was cos g, = 1, the frequency was f = 50 cycles and the 
temperature was t = 20°C, The graph on Fig, 4 is the average of 40 measurings, carried out during a period of 
200 days. The axis of abcissas give the value of the measured power (in percentages of its nominal value) and 
the axis of ordinates given the basic transformer error (in%), The maximum value of the error was 0.8%, 


The additional errors when the temperature varies by + 20°C, the supply voltage varies by + 15%, the fre- 
quency varies by + 2%, the communications line impedance varies from zero to 1500 ohms, and the distortion 
coefficient of the supply voltage varies by 5% do not exceed 1% each, 


With a two-disk meter construction, the additional error from the variations of cosy amounts to 1.5% with 
a three-disk construction it is less than 1%, 






























The time to establish the output current 
= when the input power changes instantaneously 
as a4 a5 06 Q7 08 09 W fie from zero to its nominal value, does not exceed 
0 a Sertrrt peer 2 seconds; the oscillatory process (the transient 
as icant...» response) is rapidly damped. Such a rapidly- 
damped transient response is of great value in the 


case of short-range telewattmeters, since the 
frictional error is such an instrument is reduced, 





Fig. 4, 


The transformer requirements correspond to the standards for meters, 


For different variants of the physical construction, the dc current at the transformer's output, for Uc = 
= 100 volts, l, = 5 amps and cos ¢- = 1, was taken to be from 5 to 20 milliamps. 


The transformer continued to operate satisfactorily under difficult operating conditions: for a change in 
the communications line impedance from zero to 4000 ohms, for a 40% increment of the input power from its 
nominal yalue, with a distortion coefficient up to 10%, the additional errors were of the order of 2%, 


In the case of single-phase power measurement, each transformer can be used for summing two indications, 


The necessity of stabilizing the counteracting voltage might be considered a disadvantage of the transformer, 
Since each stabilizer services four transformers, the cost of the devices is increased insignificantly, The trans- 
former can also work off a nonstabilized source, with the condition that the controlling or summing devices are 
supplied by the same source as the transformer, The errors attributable to these devices when the voltage changes 
must be equal to the transformer errors, but opposite in sign, 


The basic disadvantages of the described electronic-induction compensation transformer of power to direct 
current are the following. 


1, Simple construction, based on standard three-phase meters, The complex artifacts of exact mechanics 
are not necessary, which makes production of such transformers much more economical than the production of 
other existing instruments with moving parts, operating with closed, or with open, control circuits, 


2, Favorable conditions for mass production, since the transformer contains only easily adjustable and 
reproducible elements, which do not require complicated individual tuning, sorting or measuring, or individual 


design, 


3. The basic error (i,e., deviation from linearity) is less than 1% for the transformer of power to direct 
current, 


Such accuracy can be obtained even with the use of class 2,5 meters, because [7} 
a) removal of the metering mechanism significantly decreases the frictional errors; 


b) a mutual compensation of the nonlinearity curves of the measuring and counteracting systems is possible, 
since they are analogous; 


c) when a meter is used in a transformer of power to direct current, the error is referred, not to the given 
measured quantity, as is the case for electric energy meters, but to the nominal value of the power. 


4, Large output power, allowing the transformer to be used as a short-range telewattmeter without an 
intermediate amplifier, 


5. There is no connection between the moving and the immovable systems of the instrument, as there are 
in all other existing transformers with mobable parts, due to flexible filaments, springs, brushes, etc. 


6. Reliability in use, provided by the absence of easily maladjusted elements and by the entire physical 
configuration of the transformer. 


7. The possibility of signalling or controlling protective elements if the supply voltage fails, the counter- 
acting torque disappears, the oscilletions are disrupted or there is an external defect, for example, if the com- 
munication line ts cut, 
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CRITICISM 


REVIEW OF THE BOOK "INTRODUCTION TO AUTOMATIC 
CONTROL THEORY" BY A. YA, LERNER 


(Mashgiz, 1958) 


The book contains a short course in the fundamentals of automatic control theory, It is intended basically 
for the broad circle of engineers who work in the domains of automation of production, the development of 
automata circuitry, the installation and adjustment of automatic control systems, 


From the large number of questions contained in automatic control theory, the most basic, a knowledge 
of which is necessary for the types of readers cited above, have been included in the book, The author divides 
the book into three parts, unequal in size. 


The first part of the book acquaints the reader with the nomenclature of automatic control systems, their 
structures, with the basic technical ideas connected with the use of negative feedback, with the terminology and 
the classification of automatic controllers, The role of automatic devices in the control of productive processes 
is explained, and the basic information, necessary for the static and dynamic analysis of an automatic control 
system, is provided, Unfortunately, this portion does not contain descriptions of the standard apparatus used in 
automating production. The author seemingly assumed that descriptions of such apparatus must be contained in 
courses on the technical means of automation, but such courses do not exist, and for the reader referred to above, 
it would certainly be useful to include in this section basic information about the standard apparatus, their char 
acteristics and capabilities for use in various circuits, 


The second part is the basic portion of the book, and is devoted to the statics of automatic control systems, 
to the dynamic consideration of control systems, chiefly in linear approximations, and to certain data on the 
properties of nonlinear automatic control systems, This book usefully distinguishes itself from other textbooks 
containing approximately the same material by its detailed consideration of the questions related to the static 
relationships of the system “controlled object — controller”, 


The dynamic investigation begins with the classification of linear links and with general information about 
the algebra of linear circuits and about the different types of frequency characteristics (including logarithmic 
characteristics), The chapter devoted to stability contains an exposition of the basic stability criteria (Hurwitz, 
Mikhailov, and phase-amplitude criteria) and a brief discussion on the construction of the regions of stability by 
the method of D-decomposition, Since the Mikhailov and frequency criteria of stability are presented earlier 
than the D-decomposition criterion, the proofs of these criteria is carried out independently of the D-decomposi- 
tion method, being based, essentially, on principles which are not explicitly introduced, 


The chapter devoted to methods of analyzing transient responses is clearly inadequate, both in volume and 
in the amount of material contained in it. In particular, a considerably more detailed consideration was deserved 
by the integral estimates and by the estimates related to the distribution of zeroes and poles (some space was al- 
lotted to these questions in the book's third chapter), 


With respect to the estimates of the form of a transient response from the real characteristic, the author's 
presentation corresponds approximately to that adopted in contemporary courses, although it has recently, par- 
ticularly in the works of Blokh, been shown that one can obtain more complete informtion about the transient 
response, as well as more simply, by the use of the imaginary characteristic, 


The third part, the most interesting in the book, is entitled by the author, “Determination of the Sructure 
and Calculation of the Parameters of an Automatic Control System", Essentially, it considers the problems of 
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synthesizing linear-control systems, plus certain problems related to the creation of optimal nonlinear systems, 
The first chapter of this part (Chapter 10 of the book) is devoted to the questions of greatest importance to 
designers, and almost completely ignored in textbooks on automatic control theory — the selection of automatic 
controllers, and their adjustment to accord with the data on the controlled object and the technical data on the 
control process, The following chapter considers the more general questions, which are also treated in other 
monographs and course, of synthesizing stabilizing and correcting loops for linear systems, However, in contra- 
distinction to other courses, where only frequency methods are used for this purpose, the author also adduces the 
root hodograph method, Many works are devoted today to this method, principally in the United States, In the 
Soviet Union, unfortunately, only the work being reviewed treats of this method, and the author's attempt to 
familiarize the afore-mentioned group of reads with the basic ideas of this method should certainly be applauded, 
The last chapter is also a novel feature of the book, since the questions treated in it have, until now, been con- 
sidered only in scientific works, It is devoted to the theory of systems which have optimal speeds of response, 

In relation to this, the author introduces the basic concept of a process which is optimal as regards its speed of 
response; the reader is familiarized with the most fundamental ideas of implementing optimal processes by means 
of nonlinear connections, Two basic remarks may be made anent this section, One deals with the chapter devoted 
to the selection of controllers, The fact that only 12 pages of the book are given to this chapter speaks for itself, 
and shows that the space which this chapter occupies in the total volume of the book does not correspond to the 
tremendous importance of the questions discussed in the chapter, In designing a scheme for automating produc- 
tion, the engineer-designer has far more frequently to solve the problem of how to select the controller, and the 
engineer-adjuster, how to tune it, then either of them has to carry out a dynamic, or even static, calculation of 
the system, Certainly, it is difficult to make recommendations on choice of controller, since not enough factual 
material, unfortunately, has been accumulated for this by automatic control theory, but the book's author is 

also well known for his work in practical domains, and he could be required to have devoted greater attention 

to this question, 


The second remark has to do with optimal system theory. Essentially, the author limited himself to present- 
ing his own investigations in optimal system theory, They are interesting but, since, they were published, a 
number of first-class works have appeared, both in Soviet Union and in the United States, in the direction of the 
development of optimal system theory. Certain of them, for example, Bellman's method of dynamic program- 
ming, can be presented in sufficiently approachable form to be useful for the circle of readers for whom the book 
is intended, 


As a whole, the book is a necessary and useful addition to the courses on automatic control theory which 
had appeared previously. 


The book contains individual inaccuracies and typographical errors, Thus, for example, an objection arises 
to the classification of automatic control systems given in § 1 of Chapter 3, It is impossible to concur with the 
estimate of the effectiveness of introducing higher derivatives in integral estimates, made by the author of Chap- 
ter 11 (é. G, Uderman), Misunderstandings can arise due to the typographical errors and to the inaccuracies in 
the diagrams of which the book is not free due, apparently, to a lack of care in its editing, It may be hoped 
that these inadequacies will be eliminated in later editions of the book, which, certainly, must shortly follow, 
since the need for a book on automatic control theory which is accessible to practical workers has yet to be 
satisfied, 


V. Yu. Kazanov 





REVIEW OF THE COLLECTION “LOGICAL INVESTIGATION" 


(Academy of Science of the SSSR, Philosophy Institute, 
Izdatelstvo AN SSSR, Moscow, 1959) 


In digital technology today, particularly in automation and in remote control, broad use is made of the 
methods of mathematical logic, which turns out to be a tremendously convenient apparatus for solving the prob- 
lems of analysis and synthesis of discrete (switching) devices, In the collection reviewed here, devoted to the 
problems of contemporary logic and its scientific applications, a significant amount of space is allotted to the 
technological applications of mathematical logic, The collection contains a number of papers, a familiarity 
with which would be of interest to engineers and to scientific workers occupied with the problems of digital auto- 
mation and remote control, In particular, these papers contain valuable results in the analysis and synthesis of 
switching circuits and discrete automata, obtained by means of mathematical logic, 


The collection consists of two parts, The first part is devoted to general questions of symbolic logic, its 
history and its particular applications to the new physics, to linguistics and to other domains, We shall here con- 
sider in detail only the second part, It contains nine papers relating to “technical logic", In addition, three 
papers are contained here which are devoted to the logical problems of the foundation of mathematics, In the 
paper of A. D. Getmanova, “On the relationship of mathematics and logic in systems of the type of the'Principia 
Mathematica',” there is presented the justification, developed by the author, of the negative solution to the prob- 
lem of the reduction of mathematics to logic, in connection with which the author deals with certain character- 
istic features of logic and, in particular, discusses the question of theoretical modes of logical systems. 





In his long paper, “On potential realizability," which warrants special discussion, A. S. Esenin-Vol'pin 
presents the “outspoken”, to use the author's terminology, point of view on realizability in mathematics, This 
point of view requires that one take into account, not only the finitude of a number, but also its pratical attain- 
ability, which might turn out to be very important for the arithmetization and algorithmization of various in- 
formation-processing schemes, 





Finally, in the paper of V. S, Chernyavskii, the concept of “shuttle” algorithms is considered, and a proof 
is given of the equivalence of shuttle and normal (in the sense of A, A. Markov) algorithms and Turing machines, 





The “technicologic” series of papers is opened by the paper of G, N, Povyarov, “Logic and automation" 
which acquaints the reader with the place and value of logic in the automation of production, There are discus- 
sed in it the fundamental problems of technological logic, the history of the application of logic to technology, 
and the basic results attained by Soviet and foreign scientists in this domain, The author mentioned that the intro- 
duction of the results of technological logic into the operations of planning organizations and into the teaching 
of thecnical disciplines in the higher schools will promote technological progress in our country and increase its 
scientific-technological power, The paper, written in popular style and with the rich bibliography removed, 
could serve as a good textbook for those who might wish to enter upon the circle of questions connected with the 
technological applications of logic, 





The two other papers of G, N. Povarov are of a different character, The first, "Mathematico-logical in- 
vestigation of the synthesis of switching circuits with one input and k outputs,” is a paper given by the author 
at an international symposium on switching theory, held at Harvard University in 1957. As is well known, this 
paper was highly praised in the foreign publications, In it are investigated, by means of mathematical logic and 
combinatorics, the taxonomy and complexity of switching (1, k)-poles, and a number of practical methods for 
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synthesizing them are given, based on Boolean algebra or its graphical interpretations, in particular, the method 
of cascades, It should be mentioned, that this method underlies the machine for synthesizing switching circuits 
which was constructed at the Laboratory for Information Transmission Systems of the AN SSSR and was success- 
fully demonstrated at the Brussels Fair of 1958 and the New York Exposition in 1959, The paper devotes particular 
attention to symmetric and other ordered schemes, 


The second paper, "On the logical synthesis of electronic computers and controlling sytems,"contains sevy- 
eral methods for synthesizing electronic switching (n,k)-poles, which are subdivided into converging and diverg- 
ing. The first are analogous to series-parallel switching circuits, the second to bridge circuits, The paper con- 
tains descriptions of the tree method, the cascade method, a strengthened cascade method and a method analogous 
to the method of universal multipoles for switching circuts, The logical basis of the synthesis is the presentation 


of some function of logical (Boolean) algebra as the superposition of elementary functions, which are realized by 
cathode followers, triodes,and pentodes, 


The paper of V, I, Shestakov, “Simulation of the operations of the propositional calculus by means of switch- 
ing circuits,” consists of two parts, The fundamental content of the first part is the establishment of the isomorph- 
ism between the propositional calculus and the algebra of pi circuits, In addition, the principles of contact and 
switching circuit realizations of the formulas of the propositional calculus are described, 





In the paper's second part, the author turns from the four-terminal realizations of two-position switches to 
the realizations of three-position switches, Then, the class of three-position (three~valued) switching circuits 
is considered, one of the constructive properties of this class being used for the modeling of the operations of D, 
A, Bochvar's three-valued propositional calculus, The paper is interesting from the point of view of disclosing 
the theoretical possibilities of the mechanization and automation of logical operations, 





The paper of N, E, Kobrinskii and B, A, Trakhtenbrot, “On the construction of general theory of logical 
nets," may properly be considered one of the central papers in the logical-technological section of the collection, 
Here, the various aspects of designing logical nets are considered, a definition of the limited-deterministic oper- 
ator is given, the physical realizability of any such operator is proven, the properties of these operators are ad- 
duced and the most general ways to synthesize them effectively are investigated, It is established, with this, that 
physical realization requires two types of elements: delay elements and logical elements, which, in the language 
of switching circuits, may be called one-cycle (nonsequential) circuits, The principles of organizing nets from 


standard elements of the types enumerated are considered, wherein the Sheffer stroke is taken as the basic logical 
element, 


The paper contains theorems on the relationships between the operator weights and the number of lag ele- 
ments, A hypothesis is asserted as to the decisive influence of the complexity of the formulas of Boolean algebra 
on the difficulty of the synthesis process and the complexity of the resulting design, For this, the compiexity of 
a formula is estimated by its length, As one of the many considerations which enter into the estimate of a logical 
net, the authors choose the minimization requirement, Although significant minimization is impossible in the 


general case, one deals in practice with simple enough operators so that a significant amount of minimization 
is attainable, 


The paper also considered the interpenetration of Boolean (logical) algebra and information theory, In 
conclusion, the paper cited the general course of the solution of the problems of synthesizing logical nets, 


The character of the paper is such that an exact presentation of all the questions illuminated by it could 


turn it into an entire book, Therefore, the presentation is naturally sketchy in certain places, which in no way 
diminishes the value of the paper, 


The paper of A, D, Kharkevich, "On switching circuits and their logical essence,® is somewhat on the side 
of the traditional presentation of the applications of logic to technology, In it, the logical structure of switthing 
circuits is considered from the engineering point of view, Discrete Q-poles are considered, the class of switching 
circuits is isolated among them and then the operations on these latter are investigated, It is possible that the 
logical principle of the paper could be distinguished more prominently, 





In his paper, “Logical synthesis of switching circuits containing impedances and contacts,” B, M, Rakov 
investigated the logical properties of two particular forms of four-terminal networks, Then, by means of Boolean 
algebra, these properties were used for the synthesis of circuits with impedances and contacts, making it possible 
to reduce significantly the number of the latter, 
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In his small work, “Analytical synthesis of sequential circuits from swiching formulas," A, N, Yurasoy 
presented the essentials of this tableless logical-algebraic method which permits one, on a step-by-step basis, 
to eliminate the possible unrealizability of a primordiai switching formula by the addition of intemediate ele- 
ments, 





The concluding paper of V, F, D’yachenko and V, G, Lazarev, “The use in telephony of logical algebra 
for the analysis and synthesis of switching circuits,” is a purely engineering illustration of the capabilities of 
mathematical logic, On the basis of a concrete attempt to design a circuit in the Laboratory for Information 
Transmission Systems of the AN SSSR, the authors draw the conclusion that the mathematical methods of analyz- 
ing and synthesizing switching circuits which are based on the laws of logical algebra give a significant econom- 
ization of time, decrease the influence of the designer's intuition and allow objective solutions to be obtained, 





Thus, the technical papers of the collection on diverse problems show the importance of the technological 
applications of logic, and contain a new estimate of the results in this direction. 


The definite many-sidedness of the second part enriches the collection, For example, it should be con- 
sidered an obvious fact that an engineer or scientific worker acquainting himself with the technical applications 
of logic can only gain if, together with the “technicological” works, he becomes familiar with, say, the paper 
of Getmanova, This allows him to overcome the almost unavoidable limitations in questions of the foundations 
of mathematics, and to increase his familiarity with such questions, 


While not having the intention of discussing the first part of the collection, the author may notdeny him- 
self the pleasure of naming the brillinat paper of £, Kol'man, "The significance of symbolic logic,” the best 
paper in the collection, 


The noteworthy fact, which it is necessary to hail, is the collaboration in one collection of philosophers 
and expounders of engineering and mathematics. 


Together with the other collections on discrete automata and mathematical machines which appeared in 
1958 -1959, such as Problems of Cybernetics, Vol. 1, Questions in the Theory of Mathematical Machines, Vol, 1, 
Proceedings of the V, A. Steklov Institute of the AN SSSR, Vol, 51, and others, the present collection provides 
valuable material to the reader familiarizing himself or interested in instructing himself in the questions cited, 


A fault which should be mentioned is the unwarranted delay in publishing the collection: from the autumn 
of 1956 to February 1959. 


Yu, L, Sagalovich 
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